Nfl30-O.P> -r75(eU>5 



NASA-CR-175665 

19850015958 


LANGLEY RESEARCH CENTER 
LIBRARY NASA 
HAMPTON. VIRGINIA 


National Aeronautics and 
Space Administration 

Scientific and Technical 
Information Program Office 

Center for AeroSpace Information 


NF00723 








NOTICE 



THIS REPRODUCTION WAS MADE FROM THE BEST AVAILABLE 
COPY OF WHICH A NUMBER OF PAGES WERE OF POOR 
REPRODUCTION QUALITY. 



' ( 






. 1 , ^ V * „ « , , ^ 

Geometric Error Analysis for Shuttle 


‘r ? ;‘ 

i .-■ • .*■ ,'r-. 

t s * * , ' ’• ' , 

(H1SA-CH-175665) ‘ G2CE2TEIC BEECH AHALYSIS H85-2t269 

FOB SHUTTLE IflAGIHG SEICTBCKE12B BXP2BIHEH1 

(Jot Propulsion lab.) 172 p EC A08/MF A0 1 , . J 

CSCI 14E Onclas 

G3/35 14802 


Shyh Jong Wang , - 
Che-Hang Charles Ih 


\ j , y':. ‘ ' 

j • - • f v 


* ~ S'- 


: ✓ • / ) 


December 15, 1984 


NASA 

National Aeronautics and 
Space Administration - - 

Jet Propulsion Laboratory - 
California Institute of Technology 
Pasadena, California 




■5 «a ■>. 

^ CP , 

tfl O’ * 

O C 




' V >' 




JPL PUBLICATION 85-2 


Geometric Error Analysis for Shuttle 

Imaging Spectrometer Experiment 

• - 

Shyh Jong Wang 
Che-Hang Charles Ih 


December 15, 1984 


IWNSA 

National Aeronautics and 
Space Administration 

Jet Propulsion Laboratory 
California Institute of Technology 
Pasadena, California 


The research described in this publication was earned out by the Jet Propulsion 
Laboratory, California Institute of Technology, under a contract with the National 
Aeronautics and Space" Administration 

Reference herein to any specific commercial product, process, or service by trade 
name, trademark, manufacturer, or otherwise, does not constitute or imply its 
endorsement by the United States Government or the Jet Propulsion Laboratory, 
California Institute of Technology 


Acknowledgements 


The authors wish to express their gratitude to Y. B. Lin and J. M. Cameron 
for their contributions at the early stage of this work, and to M. Herring, 

P. N. Kupferman, and J. B. Wellman for their support. 


ill 







ABSTRACT 


The demand of more powerful tools for remote sensing and management of 
earth resources has been steadily increasing over the last decade. Uith the 
recent advancement of area array detectors, high resolution multichannel 
imaging spectrometers can be realistically constructed. 

This report documents the error analysis study for the Shuttle Imaging 
Spectrometer Experiment system for the purpose of providing Information for 
design, tradeoff, and performance prediction. 

Error sources including the Shuttle attitude determination and control 
system, instrument pointing and misalignment, disturbances, ephemerls, earth 
rotation, etc., have been investigated. Geometric error mapping functions were 
developed, characterized, and Illustrated extensively with tables and charts. 
Selected ground patterns and the corresponding image distortions have been 
generated for direct visual inspection of how the various error sources affect 
the appearance of the ground object images. 
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I. INTRODUCTION 


A. BACKGROUND 

Earth resource management and utilization have experienced great-success 
over the last decade through the Landsat programs. In more recent years, both 

NASA and user communities have envisioned the need for development of better 

•• 

and more powerful instrument's for surveying and managing earth resources. The 
Landsat D's new sensor. Thematic Mapper (TM) , the preposed utilization of 
Tracking and Data Relay Satellite System (TDRSS) , and a more ? j.'anced ground 
system represent an advancement in the earth resource satellite development [l]. 

The Thematic Mapper of the Landsat D (launched in 1982) has seven 
spectral bands, two more than those of the Multispectral Scanner associated 
with the °r.rlier Landsats. However, study shows that the reflectance spectrum 
of earth surface materials contains a significant amount of information which 
can only be identified with spectral resolution much finer than those of the 
Thematic Mapper [2l. With the advancement of area array detectors, a push 
broom imaging spectrometer can be realistically constructed for simultaneous 
Imaging and registration of hundreds of spectral bands. For the case of 
Shuttle Imaging Spectrometer Experiment, 12S spectral channels have been 
proposed to cover the spectral range of 0.4 to 1.0 pm for VN1R (visible and 
near infrared) and 1.0 to 2.5 pm for SWIR (short-wavelength infrared) with 
instantaneous field of view of 30m. Table 1 shows the required sensor perfor- 
mance [2]. 

The purpose of this report is to document the error analysis stud/ for 
the imaging spectrometer experiment. Error analysis is an important aspect of 
the overall remote sensing system, since errors from many sources, including 
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the spectrometer Itself, the spacecraft that carries the instrument, knowledge 
limitations on the true spacecraft attitude and locations, earth rotation, 
curvature, and terrain variations, etc., will all contribute to image distortion, 
shift, rotation, and misregistration. Error correction or compensation are 


necessary and are an Integral part of the image processing. This work covers 
the analysis of fundamental and geometric errors and error sensitivities, the 
development of geometric mapping functions, and the computation of ground 

Table 1. Sensor Performance Requirements 

Parameter 

Value 

Comments 

Spectral Coverage 

0.4 to 2.5 pm 

Although the entire spectrum 
is probably not required for 
any one discipline, in the 
aggregate of all remote 
sensing disciplines, the 
entire region is required 

Spectral Sampling Interval 

VNIR (0.4 to 1.0 urn) 

0.01 pm or better 


SWIR (1.0 to 2.5 pm) 

0.02 pm or better 


Instantaneous Field of View 

30 m 

Adequate for most research 
topics 

Swath Width 

at least 10-12 km 

This is adequate for 
research if pointing capa- 
bility is provided to assure 
target access 

Pointing Mirror Range 

Along track 

at least + 45 deg 

Essential for atmospheric 
and BRDF (Bidirectional 

Cjcss track 

not more than + 

Reflectance Distribution 


25 deg 

Function) Studies 

Radiometric Performance (NEdR) 


VNIR 

< 0.5Z 


SWIR 

< 1.0Z 
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pattern distortions. The results are translated into many tables, plots, and 
patterns for visual apprehension. It is believed that the results reported 
here are important for design, trade-off, and performance prediction. 

B. APPROACH 

This study has been performed in four progressive stages as shown in 
block diagram form in Fig. 1. In Stage I, the error sources were identified, 
dynamic disturbances and the Space Shuttle error dynamics were modeled, and 
the error power spectral densities for two in-orbit configurations w£~re 
developed. Stage II of this study concentrated on the development of geometric 
error mapping functions and error sensitivity analysis. Geometric errors due 
to ephemerls uncertainties, attitude deviations, earth rotation, etc., were 
studied. In Stage III, ground pattern Image distortions caused by various 
error effects and forward and side looking angular offsets, altitude change, 
and effect of earth rotation were generated. Stage IV consists of the analysis 
of the imaging spectrometer instrument errors. These errors include optical 
jitter, nonlinearities, processing errors, and repeatability. The study of 
Stages I, II, and III has been completed and the results are Included in this 
report. The study called for by Stage IV has not been planned. It is emphasized 
here that the imaging spectrometer error model development is an important step 
for the overall performance prediction and design of the imaging spectrometer 
system. 

Major findings of this work are summarized in the following subsection. 

In Section II, the orbital and imaging spectrometer configurations are described. 
The attitude dynamics and error power spectral densities are documented in 
Section III; and parametric analysis of geometric errors are treated in Section 
IV. Section V deals with the ground pattern image distortions. Conclusions 


are summarized in Section VI. 
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Figure 1. IS Error Analysis System Block Diagram 1 
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SUMMARY OF MAJOR FINDINGS 


The following Is a summary of major findings. The details of these are 
treated in the sections to follow. 

1. The results show that the IS Experiment with Image pickup period of 

20 seconds at a time is feasible with the shuttle properly phased 

inside the control deadband. The error PSD (power spectral density) 

characteristic reveals that the system resonates at very low frequen- 

• ~ 

•*5 . 2 

cles (in the 10 to 10 A Hz region) . Excitations at these frequen- 
cies must be avoided through design precautions. The analysis also 
showed that errors below 0.01 Hz are dominated by the shuttle 
dynamics reacting to disturbances, whereas those above 0.01 Hz are 
dominated by the shuttle inertial measurement system uncertainties 
and its Inherent noise. These high frequency errors limit pointing 
performance and result in a one-sigma ground track error of 54.8 
meters per axis. One-sigma rate errors are shown to be less than 

4 meters per second per axis in the frequency range of 10 ^ to 
_2 

4 x 10 Hz. The image smear is not significant because of the short 
millisecond-level line time. 

2. The effects of earth curvature are very small for the application 
here (see Fig. 42). 

3. Altitude uncertainties cause only moderate geometric errors. The 
worst la geometric errors are 11.71 m in position and 0.133 m/sec 
in rate with STDN and large unmodeled perturbations at 200 km orbit. 
The performance Improves with TDRSS. For the 300 km orbit and with 
small unmodeled perturbations, the la geometric errors will reduce to 
0.22m and 0.0032 m/sec (see Table 11). 
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4. The effects of other navigation errors are significantly greater. 

The lo downtrack errors range from 203m (300 km orbit) to 8128m 
(200 km orbit) ; and those for the crosstrack are 152m to 508m 
(see Table 12) . 

5. The effects of roll and pitch attitude errors are relatively large 

compared with, for instance, those caused by yaw errors and altitude 

« - 

uncertainties. The error sensitivity is 1.94 m/arc sec or approxi- 
mately 7000 m/degree. The yaw sensitivity is 0 for 0° view angle 
and 0.028 m/arc sec for the maximum view angle of jjO.825 0 (see 
Table 13). 

6. The error sensitivities of attitude errors increase significantly 
for large attitude offsets. For 20° side looking, the sensitivity 
increases to 2.28 m/arc sec for roll errors and to 0.75 m/arc sec 
for yaw errors and the pitch error sensitivity is almost unaffected. 
For the 45° forward looking case, the sensitivities for the pitch, 
roll, and yaw errors increase to 2.83, 4.36, and 1.97 m/arc sec, 
respectively (see Table 13) . 

7. The performance of the Shuttle Imaging Spectrometer is limited by 
the Shuttle IMU (Inertial Measuring Unit) accuracy, instrument mis- 
alignment, etc.. Shuttle RCS (Reaction Control Subsystem) deadband, 
etc. unless some means of error reduction are employed. For Instance 
ground control points may be used to reduce navigation prediction 
errors; and precision point mounts, such as AGS (ASPS* Gimbal System) 
and IPS (Instrument Pointing System) , may be used to reduce the 

♦Annular Suspension Pointing System 
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attitude errors. 

The single axis geometric errors due to the combined Shuttle/IS mis- 
alignment, for instance, for normal nadir pointing, 20° side looking, 
and 45° forward looking are 169m, 198m, and 379m, respectively 
(refer to Table 10 for detailed breakdowns) . 

8. Earth rotation causes shifts of images toward the direction of 

rotation. The magnitude of this shift depends on the latitude of •• 
the object. For instance, at the equator the object moves approxi- 
mately 462m in one second (for 400 Vm orbit) , and at 60° latitude 
it moves only 231m in one second. 


; i " 

i < 

i 

« 

1 




X 


* 



II. ORBITAL AND IMAGING SPECTROMETER CONFIGURATIONS 

A. SHUTTLE NOMINAL FLIGHT CONFIGURATIONS 

Two shuttle in-orbit configurations have been considered, the Payload- 
Bay Nadir and the Nose-Down Nadir, as illustrated in Figures 2 and 3, respec- 
tively. It will be shown later that the Nose-Down Configuration is gravity 
gradient stabilized and the Payload-Bay Nadir Configuration is not. However, *- 
the Payload-Bay Nadir Configuration offers simpler instrument mounting and less 
aerodynamic drag. Besides, for certain experiments that require large forward 
looking angles, the Nose-Down Configuration will be unsuitable. 

A circular orbit of 400 km altitude has been selected for this analysis. 

* 

This selection is consistent with the SIS-B parameters. 

There are a number of possible instrument mounting options that will 
affect the pointing and the geometric errors of the instrument. These options 
include : 

a) Direct mount 

b) AGS (ASPS Glmbal System) mount 

c) IPS (Instrument Pointing System) mount 

Since both the AGS and IPS systems are capable of providing precision payload 
pointing and measurements, the system performance will be improved at the expense 
of significant extra cost. The direct mount approach is the least complex and 
most cost effective* provided that the errors are within the tolerance. In this 
report the direct mount approach is considered. Furthermore, the shuttle IMU 

^During the period when this part of the work was performed, SIS-B (Shuttle 
Imaging Spectrometer-Configuration B) was considered. However, the methods 
used here are mostly generic, and hence, can be applied to systems of similar 
configurations . 
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Earth 

Figure 3. Nose-Down Nadir Pointing Configuration (B) 


9 




T i 


'a 





(Inertial Measurement Unit) and the shuttle state estimator unit are used for 
obtaining attitude and rate Information, without additional Instrumentation. 

Other options may be studied In the future if necessary. 

B. IMAGING SPECTROMETER SYSTEM DESCRIPTION 

In order to correlate the analysis reported here to the IS applications, 
it is desirable to understand the basic operating principle of the IS instru- 
ment. Figure 4 shows a sketch of the basic elements of the IS, except the 
memory banks, registers, and the processors. Some of the IS parameters that 
are relevant to the geometric error analysis, along with definitions of the IS 
terminologies, are also listed in Fig, 4. It is understood that the sketch and 
the parameter values used are for illustrative purposes only since the para- 
meters may change as the development of the IS system is finalized. 

Referring to Fig. 4, as the reflection from the ground objects passes j - 

through a narrow slit, it strikes the incident surface of the prism. The prism I - 

separates the incident light into spectral images projected onto an area array 
of light detectors. The area array consists of 384 linear arrays, corresponding 
to 384 spectral channels. Each linear array consists of 384 detector elements. 

Each detector element corresponds to an image area of 30m x 30m on the ground. 

Hence, each linear array corresponds to an image of a specific spectral channel 
of the same ground area of 30m (along-track) x 11520m (cross-track). This 11.52 
km cross-track measure, referred to as the swath width, defines the IS coverage 
for each orbit pass. The 30m x 30m area is referred to as the pixel (picture 
element) which defines the resolution of the image, i.e., within this element 
no features can be resolved. The IS is operated based on a push broom principle. 

That is, a specific detector element on the instrument (moving with the space- 
craft) collects photons from the 30m x 30m moving window for a specific period 
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30 m x 30 m 

DETECTOR 

ELEMENT 


AREA ARRAY 
(CONSISTS OF 36 
MODULES OF 
64 x 64 DETECTORS) 


! 

i 


KEY PARAMETERS 
• LINE TIME = 4 m 


• SPECTRUM: 

VISIBLE REGION 
SHORT WAVELENGTH IR 


0.4 - 1.0 
1.0 - 2.5 


• FIELD OF VIEW: ± 0.825° 


ftm 
fi m 


• SWATH WIDTH: 11.52 km 

• POINTING ANGLES (RANGE): 

± 45° (PITCH) 

± 20P (ROLL) 

• ORBIT: 400 Ian CIRCULAR 


DEFINITIONS 

GIFOV (GROUND INSTANTANEOUS FIELD 
OF VIEW): PROJECTION ON THE GROUND 
OF EACH SQUARE DETECTOR 

LINE TIME: TIME TO MOVE ALONG THE 
GROUND A DISTANCE OF 1 GIFOV 

DN: DATA NUMBER, THE -BRIGHTNESS" OF 
OF THE ASSOCIATED PIXEL 

PIXEL: PICTURE ELEMENT (30 m x 30 n. 

HERE, IT DEFINES THE RESOLUTION OF 
THE IMAGE) 

SWIR: SHORT WAVE LENGTH IR 
VIS: VISIBLE WAVE LENGTH 


Figure 4. Imaging Spectrometer Basic Operating Principles 

i 
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of time called the line time. The line time in this case is the time required 
to move 30m along-track, which is determined by the orbit. For a 400 km circular 
orbit, it ia about 4 ms. At the end of each line time, the total number of 
photons collected by each detector is recorded and processed and the detector/ 
register is reset and a new 4 ms cycle is repeated. For digital processing, 
the amount of light collected by a collector is assigned a number called DN 
(data number), which is proportional to the number of photons accumulated. 
Therefore, the processor has to record and process 334 x 384, or approximately 
1.475 x 10 5 DN's every 4 ms. 

As the shuttle flies over an area, banks after banks of DN's are 
collected. By spacing the banks of DN's 30m apart, the features of the ground 
image emerge. The ratios of the DN's of various channels for the same area are 
of special interest, as these ratios are closely correlated to geological and 
ecological states of the earth including mineral deposits, forestry, crops, 



Figure 5. Shuttle Imaging Spectrometer 
Optical System Configuration 
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disease and insect infestations, land and soil erosion, precipitation in 
the form of snow and ice, air and water quality, etc. 

The Imaging Spectrometer optical system configuration and the arrange- ' 

i 

l 

ment of lenses, mirrors, slit, prisms, focal plane detector, etc. are , 

1 

! 

shown in Fig. 5 [2]. Additional information on the instrument design and require- 


ments can be found in Refs. 2 and 3 


1 


% 










III. ERROR SOURCES AND SHUTTLE ATTITUDE DYNAMICS 
A. ERROR SOURCES 

There are two types of errors that are important to the Imaging 
spectrometer experiment — the direct errors and the derived errors. The 
direct errors are those pertaining to the spacecraft and instrument pointing, 
ephemeris, and instrument errors as shown in Fig. 6. The derived errors are 
the geometric errors and pattern distortions of ground objects which result 
from the direct errors and the effects of earth rotation, curvature, 
oblateness, and local vertical uncertainties. Fig. 1 shows relationships of 
the error sources and the system dynamics. 

Attitude deviations, rate errors and structural vibration 
Measurement noise and drift 

Misalignment of shuttle imu and attitude reference frame 
Misalignment between: shuttle and IS instrument references 
Ephemeris prediction errors 
Earth rotation, curvature, and oblateness 

IS INSTRUMENT ERRORS INCLUDING 
OPTICAL JITTER 
NONLINEARITIES 
PROCESSING 
REPEATABILITY 


Figure 6. Error Sources 
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In this section, the steady state analysis is performed In the 
frequency domain. The PSD (power spectral density) for the pointing errors 
and the rate errors were obtained by considering the dynamics cf the Space 
Shuttle Orbiter, the DOJ (Inertial Measurement Unit) , the attitude control 
state estimator, the measurement noise, the misalignment of reference frames, 
and the disturbances including gravity gradient, gyroscopic torques, and aero- 
dynamic drag torques. 

B. SHUTTLE MASS PROPERTY 

The shuttle mass properties employed here were obtained f.om the 
Shuttle Operational Data Book [4] for 0V-102/STS-3. The mass and the c.m. 
are, respectively* 

Shuttle mass: 102,153.73 kg (224,738.21 lbs) 

Shuttle c.m.: X » 1105.5", Y =■ 0 , Z ** 374.3" 

o o o 

where the coordinates X q , Y q , and Z q are the Orbiter Coordinates [5] defined in 

2 

Figure 7. The moment of inertia matrix (kg-m ) is, in the shuttle body frame 
(refer to Fig. 8), 


1.36 x 10 6 -4.69 x 10 3 -3.49 x 10 5 


I„ - -4.69 x 10 J 


1.00 x 10 7 -3.32 x 10 3 


-3.49 x 10 5 -3.32 x 10 3 


1.05 x 10' 


1.36 x 10 u 


-3.49 x 10 J 


1.00 x 10 


-3.49 x 10 J 


1.05 x 10' 


The magnitude of the off-diagonal terms of the inertia matrix I 

D 

suggests strong couplings exist especially between the X^- end the Z^-axls. To 


simplify the dynamic equations, it is convenient to cousider principal-axis 
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ORBITER COORDINATES 

TYPE: Rotating, Orbiter referenced 

ORIGIN: Approximately 200 inches (5.1m) ahead of the 

nose and approximately 400 inches (10.2m) 
below the centerline of the cargo bay 

ORIENTATION AND LABELING: 

The X-axis is parallel to the centerline of 
the cargo bay, negative in the direction of 
launch 

The Z-axis is positive upward in landing attitude 

The Y-axis completes the right-handed system 
The standard subscript is 0 

Figure 7- Orbiter Coordinate System 




pointing Instead of body-axis pointing, e.g., pointing -Z rather than -Z_ as 

P B 

illustrated in Figure 8. 



Earth 

Surface 


Earth 

Center 

(Xg, Yg, Zg) — Shuttle Body Coordinates 
(Xp, Yp, Zp) — Shuttle Principal Axes 

Figure 8, Shuttle Body and Principal Coordinates 
(Payload-Bay Nadir Pointing Shown) 

The orientations of the principal axes can be determined by rotating 
Che tody frame an angle a about the Yg-axis, i.e., let B: Xg -*• X^, then 


B - 


cosa 

0 

since 


0 

1 

0 


-sina ~i 
0 

cosa 


and 




(2) - 


( 3 ) 


It can be shown that 
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1 ► - 1 
a m 2 Can 


r 3ZZ " I BXX 


provided that I m - Ig YZ - 0. 

FOr T BXX = 1,36 x 1()6 kg-” 2 . \zz a 1 * 05 x 10 ? kg-m 2 , and I BXZ = 

5 2 o 

3.49 x 10 kg-m , the angle a = -2.18°. Therefore, in order to do principal- 
axis nadir pointing, the shuttle has to rotate 2.18° about the orbit normal (see 
Figure 8) . The moments of inertia about the principal axes are 

I p - diag (1.38 x 10 6 , 1.00 x 10 7 , 1.05 x 10 7 ), kg-m 2 . 


ASSESSMENT OF DISTURBANCES 


The shuttle motion is characterized in part by the environmental 
disturbances, the major sources of which are the aerodynamic drag, gravity 
gradient, gyroscopic effect, solar radiation, and on-board causes such as 
astronaut activities, equipment vibrations, and venting. On-board activities 
may be partially eliminated or reduced through mission planning and their 
effects will be assessed in the future. In this subsection, the gyroscopic 
torques, the gravity gradient torques* and the aerodynamic drag torques are 
estimated. The solar pressure is at least one order of magnitude less than 
the aerodynamic drag forces and*hence, it is not included in this analysis. 

C.l The Aerodynamic Drag Torques 

To estimate the aerodynamic drag torques, an approach referred to 
as the three-plate model [6] has been used. Referring to Figure 9, let n^ n 2 , 
and a 3 be the unit normal vectors for the equivalent plates 1, 2, and 3. 

Vhere n^, n^, and n^ are in the direction of Xg~, Yg- and Zg-axis, respectively. 
To express n^ in body frame, 





( 5 ) 



Figure 9. The Three-Plate Aerodynamic Model in Body Frame ; 

I 

Let A^ f A 2 , and A-j be the corresponding areas on which the aerodynamic 
pressure applies, and C pl , C p2 , and Cp 3 the three corresponding centers of 
pressure, respectively. The model assumes that the aerodynamic force applies 
to each area in the direction opposite to the vehicle velocity vector and 
assumes no shading among the plates. It is further assumed that the drag 
coefficient C Q is constant and the lift coefficient is zero. 

_ X 

Let v g - ( v Bl * "jj 2 » v b 3 * be thc velocity vector in the 

body frame, with magnitude v. The force and torque applied on A^ are 

? Bi “ ~ ( T V ' /2 V 1 ( V V I “B <6) 

i 

? Bi " F Bi X ? 3i ' " ( I S P I (F i’ V I *Bi X “b (7) 

where P is the atmospheric density at the orbital altitude, u • v /v the 

B B 

unit velocity vector, and r gi the vectors of the center of pressure of A^ relative 
to the vehicle center of mass. 
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( 8 ) 


The total aerodynamic forces and torques on the vehicle are 
3 

? B “ £ “ ( I C D P v2 A i> I ( V V I “b 
i-1 

3 

? B “ £ " 4 C D P v2 V I ( V V l 7 Bi x (9 > 

i-1 

The value of C D depends on the shape of the vehicle. In Ref. 7, 
the values of 2.5 to 3.0 were suggested. The value of - 2.0 is used here 
as it was used in Ref. 6 for Space Shuttle Simulation. 

The atmospheric mass density p can be found in a JPL internal memo- 
randum. For 1935 mission time, the peak density is expected to occur in April 
for the 400 km orbit; the densities are 

2.64 x 10 * 2 kg/m^ (predicted) 

and 

3.77 x 10 * 2 kg/m 2 (97.7 percentile) 

The orbital velocity v for a 400 km circular orbit is computed as 

2 7 2 

v - 7669.60 m/sec or v - 5.882 x 10 (m/sec) . 

The plate areas, according to the attachment (SSFS On-Orbit Aero 

Data, 7/24/80) to Ref. 6, with Cargo Bay doors open, are 

A x - 119.45 m 2 

A 2 - 229.92 m 2 (10) 

2 

A^ - 454.46 m 

However, a different set of values are given in Ref. 8. 

A x - 64.1 m 2 

A 2 - 212.7 n 2 (11) 

A^ - 367.0 
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f 


The values in (10) were used In this work. 

The position vectors r Q ^ nay be obtained using data given in Ref. 4 
and Ref. 6. The values are. 


B1 


3.797 

0 

- .231 


L 


m 


B2 


T * 

r B3 


.927 

0 

- .742 

1.166 

0 

-1.074 


m 


The torques T^ for the predicted density, are 


T B i - - i.855 x it) |u Bl | 


" r B13 U B2 + r B12 U B3 


+ r Bll u B3 


C B13 U B1 


" r Bl2 U Bl + r Bll u B2 


( 12 ) 


U 


“ “ 3.570 x 10 -2 |u_-| 


B2 


B2' 


" r B23 U B2 + r B22 U B3 
r B23 U Bl " r B21 u B3 


' r B22 U Bl + r B21 U B2 


(13) 


T 33 “ ‘ 7 * 057 X 10 " 2 | U B3' 


" r B33 U B2 + r B32 U B3 
r 333 U 3l " r B31 U B3 
• r B32 U B1 + r B31 U B2 
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In Eq. (13), Che torques are known once the unit velocity vector u^ Is 
specified! u" B varies with the pointing configuration and spacecraft attitude. 

C. 1.1. Drag Torques for Configuration A 

Referring to Fig. 2, since this is ”-Z " pointing, for nominal 
attitude, the vehicle moves in the Xp direction. Let u p be the unit velocity vector 
in the P-frame and let be the rotation matrix due to a, then 


U B " A a U P 


cosa 

0 

sina 

0 

1 

0 

-sina 

0 

cosa 

cosa 



0 



-sina 




Since a(* - 2.18 - - .038 rad) is small, 
l' 

-B ■ ° 

Using Eqs. (12) ,(13), and (15), and 
3 


■E 


A a T Bi 


1 

0 

0 


i-1 


the tcrqucj in the P-Frame, for nominal attitude, are 


(14) 


(15) 



(16) 



9.84 x 10 
0 


-3 


N-m 


(17) 


For small attitude errors, $, 0, and ty, from the rotating orbital frame. 


1 

0 

C 


(18) 
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where 


and 


1 

-* 

e 


T 

A A 


* 

1 


1 

0 

u°j 


-0 


1 J 


( 20 ) 


Using Eqs. (12), (13), and (20), the predicted aerodynamic torques become, by 
retaining only the first order terms, in N-m, 


•£ 

i-i 


a Bi 


( 21 ) 


4.29 x 10" 3 * 


4.29 x 10" J + 7.04 x 10 -2 (9-a) + 7.58 x 10 -,t |e-aJ + 2.65 x 10 _i£ |i|<| 


- 2 , 


.- 2 . 


7.04 x 10 -2 $ 


C. 1.2. Drag Torques for Configuration B 

Referring to Figure 3, under this configuration, is in the Nadir 
direction and Y_ is in the direction of motion for nominal attitude. In this case. 

To! 


U B " A a U P 


( 22 ) 


( 23 ) 


and from Eqs. (12), (13), and (23), 
3 

vL v.i 

i-l 



vs 


-3.96 x 10 
0 

-4.53 x 10 


-2 1 


-2 


N-m 


(24) 


In the presence of small attitude errors, $, 9, and $, the aero- 
dynamic drag torques are 

-2.45 x 10~ 2 - 7.58 x 10 _2 |$|-4.29 x 10 _3 j $1+3.31 x 10 _2 a 


-3.31 x 10 -2 $ + 2.65 x 10 " 2 $ 


(25) 


|_— (3 . 3 1 x 10 -2 + 8.23 x 10~ 2 1 <j>|+7.04 x 10~ 2 |$|+2.65 x 10 _2 a)J 


(25) 


C.1.3. Estimation of Disturbance PSD 

Before the PSD's are estimated, the uncertainty part of the distur- 
bance torques has to be determined. The torques of Eqs. (21) and (25) consist 
of static parts and the dynamic parts. The dynamic parts are functions of the 
attitude errors $, 9, and $. The attitude errors are assumed to be random and time- 
varying with standard deviation of 1° (.01745 rad.) per axis. Therefore, the 

estimated values of the random disturbance torques are, for Configuration A 

,-5 T 


PA 


7.48 x 10 
1.86 x 10 
1.23 x 10 


-3 

-3 


N-m 


(26) 


■ ^ i 


and that for Configuration B are. 


PB 


1.33 x 10 


7.40 x 10 
1.89 x 10 


-3 

-4 

-3 


N-m 


(27) 


and the corresponding PSD's are obtained by the following approximation with 
the correlation time of T = 180 seconds. 


PA 


2T (o 2 ) 

PA 


24 


T 



/ 


2.01 x 10 


-6 


1.25 x 10 


-3 


5.44 x 10 


-4 


(N-m) - sec 


(28) 


PB 


2T (<£ ) 

PB 


6.32 x 10 


-4 


1.97 x 10 


1.29 x 10 


-3 


(N-m) - sec 


(29) 


C.2 The Gravity Gradient Torques and Gyroscopic Torques 

The gravity gradient torques and the gyroscopic torques can be 
estimated using the following equations, respectively ' 


Where u^ and &> o p are the unit earth vector and the spin velocity vector, 
respectively, in, principal frame, and u^ is the skew symmetric matrix of the 


vector u . For Configuration A, 


“up 


- A 


0 

0 

LU 


(32a) 


“oP " 


Uj 


(32b) 


and for Configuration B, 
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and 


T gP ” ^o “rp x p “rp 


t grp " M 0 P x np * “op *p “op 


(30) 

(31) 



(33a) 


U RP 


A 


oP 



(33b) 


the corresponding torques, to include only the first order effect, are, for 


Configuration A 



T 

CRP 



__ m 


“ m 

^PZZ “ 1 PYY^ 

(I PZZ " i pxx )9 

s 

1.92 <j> 
33.11 9 

0 

m « 


0 

“^ X PZZ “ I PYY^ 
G 

S3 

-.64 0 
0 

r 

M 

*0 

a 

i 

M 

. -o- 


11.04 <|» 


and for Configuration B 


(34a) 


(34b) 




0 


0 

(I PZZ ~ I PXX^ 9 

S3 

- 35.03 0 

- (I PXX " I PYY ) ^ - 


43.71 ip 

'^PZZ “ I PYY^° 


-.64 ip 

(i pzz ~ X pxx^ e 

- 

11.68 0 

0 


0 


(35a) 


(35b) 


Since the gravity gradient torques and gyroscopic torques are propor- 
tional to the attitude errors, they may be included in the equations of motion 
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as pare of the vehicle dynamics rather than disturbances to the plant. 

D. THE MISALIGNMENT AND MEASUREMENT NOISE 

In this subsection, the errors that will contribute to pointing 

uncertainties are estimated. These error sources Include the misalignment 

errors of the IMU and the IS Instrument, and the sensor noise. 

D.l The Misalignment Errors 

Let b„_ be the misalignment between the IMU and the shuttle reference 
MS 

frame; let b__ be the misalignment between the shuttle and the imaging spectro- 
meter reference frame. 

Based on the space shuttle performance requirement [5], the IMU 3o 
misalignment uncertainty is + . 133°/axis, hence. 


a - .044 /axis “ 159.6 arc-sec/axis 
bMS 


(36a) 


However, based on a JPL internal memorandum, the shuttle flight test performance 
is much better. 


a bMS = 82 arc “ sec / axis 


(36b) 


The performance data Eq. (36b) was used in this report. 

The misalignment data for the IS instrument is not directly available. 

Based on Ref. 9, the estimated LANDSAT U initial alignment bias between the sensor 

optical axis and the vehicle pointing vector was + 200 arc-sec. The alignment 

bias can be measured before launch and can be removed from the image data. The 

variation part that cannot be removed without ground control points was estimated 

as + 30 arc-sec. This latte > number was used in this work, 

a, - 30 arc-s- — *s (27) 

b SI 

Therefore, the total misalignment error becomes 

2.2 . 1/2 
a. “ (°i, + a b ' 

b o °MS b SI 
- 87.32 arc-sec. 
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D.2 The IHU Model 

The sensor dynamics for Che shuttle IMU was modeled approximately 
by the first order low pass filter as shown in Figure 10. 



Figure 10. Simplified IMO Model 

Where t is the time constant and v the white gaussian noise. The time constant 
was estimated from a DRIRU II Model, 

T " " T4x “ * 023 sec (39) 

s 

That is, this measurement is assumed to have a 7 Hz bandwidth. 

To determine v, from Ref. 5, the 3a IMU readout error is + .G73 c /axis. 

Assume that the measurement noise is the readout error, then 

a m .0243°/axis = 87.6 arc-sec/axis (40) 

v 

However, the actual performance of Che shuttle IMU was much better; it has a 
la gyro resolution error ( CT ££ 5 q) of 20 arc-sec/axis. If we assume the random 
noise has the same amplitude as tnat of the resolution noise, then 

°» ‘ * °RES0 (41 > 


The corresponding measurement noise PSD, R, is estimated as, 

R : 2 (k) °v 2 

» 8.55 x 10 (rad) 2 -sec 
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E. SHUTTLE DYNAMICS AND INSTRUMENT POINTING ERRORS 

E.l The Equations of Moclon 

Consider Configuration A. For simplicity, the subscript P for 

principal axes is dropped from all equations. Again, let <?, 9, aad be the 

— T 

small roll, pitch, and yaw angles, respectively. Let w - u. ; , u z ) be the 

angular velocity vector of the shuttle, then for small attitude errors, 

«> x * i + <|> u Q 

u Y - 9 + « o (43) 

u» z * i> - t io Q 


and 


• • 

oj v “ $ + V u 
X ( 

Wy “ 9 


(44) 


6» 

Z o 

With this simplified relation between the inertial rates and the attitude error 
rates, one can 3how that the equations of motion may be summarized as follows, 
accounting only for first order effects! 

2. T 


V + (1 x - X Y + V “o* + 4 C1 Y - V “o * 

V + 3 (I X " I Z ) U o 2 ® 

- d v - i Y + i 2 ) - d x - i Y ) u 0 \ 


T . + T 
d c 


(45) 


Z T 'X 

The left-hand side of Eq. (45) accounts for the body dynamics ana the gravity 
gradient and the gyroscopic torques; and on the right-hand side of Eq. (45) T, 


and T are the disturbance torques and the control torques, respectively, 
c 

E.2 The Space Shuttle State Estimator 

The Shuttle State Estimator consists of two parallel Kalman-type 
filters, one for acceleration estimation and one for rate estimation [lO], 

The attituae estimate is the extrapolation of the measured attitude and the rate 
estimate may be approximated by using a second order filter with parameters 
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determined by the filter gains. With the current baseline filter data, the 
equivalent damping coefficient and the corner frequency for the rate filter are 
.8 and .04 Hz (for vernier rate filter), respectively £ll]. 


E.3 The System Pointing Errors 

E.3.1. The Fitch Loop (Configuration A) 

From Eq. (45) and ignoring the initial conditions, one has the 
following output function. 


6(S) 


v 2 + 3 “o 2 ' vv 


(I dV * V 


(46) 


Figure 11 8 hows the open-loop block diagram for the dynamics of the instrument 
pointing error excited by the random disturbances. Included in the diagram are 
the dynamics of the vehicle, the IMU, and the rate filter. Measurement error 
and the misalignment error are also included. 

E. 3.7. The Roll and Yaw Loops (Configuration A) 

The output equations for the coupled roll and yaw dynamics are, from 

Eq. (45) 


%(S)- 

1 

*(S) 

D(S) 


l Z S + “o ( ' l Y~ 1 ^} "“o^W^^ 

x : 


“o ( VW S 


I„S 2 4- 4w q 2 (W 


(47) 


where 


and 


D(S) - I x i/*o o 2 (48) 


T x - T cX + T dX 
T Z “ T cZ + T dZ 


(49a) 

(49b) 


Figure 12 shows the block diagram for the roll and yaw instrument error dynamics. 
The block diagrams for the pitch and the roll-yav loops are similar except for 



the coupling terms for the roll and yaw axes. 


Random 

Disturbance 


Measurement Misalignment 
Error Error 



Figure 11. Pitch-Axis Block Diagram (Configuration A) 


Simplified 
Rate Filter 
Model 


E.4 The Instrument Pointing Error PSD 

E.4.1 The Pitch Error PSD (Configuration A) 

The PSD' s for the pitch-axis instrument pointing error and the rate 

error are, 

? e («) - F e *(u) F 0 ( u)Q Y + R + °bo 6(w) (50) 


P-(u) - (Fg *(<!)) F q (w)Q y + R) H r (oj) H>) (51) 

where Q v •• Q_ Y and Q „ is the second component of Eq. (28), and where R, and 
Y ^PA *PA 

a, are defined in preceding sections, and 5(u) is the Dirac delta function; 

DO 


and where. 


F 9 (u) - G 9 (m) H(u) 


G e <c) -J- = ” 2 '1 

Lv + 3 "o l vv J 
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Random 

Disturbance 


dx 




cX 
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Dynamics 
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D(S) 



Dynamic Coupling 


D(S) 
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Figure 12 


Roll and Yaw Axes Block Diagram 
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o 

Misalignment 

Error 


(Configuration A) 








(52c) 


H (u) ” C T s + i ^ s _j U 


F 0 *(u) - F 0 (-u) 


r - 2 •> 

S a) 

H(») * — " — 

r S + 2{ u S + m S=ju 

r r r _j J 


The output power spectral density for the system is closely related 

to the frequency response of the system,vhich characterizes the steady state 

dynamics of the system, and,hence, it is meaningful only if the system is 

3table. Unfortunately, in Eq. (52b), 1^ < which implies that the system is 

unstable. The destabilizing term comes from the gravity gradient because the 

nadir pointing axis is not the axis of minimum inertia. To proceed, one has to 

consider the gravity gradient as external disturbance rather than a part of the 

2 

dynamics. Figure 13 shows the instrument pointing error PSD in (rad) -sec as 

2 

a function of frequency in Hz. Figure 14 shows the rate error PSD in (rad/sec) -sec 
as a function of frequency. 

E.4.2. The Roll and Yaw Error PSD (Configuration A) 

The instrument pointing error PSD's for the roll and yaw axes are, 

V“> - F $<“> y w > Ox + % ^ y < u > 

+ R + o£ o 6(w) (53a) 

y«) - F |< u > y«> + % <"> 

+ R + a. 6 (to) (53b) 

DO 

and the instrument pointing rate error PSD's for the roll and yaw axes are, 

?•(«) - [F*(ui) F^(u) + F*^ («) F^(u) Q z + R] H* (w) H.(u) (54a) 




[F*( U ) y«> Q z + F*^ (co) F^ (u) Q x + R] H*(u.) H>) (54b) 


*$<"»> “ 


where 


F,(u) - Gi(w) H(oi> 


G . (id) 


I S 2 + u 2 (I -I ) 

Z O V X 

D(S) 


S-Ju 


V o) " V u) H(u) 


V M> 


Id (I -I +1 ) 
O X y z 

D(S) 


S**ju» 


F (u) - G (u) H(u) 


V w) - 


I S 2 + 4u 2 (I -I ) 
X o y z 


D(S) 


S-jid 


(55a) 

(55b) 

(56a) 

(56b) 

(57a) 

(57b) 


and Q and Q are the first and the third components of Q in Eq. (28). 

X Z T pA 

For the same reason discussed earlier, the roll and yaw dynamics are 
also unstable (referring to Eq. (48), D(e; has roots in the right-half complex 
plane). The PSD's for this case are obtained again by treating gravity gradient 
and gyroscopic torques as disturbances. Figures 13 and 14 show the roil and yaw 
power spectral densities. 

E.4.3. The Error PSD's for Confiscation B 

The stability problem may be resolved by reorienting the shuttle 
from payload-bay nadir (Configuration A) to nose-down nadir (Configuration B) . 

The advantage of this new configuration is chat it is a gravity gradient stabi- 
lized system. However, there are drawbacks for this configuration. First, it 
will require a larger support-tower for the IS instrument, and the second drawback 
is that the aerodynamic forces and torques have increased significantly over the 
payload-bay nadir case. 
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To obtain the error power soectral density for Configuration B, it 

A 

is only necessary to replace Y and the corresponding notations of the equa- 

r yUzJa 

tions in this section byl-Z I . That is, for instance, Eq. (52b) becomes. 


V“) - 


I,S 2 + 3 0) 2 il Y -ly L ) 


S“ju 


(58) 


Z ~ o 

where G g (‘) here still represents the pitch dynamics. 3ince > I x> Eq. (58) 
is stable. 


The instrument pointing error PSD’s are shown in Figure 15 and the 
race error FSD's are illustrated in Figure 16. 

The computer programs that are used for generating these results are 
included in Appendices C and D. 

F. SUMMARY OF MAJOR FINDINGS 


F.l The Major Environmental Disturbances 

. v 

The major disturbance sources that are modeled are the aerodynamic 

drag torques, the gravity gradient torques, and the gyroscopic torques. The 

unmodeled disturbances are the solar pressure torques, the on-board equipment 

vibrations, crew motions, and venting. Table 2 shows the static disturbance 

2 

torques in N-m and the stochastic torque PSD's in (N-m) -sec for both Configura- 
tions A and B witn a circular orbit of 400 km altitude. 

F.2 The Measurement Uncertainties 

The modeled measurement uncertainties are summarized in Table 3. 

F.3 Ground Track. Errors and Navigation Uncertainties 

The power spectral densities for the ground track errors aud the rate 
errors due to instrument pointing uncertainties for Configurations A and B are 
shown in Figure 17. It is important to note that strong resonances occur within 




—3 

the frequency band of 10 Hz to 10 Hz. The ground track errors at higher 
frequencies, .01 Hz and above, are dominated by the measurement noise. Recall 
that the la measurement noise is 28.28 arc-sec/axls . The corresponding ground 
track error is 54.84 m/axis or 77.56 m/lateral motion. 

The 3a ground track errors due to navigational uncertainties are 
about one order of magnitude greater than the attitude errors with the aid of 
TORSS’; and the error is greater with STDN as indicated in Table 4 [12]. Note 
that the values tabulated in Table 3 are the 3a rms values. 
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Figure 15. Instrument Pointing Error PSD, (rad) -sec 
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Table 2. Modeled Environmental Disturbances 


T5 


u> 


i 


i 


i 



- 

Static, N-m 

2 

Stochastic — PSD, (N-m) -sec 

Sources 

X P 

Y P 

Z P 

X P 

Y P 1 

Z P 

Aerodynamic 

Drag 

Torques 

Config. A 

0 

-.68xl0~ 2 

0 

2xl0 -6 

1.2x10“ 2 

5.4XlO“ 4 

Config. B 

-2.78xl0~ 2 

0 

-3.21xl0~ 2 

6.3xl0 -4 

1.9xlO~ 4 

1.2xl0 -3 

Gravity 
Gradient and 
Gyroscopic 
Torques 

0 

0 

1 

0 

Proportional to attitude errors 


Table 3. Modeled Measurement Uncertainties 


Sources 

Performance (lo) 

Requirement (lo) 

IMU/Shuttle Misalignment 

82 arc-sec 

160 arc-sec 

IMU Resolution 

20 arc-sec 

20 arc-sec 

1MU Noise 

(20) arc-sec 

- 

IMU Derived Rate 

See filter 

dynamic model 

Rate Gyro 

- 

60 arc-sec/sec 

Shuttle/IS Misalignment (lo) 

30 arc-sec 
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Figure 17. IS Ground Track Error and Rate Error PSD 
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Figure 17. (Continued) 
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Table 4. Expected On-Orbit Navigation Accuracies (3o) 


■e- 

ro 


DESCRIPTION 


c 

NAVIGATION 

TRACKING 

UNMODELED 

MINIMUM 

PERIGEE 


POSITION, FEET 



VELOCITY. FPS 


CASE 

SYSTEM 

PERTURBATION 

n.mi. 

RADIAL 

DOWNTRACK 

CROSSTRACK 

RADIAL 

DOWNTRACK 

CROSSTRACK 

1 

STDN* 

NOMINAL 

105 

5,000 

34,500 

5,000 

38.4 

5.9 

9.8 


TDRSS** 


105 

3,600 

18,500 

4,000 

22.1 

3.3 

4.0 

2 

STDN 

SMALL 

105 

3,000 

22,000 

1,500 

' 24.6 

, 3.3 

2.0 




150 

1,200 

5,500 

1,000 

5.8 

1 

, 1.3 

1.8 


TDRSS 


105 

1,800 

10,000 

2,000 

11.0 

■ 1.8 

2.9 




150 

800 

2,000 

1,500 

2.2 

0.9 

2.1 

3 

* 

STDN 

LARGE 

105 

8,000 

80,000 

5,000 

91.2 

8.0 

8.0 

i 

TDRSS 


105 

6,008 

50,000 

5,000 

56.0 

6.2 

7.8 



i 


\ 

\ 


\ 


NOTE: 


The correlation between downtrack position and radial velocity is -0.95 


The correlation between radial position and downtrack velocity is -0.80 
* Spaceflight Tracking and Data Network 
** Tracking Data Relay Satellite System 


! 



-A 


w 




X ' 


I 






IV. PARAMETRIC ANALYSIS OF GEOMETRIC ERRORS 

In section III, the steady state error dynamics of the Shuttle Imaging 
Spectrometer system has been analyzed with the major error sources and disturbance 
effects estimated. In this section, the emphasis is on the geometric error analysis. 
Geometric errors are consequences- of the more direct errors including attitude and 
rate errors, ephemeris uncertainties, misalignment errors, earth rotation and 
curvature, etc. Figure 1 shows how the various errors propagate and how the geo- 
metric errors and the ground pattern distortions can be generated through dynamic 
analysis and simulation. The block that is considered in this section is II. Due 
to mounting and other practical considerations, only the payload-bay principal- 
axis nadir pointing configuration is considered (see Fig. 2). 

The geometric error mapping functions due to the individual error 
sources as well as the aggregated errors are derived. The earth curvature effect is 
incorporated in all of' the results. For the purpose of quick, reference, the key' 
mapping functions are tabulated in Appendix A. - A list of the source code of the 
computer program that has been used for generating the geometric error characteris- 
tic curves is given in Appendix E. 

A. COORDINATE CONFIGURATIONS 

In Fig. 18, the coordinate frame (Xp,Y ? ,Zp) on the Shuttle c.m. con- 
sists of the principal body axes. For the purpose of geometric analysis, another set 
of coordinates is used, i.e. , the (X,Y,Z) frame centered at the nadir point on the 
ground. This frame is the nadir projection of the orbital rotating frame centered 
at the Shuttle c.m. Specifically, the X-axis is in the direction of the projected 
motion on the ground, or the along-track direction, the Z-axis is in the nadir direc- 
tion, and the Y-axis is in the cross-track direction, so that a right-hand coordinate 
systdm is formed. ' 4 ** J 

For the purpose of this report', it Is assumed that the Imaging Spectro- 
meter is attached to the payload bay with its optical axis aligned with the body 
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Zp-axia . For normal IS operations, this Is the desired configuration. However, 
for some IS experiments, such as those for assessing the atmospheric effects on 
the Images, the IS Instrument is required to point up to +45° about the pitch 
axis (Y) from nadir (forward and backward looking) or up to +20° about the roll 
axis (X) (side looking). For those cases, it is assumed that the instrument is 
gimbal mounted. However, since the analysis performed here is not primarily con- 
cerned with dynamics, no specif ic -details are made at this time regarding mounting 
configurations . 

B. GEOMETRIC ERRORS INDUCED BY EPHEMERIS UNCERTAINTIES 

’ Ephemeris uncertainties include radial, along-track, and cross-track 

prediction errors. Tha geometric errors induced by ephemeris uncertainties are 
discussed in the following subsections. 

B.l Geometric Errors Due to Altitude Uncertainties 

When the altitude of the shuttle varies, the ground point will shift in 
both the X and Y directions accordingly. This can be investigated in the following 
two ways . 

B.1.1. With Fixed Viewing Angles 

Referring to Figure 19, when the shuttle flies at nominal altitude h 
(position A), a ground point B with coordinates (Xo, Yo) corresponding to a view 
angle A from the shuttle IS is located. After the shuttle elevates Ah to position 
A', point B' with coordinates (Xo', Yo') corresponding to the same view angle A 
is located. The problem is to determine the shift of image due to the altitude 
change. Consider that when the altitude increases, the image of ground objects 
tends to move toward the nadir point, which causes reduction in Image size and in- 
crease in field of view. Mathematically, the shift of image may be characterized 
by computing AX and AY, where 
AX - Xo - Xo' 

AY - Y - Yo' 


45 


a. 




Starting by assuming (Xo, Yo) known, the task is to determine (Xo*. Yo' 
by first computing the view angle X. Figure 19b is generated by taking a side 
view from Fig. 19a in the direction perpendicular to Plane A’AOC^’. Apparently, 


\fjLo 1 + Yo 2 > 


R sing 


X 



IN FIGURE 4 (b)) 


LEGEND: 

B = A GROUND POINT CORRESPONDING TO 
VIEW ANGLE X WHEN SHUTTLE IS AT 
NOMINAL ALTITUDE {POSITION A IN FIG. 4<b)) 

B 1 = THE GROUND POINT CORRESPONDING TO 
THE SAME VIEW ANGLE AFTER THE SHUTTLE 
ELEVATED Ah (POSITION A* IN FIG. 4(b)) 

(X Y ) = COORDINATES OF POINT B 

Ot o 

(X * Y ') = COORDINATES OF POINT B’ 

O a O 

AX, AY = X AND Y COMPONENTS OF THE POSITION CHANGE 


(a) TOP VIEW (LOOKING TOWARD NADIR) 


. Figure 19; Ground Point Shift in Z and Y Direction due to Altitude Error 







h = NOMINAL ALTITUDE 
Ah - ALTITUDE ERROR 
X = VIEW ANGLE 

A = SHUTTLE'S POSITION AT NOMINAL ALTITUDE 
A' = SHUTTLE'S POSITION AFTER BEING ELEVATED Ah 


Figure 19. (Continued) 
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Substituting Eq. (63) into (64) 




L . h + Ahr .2.1. . 
II -i I sin X sin \ 

-* C66) 

line OB, the angle p is 








When point A' is below point A, the negative value for Ah should be used in 
the above equations. 

However, if one starts with given view angle land p, the values of Xo, To, 

Xo', and Yo* can be readily obtained using Eqs. (68) and (69) directly (set Ah 
to 0 for Xo and Yo) . 

B.l.Z With Fixed Ground Points 

Referring to Fig. 20, the view angle for the fixed ground point B will 
change as the shuttle altitude varies. Follow the same approach of paragraph B.1.1, 
the view angles ft and ft' can be found. 


ft “ tan 


' R + h - 


Jxo 2 + Yo 2 


J 


R 2 - Xo 2 - Yo 2 / 


and 


tan 




J Xo 2 + Yo 2 


(70) 


(71) 


JR 2 - Xo 2 - Yo 2 ^ 


‘ R + h + Ah - Jl 

Using Eqs. (70) and (71), the value of ifi = Q 1 - Q can be obtained. Again if 
point A' is below point A, the negative value for Ah should be used in the above 
equations. 


B.2. Geometric Errors Due to Intrack and Crosstrack Prediction Errors 

The intrack (along X-direction or direction of orbit) and the crosstrack 
(along Y-direction or orbit normal direction) prediction errors will cause the 
ground objects to shift along the X and Y directions, respectively, as shown in 
Fig. 21. To determine the overall geometric errors associated with ephemeris 
uncertainties, the intrack error AX* and the crosstrack error AY* can be incor- 
porated in AX and AY found in Section B.1.1, respectively. 
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W\ 


.r'/l 


R< 


I / 

I / 


/ 


/ 


/ 


EARTH 

CENTER 



LEGEND: 

ft = VIEW ANGLE FROM SHUTTLE TO POINT B 
WHEN SHUTTLE IS AT NOMINAL ALTITUDE h 

a' = VIEW ANGLE ROM SHUTTLE TO THE SAME 
POINT B AFTER SHUTTLE ELEVATES Ah 


Figure 20. 


Change of View Angle of a Ground Point due to Altitude Error 
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( b ) CROSSTRACK ERROR (AY*) 


Figure 21. Geometric Errors Induced by Intrack and Crosstrack 
Prediction Errors 
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C. GEOMETRIC ERRORS INDUCED BY ATTITUDE UNCERTAINTIES 

Attitude uncertainties refer to the angular errors with respect to 
the nominal roll ($) , pitch (6) , and yaw ($) axes. The effects of these errors 
to the images of ground objects are studied in this subsection first on the indivi- 
dual error basis and then on the aggregated basis. For convenience, in the subsequent 
discussions, the term nominal flight condition will be used to signify the condi- 
tion of the shuttle flying in a 400 km circular orbit with the attitude of payload- 
bay nadir pointing. 

C.l Get *etric Errors Induced by Roll Error 

Referring to Fig. 22, consider the nominal flight condition with 

offset angle about the roll axis. The IS slit in this case is directed at the 

Y-axis on the ground. For a given view angle A, the image point (or pixel) on the 

ground is P Q with coordinates (0 ,Y*) . For the same roll offset and view angle the 

ground point P^ is located after a roll error 4 is introduced. The coordinates 

of P' are (0, Y*') . Equivalently, by rotating the optical instrument to the 'left," 
o 

the image recorded on the film appears to move to the "right." Therefore, the 
geometric error may be defined as the displacement of the object image after error 
is introduced relative to the image before the error is Introduced. Mathematically, 
this is F -? or p. , - P. . , . Since, in the case being considered, the 

o o $ o> a <p Q -Hp , A 

image is on the Y-axis, one has 

AX - 0 

AY - Y* - Y*' 
where Y* and Y*' are, 


(72) 

(73) 


Y* - -R sin (A + $ ) | ( 1 + 


(l + t)co» <* + * 0 > -V 1 "( X + |) 


sin (A + $ ) 
(74) 
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LEGEND: 

X = VIEW ANGLE 
+o = ROLL OFFSET ANGLE 
4> = ROLL ERROR 

P = THE GROUND POINT CORRESPONDING TO 
° VIEW ANGLE X AND ROLL OFFSET 

p o = THE GROUND POINT CORRESPONDING TO 
THE SAME X AND <t> AFTER ROLL ERROR 
4> IS INTRODUCED ° 


' ' - * 
Figure 22. Shift of Ground Point Induced by Roll Error 
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f.-.3eg5>e 



M Sii 


Y* # * -R sin (X + <f> + *) 

o 


[('+!) 

-J'- (' * if sl " 2 (1 ♦ K + *>_ 


cos (X + $ + $) 
o 


(75) 


C.2 Geometric Errors Induced by Pitch Error 

Consider the case of pitch attitude error. Let 0 q be the offset angle about 

the pitch axis (desired attitude), and 0 be the pitch error. Referring to Fig. 23, 

• - 

for a given view angle X, projecting the slit onto the XY-plane forms two line 
segments: BC, corresponding to null pitch error (0=0), and B'C', correspond- 

ing to arbitrary pitch error 0. Projecting along the line of view AC onto the 
ground, P(X*,Y*) is obtained; similarly, P , (X* , ,Y*') is obtained, corresponding to 
line of view AC' . The geometric error here is defined the same way as that for 

the roll case, i.e., ? 0 - P 0 +0>x> or 

o o 

AX = X* - X*’ (76) 

AY = Y* - Y*’ (77) 

The formulae for computing X*, Y*, X*’, and Y*’ are derived as follows. In 
Fig. 23, assume that h, 0, 0 q , and \ are known. The lengths of the line segments 


OB and OB' are 

OB = h tan 0 

and 

OB' 

and that of AB and AB' are 


h tan (0 + 3) 

o 


AB = h 


cos 6 


AB’ 


cos (6 + 0) 

o 


(78a) 

(78o) 

(79a) 

(79b) 
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LEGEND-. 

X = VIEW ANGLE 
e o = PITCH OFFSET ANGLE 
9 = PITCH ERROR 

P - THE GROUND POINT CORRESPONDING 
TO VIEW ANGLE XAND PITCH OFFSET 9 0 

p ‘ = THE GROUND POINT CORRESPONDING 
TO THE SAME XANO 9 0 AFTER 
PITCH ERROR 9 IS INTRODUCED 


Figure 23. Shift of Ground Point Induced by Pitch Error 
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I 


BC and B'C' , with respect to Y-axis, can be obtained as 


BC = -AB can X 


-h tan X 
cos 9 


B’C' - -AB* tan X 


-h tan X 

cos(9 +9) 
o 


implies that the signs of BC and B'C’ are opposite to that of X. I and l' are. 


referring to Fig. 23, 


I — ? - — 2 1 2 2 

2 =JOB + BC = h sec 9 q yjsln 9° + tan ,\ 


(81a) 


l' -Job' 2 + B'C' 2 - h sec (9 + 9) J sin 2 (9 + 9) + tan 2 


the angles t and £ (see Fig. 23) can then be determined. 


-1 ( l\ -1 / 

t = tan 1—1 = tan \sec o q 


12 2 
f sin 9 + tan X 
1 o ) 


|-~*j - tan ^ £sec (9 q + 9) -/sin 2 (6 q 9 9)+ can 2 X^ 


Finally, 


X* = Jx * 2 Y * 2 ™ 


R jjl + |) cos r - J 1 - (l + f) S m 2 t] — = 


t sin 3 


Y* = Jx * 2 + Y* 2 


-R ^1 + ^ j COS T - . 


sin 5 + tan X 
o 


fl - (l t|) sin 2 T J - sin T Can ^ (83b) 

■^sin 2 9 q + tan 2 X 


c*\ 
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-1 11 t 



X*' - -Jx*' 2 + Y t2 -“ 

• R [(> + 1) “* « -J 1 - (‘ + lf sl » 2 J 


Y** - ,/ x *’ 2 + Y* ,2 ~T“ 


sin C sin (9 + 6) 

o 

^sin 2 (0^ + 9) + can 2 X 

(83c) 

sin g tan X 

V sin 2 (0 q + 9) + tan 2 X 
(83d) 

However, if one starts with (X*, Y*) given, then X and 9 q would have to be coa- 


-R 


( I+ !) "■ 5 - J‘ - (> * ii ) 2 


sin 2 5 


puted. In that case. 


X =» -tan 


-1 


Y* 


(r + h - J 


2 2 2 
- X* - Y** 




X* 


J 


R 2 - X* 2 - Y* 2 / 


and 


9 = tan 

o 


•1 /____X*_ 

y r + h - Jr 2 - 


(84a) 


(84b) 


9 2 

X*^ - Y* 


Once X and 6 are determined, t and E can be computed using (82), and X*' 
o 

and Y*' are determined using (83) . 


C.3 Geometric Errors Induced by Yaw Error 

For a given yaw offset angle, \ l> o , and the view angle, X, the ground point 
p(X*, Y*) is located. A new point P'(X*', Y*') is found as yaw attitude error, 
<Ji, is introduced (see Fig. 24) . The values for X*, Y*, X*', and Y*' can be com- 


puted as follows: 


X* 


(l + |) cos X -Jl - (l + slo 2 


R sin X sin ’1/ 


(85a) 
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LEGEND; 

* = YAW OFFSET ANGLE 
o 

+ = YAW ERROR 

P = THE GROUND POINT CORRESPONDING TO VIEW ANGLE LAND YAW OFFSET ♦ 

o 

p* = THE GROUND POINT CORRESPONDING TO THE SAME X AND * AFTER 
YAW ERROR + IS INTRODUCED 


Figure 24. Shift of Ground Point Induced by Yaw Error 
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A 




i - 



1 

i 

* 


l 


✓ 


a. 



and 


Y* - - 


X** 


(l + |) cos X -J 1 - (l + |) 2 
(l + |) cos X -Jl - (l + | f 


sin ^ X 


sin ^ X 


R sin X cos (85b) 


R sin X sin (ii + ty) 
o 


(86a) 


Y*' 


- ^|l + ~j COS X -^1 ~ ^ + R j sin ^ A R Sin * COS (t^ o + tj;) 


(86b) 


From Eqs. (85) and (86), we can obtain 
AX =» X* - X*’ 

AY = Y* - Y*’ 


C.4 Geometric Errors Induced by Roll, Pitch, and Yaw Attitude Errors 

C.4.1. The Shift of IS Line-of-Sight due to Attitude Error 

Referring to Fig. 25, the shuttle is flying at an arbitrary attitude 
(A, 9, <l>) - Consider the 3-2-1 sequence, i.e., the yaw, pitch, and roll rotation 
sequence. After yaw and pitch rotations, the Imaging Spectrometer line-of-sight 
will shift to AB, where B(Xj, Y^) is its intersection with the XY-plane. The 
^1’ coordinates are determined as follows: 

(87a) 
(87b) 


X, = X + h ran S cos ^ 
1 o 


Y = Y + h tan 9 sin 
l o 


Through a roll angle rotation, point B will move to CCXj, Y^) on the XY-plane, 


where 


X_ = X + h tan 9 cos ^ + — 

2 o cos 0 


Y_ = Y + h tan 9 sm ^ r 

2 o cos 9, 


tan $ sin v 


tan <{> cos ’> 


(88a) 

(88b) 
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(TANGENT TO 
THE PROJECTION 
OF THE ORBIT) 


/ PROJECTION 
OF THE ORBIT 


Figure 25. 


The Shift of Ground Projection of the IS Line-of-Sight due to 
Attitude Error , 
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The Interest here is to determine the coordinates of O’(X.Y) corresponding 
to the final line-of-sight projected onto the ground. By taking a view normal to 
the ACO-plane, Figure 26 is obtained. The quantities e, i, and P are obtained as 


follows. 


l *» P. sin C 


rl (gp) 

sin c ^(l + cos e ~J l ~ (* + |) si ° 2 c ^ 


p => sin 


-l/ Y 2 


- 1 / __2 


14 + 4 1 


Hence, the coordinates of X and Y are. 


X = l cos P ** R sin c 




Y=lsinP=Rsinc 


(> * 1) “ * -f- M)* s " £ ] gff 2 


where X 2 and Y 2 are given in Eq. (88). 

Soto that this derivation 1. general enough so that the tor-ulae are valid 
arritude errors with errors plus offsets, i.e., replacing 


if one replaces the attitude errors with errors P ius 
*, 0,and * by + *. * Q * and respectively. 




• V\AV< U. 'K’fc.JXX 


■« *1 L.- 







f 


*v 


i 

C.4.2. The Shift of Ground Image Due to Attitude Error 

For simplicity, $, 8 t and <f> are used to represent the attitude angles I 

1 

1 

as sums of attitude offsets and attitude errors. To fix the derivation, the i 

i 

yaw ('ji), pitch (8), and roll ($) Euler sequence is assumed. 

In this subsection, the results of subsection C.4.1 are extended to cover the 
entire field of view of the IS slit. First consider the case of a nominal flight 
with zero attitude offset. Referring to Fig. 27 , let 0 be the nadir point wh'ich 
is the origin of the XY-Frame and the X'Y'-Frame. The X'Y'-Frame is formed by 
rotating the XY-Frame through an angle tp about the yaw axis. Let P be a point 
corresponding to a view angle X (negative value shown in Fig. 27 ) on the Y-axis , 

before any rotation occurs. After a if; rotation, for the same view angle, the IS 
is sighting point P^. Similarly, after a pitch (8) and then roll (if)) rotation, ,~ 

the points ? 2 and then P 3 are sighted, respectively. - 

The coordinates of F_ in the X'Y'-Frame are, from the results of subsections C.l 

2 n 

and C.2 and Fig. 23, 


R sin 8 sin 


x; = 


h 1 ‘ * 1 ) cos h - t / 1 - (‘ + 1 ) “° 2 


I 


2 2 
sin 8 + tan X 


(91a) 


R can > sin C, 


Y 2 


_(l + %) cos c 2 ( l + r) s1r2 C 2, 


/— 2— 

•Jsm 9 + 


(91b) 


2 

tan X 


where the angle ^ is the an gl e £ defined in Fig. 23 corresponding to point ? 2 
here, and 


S 2 - 


tan 


1 |sec 9 Jsi 


sin 9 + tan X 


(91c) 
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J 


The coordinates of in the X*_Y f -Frame are 


R sin 0 sin £- 


x f 

X 3 


Y' 


i 1+ i ) cos £ 3 -J 1 - l 1 + r) 3i ° 2 h. 


~ 2 


(92a) 


/sin 0 + tan (X + <)>) 



/ \ / / .2 

R tan (X + 4) sin $ 3 

_( 1+ l)cos ( 1+ l) sl ° 2 5 3. 


Jl - 2 - 2 


in 0 + tan (X + $) 


where ? 3 is similarly defined as ^ and 

£ 3 ■ tan * ^sec 0 Jsln 2 0 + tan^ (X + 4>)j 
or, finally, in the XY-Frame 


X 3 = X^ cos ip - Y' 3 sin ip 

Y 3 = X^ sin p + Y^ cos ip 

where and Yj are given in (92). 

The geometric error caused by the attitude errors is therefore, 
AX = X - X 3 

AY = Y - Y, 


(92b) 

(92c) 

(93a) 

(93b) 


(94a) 

(94b) 


where 


X = 0 


Y = -R sin X 


[(i + r) cos x - J 1 - ( 1 + r) sin ^\ 


Now, consider the case of the shuttle (or the IS) flying at a nominal attitude 

(* , 0 , * ) with attitude uncertainties of ($, 0, i|<). Let P (X ,Y ) be a ground 
o o o o o o 

point that is sighted by the IS with a view angle X and attitude (d> 0 , 9 q , ’P q ) ■ 
After (<$, 9, \ P ) attitude rotations from the nominal, the ground point P 3 (X 3 ,Y 3 ) 
is sighted with the same view angle. In this case, the geometric error is 


' i 

, I 

n 
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* 


where 


and 


where 


and 


AX = X - X. 
o 3 

AY = Y - Y, 
o 3 


X • X* cos - Y* sin ili 

0 0 0 0 o 


Y ■ X* sin ^ + Y* cos 
o o oo o 


X- * Xi cos (it> + i>) - Y' sin + !^) 
3 3 0 3 o 

Y- ■ X' sin ('i; + i{i) + Y' cos Ui + i>) 
3 3 0 3 o 


X’ 

o 


Y’ 

o 


R sin 0 

sin 5 

(l + cos 5 --/l - I 

i + |V 

' sin 2 C 

o 

o 

_\ R/ o \ ' 


0. 


J sin 2 0 + tan 2 (X + * ) 


(95a) 

(95b) 

(96a) 

(96b) 

(96c) 

(96d) 

(97a) 


R tan (X + (J ) sin C 
o o 


f 1 + 1) cos s o -j l - ( 1 + i) sln2 g o. 


Jain 2 0 + tan 2 (X + ) 

1 n n 


E = tan 
o 


i" 1 (sec 0 sin 2 6 + tan (X + $ )) 
\ o o o / 

V ^ 


(97b) 


(97c) 


X 3 


Y 3 



/ h\ i 

, | 

h\ 2 2 

R sin (9 + 9) sin 

o J 

1 1+ r) C0S 5 3’V 

1 - 

. 1 + r) sin ? 3. 


J sin 2 (9 + 0) J- tan 2 (X + $ + $) 

o o 


R tan (X 


+ ♦<, + «) si * Lt 1 + 1) c ° s h-J 1 - ( i+ 1) 2 


(98a) 


sin 2 5- 


J: ■ 2 * ~-- 2 


sin (0 + 0) + tan (X + $ + $) 

O o 


' 1 
> -’-i 


(98b) 


and where 


m tan * £sec (0 q + 0) ^sin 2 (9^ + 9) + tan 2 (X + 4> o + 


fca rl 
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« 


To include errors caused by altitude uncertainties, one can replace 

h by h Q in Eq. (97) , and h by h Q + Ah in Eq. (98) , where h Q is the nominal or 

estimated altitude and h + Ah is the actual altitude. 

o 

D. GEOMETRIC ERRORS INDUCED BY ATTITUDE RATE ERRORS 


Attitude rates cause attitude angle changes which in turn cause geometric 
errors. To account for attitude rates, one integrates the rates to obtain the 
attitude angles* and then uses the formulae derived to obtain the corresponding 
time-varying geometric errors. The instantaneous attitude angles are: 

r t 


it ' = 4> + 

O 


Jdx 


(99 a) 



(99b) 


(99 c) 


E. GEOMETRIC ERRORS INDUCED BY MISALIGNMENT ERRORS 


The misalignment errors can be incorporated in the attitude errors. Let 

(* , 9. , U. ) be the attitude bias or misalignment, then the effective IS attitude 
b b b 

angles will be 


: = t + j. + t (100a) 

o b 

9' = 9 + 3. + 9 • - ' (100b) 

o b 

J = > + C-. + J (100c) 

O D 


) 

i 


) 


Gc 



The corresponding geometric errors can be determined using the formulae derived 
in subsection C. 

F. GEOMETRIC ERRORS INDUCED BY EARTH ROTATION 


The effect of earth rotation on images varies with the position of the Shuttle 

relative to the earth and the Shuttle orbital elements. In general, ground 

images become skewed as a result of earth rotation. 

Referring to Fig. 28. Lee 0 be the nadir point at which the latitude is 5 . 

Let a> e be the spin rate of the earth, where * ~ 8 ^q ' q * 7.2722 x 10 ^ rad/sec. 

Let be the linear velocity of the earth at 0, 

V - ui R cos c ( 101) 

e e 

Let T be the orbital inclination at the equator, and let A be the orbital inclina- 
tion to the local meridian at 0. From rig. 29 and Appendix B, one has the follow- 
ing spherical geometrical relation, 
cos f = cos C sin A 


or 


sin A = 


cos r 
cos C 


(102a) 


and 


cos A = J l - sin 2 A - (~-0" 


(102b) 


From Fig. 29 again, \ is also the angle between the linear velocity vector and 


the Y-axis. Hence, velocity components V eX and are, using (102), 


and 


V „ = V sin A 
eX e 


V „ = V cos A 
eY e 


v cos r 
e 

cos C 


i-(Sr) 


(103a) 


(103b) 


*1 
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LEGEND: 

= ANGULAR VELOCITY OF EARTH SPIN 
C = LATITUDE A1 POINT O 
v e = LINEAR VELOCITY OF THE EARTH AT POINT O 


Figure 28. Linear Velocity of the Earth at the Nadir Point 
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Figure 29. Earth Sperical Geometry 
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During the time interval from t^ to the ground point has moved from 
P(X, Y) to P'(X* , Y'), the corresponding position changes are* 

'•t^ 


AX 


AY 


. P 2 Ve~a_L dt 

X I cos £ 
x 

f 2 v*- f 2 W 1 - (fftr)‘ 

Jt i J h 


dt 


(104a) 


(104b) 


Note that for small time intervals, 5 may be considered constant and (104) may be 
approximated as 

(105a) 


AX s - £22 -7 V e (t, - t ) 
cos ? e 2 1 


* Y V ■ (^H) 2 


and 


X' ■* X + AX 
Y' - Y + AY 


(105b) 

(106a) 

(106b) 


G. NUMERICAL RESULTS: GEOMETRIC ERRORS AND ERROR SENSITIVITIES 

Numerical results of error sensitivities and geometric errors due to various 
direct errors have been generated. These data are summarized in Tables 5-10 and 
discussed in the following subsections. 

G.l Geometric Errors Due to Altitude Uncertainties 

Table 5 shows the geometric error for a 1 km altitude change and error sensi- 
tivity caused by altitude uncertainty. The nominal altitude of 400 km was used 
tc generate these data. This result is also plotted in Fig. 30 as a function of 
view angles for altitude changes ranging from 1 km to 4 km. As indicated in this 
figure, the geometric errors are quite linear with the view angles and the altitude 
errors for the range shown. 
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For Che Imaging Spectrometer illustrated in Fig. 4, the view angles 

o 

are limited by the field of view of +.825 . Using the various expected on-orbit 
navigation accuracies shown in Table 4 , and the sensitivity data of Table 5, the 
corresponding geometric error can be determined. For instance, from Table 4, the 
3a altitude uncertainty, with TDRSS in the tracking system for the 150 lyni. orbit, 
if the small unmodeled perturbation is used, is 800 feet (or ^244 m) and with the 
STDN system and large unmodeled perturbation and 105 nmi. orbit, the 3a altitude 
uncertainty is 8000 feet (or ^2440 m) . From Table 5 , the corresponding 3a geo- 
metric errors for the largest view angle (.825°) are 3.51 m and 35.1 m, respectively. 
Kote that the nominal pixel size is 30 m. 

G.2 Geometric Errors Due to Roll Uncertainties 

The geometric errors for view angles ranging from -9° to +9° due to a 
roll error of 1° are tabulated in Table 6. The error sensitivity at the nominal 
attitude is also tabulated in Table 6. By comparing the errors for 1° and the 
sensitivity for the same view angle, it is found that the differences are.2Z or 
less. That is, for small angles, the error sensitivity data can be used to com- 
pute errors. Fig. 31 shows the plots of geometric errors for roll errors of up 
to 5°. Unlike the errors corresponding to altitude uncertainties, these curves 
are relatively flat. The geometric error is about 7 km per 1° of roll error. 

It is noted from Fig. 31, that for small roll errors, the geometric 
errors are nearly symmetrical with the view angles. For larger roll angles, the 
geometric errors are greater for view angles that have the same sign of This 
is more visible from Fig. 32 which shows the plots of geometric errors corresponding 
to roll errors of 1° and 2° from a roll offset of 20°, 



G.3 


Geometric Errors Due to Fitch Uncertainties 


Table 7 shows the error sensitivity and the geometric errors corres- 
ponding to the 1° pitch error for the view angles of up to +9°. The geometric 
errors caused by pitch errors are about the same in amplitude as those caused by 
roll errors. Fig. 33 shows that pitch curves are rather flat, unlike the corres- 
ponding roll curves. Figs. 34 and 35 show the 1° and 2° pitch uncertainty induced 
geometric error curves with 22.5° and 45° pitch offsets, respectively. These 
curves, especially those in Fig. 35, are more similar to the roll curves in shape. 
G.4 Geometric Errors Due to Yaw Uncertainties 

Geometric errors due to rotations about the yaw or the local vertical 
axis (assuming zero offsets) are tabulated in Table 8 for 1° and plotted in Fig. 36 
for up to 5°. Linearity applies for both the view angles and the yaw angles for 
the ranges covered here. Unlike roll and pitch errors, yaw errors induce much 
smaller geometric errors. 

G.5 Geometric Errors Due to Attitude Uncertainties 

Table 9 shows the geometric errors due to 1° error in each of the roll, 
pitch, and yaw axes. The error magnitude for a given view angle is approximately 
the vector sum of the roll and pitch induced errors. 

Fig. 37 shows the plots of geometric errors caused by up to 5° roll, 
pitch, and yaw errors. Fig. 38 shows two curves, one corresponding to 1° pitch 
error and the other, 1° yaw error, with the roll offset of 20°. 

Comparing pitch curve with that of Fig. 33, the 20° roll offset has 
no noticeable effects on the geometric errors. On the other hand, the effects on 
the yaw curve are quite pronounced, as a large bias of approximately 2.5 km results 
due to the 20° roll. 
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Fig. 39 shows the l w roll and I s yaw curves for a pitch offset of 
22.5°. Fig. 40 shows the same curves for 45° pitch offset. 

With the large angular pitch offsets, the yaw curves appear quite 
different from those without pitch offset, as noted in Figs. 36, 39, and 40. With 
these offsets, yaw-induced error curves look quite similar to those Induced by 
roll errors in shape. By offsetting the pitch angle to 45°, the yaw— induced 
geometric errors increase to about sevenfold from zero pitch offset, and more than 
twofold over the 22.5° pitch offset. 

The roll-induced geometric errors were lest, drastically affected by 
pitch angular offsets. For instai.ee, for 22.5° pitch offset, the error has in- 
creased about 7% from that with zero offset, and about 46% for 45° pitch offset. 

G. 6 Ground Point Shift Due to Shuttle Motion and Earth Rotation 

The earth rotation effect on the shift of the ground image varies 
with the latitude of the shuttle, and the combined effects of the earth rotation 
and shuttle motion vary with the orbital inclination in addition to the latitude. 
Table 10 shows the tabulation with latitude ranging from 0° to 80°. The largest 
earth rotation effect occurs at the equator, .46 km in one second, and diminishes 
to .08 km in one second at the latitude of 80° as shown in Fig. 41. 

The combined effects are dominated by the motion of the shuttle since 
,at 400 km altitude the projected ground speed is much greater than the earth rota- 
tion speed. In a one-second period, the image will shift more than 7 km for any 
latitude. 
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Tuble 5 


Geometric Errors and Error Sensitivity 
Due to Altitude Uncertainties 



• Nominal Earth Radius - 6356.785 km (Polar) 

• Nominal Altitude 0 400 km 


View Angle, deg. 

0 

l 

2 

3 

4 

5 

6 

7 

8 

9 

Geometric Error Induced 
by l km Altitude Error, 
km 

0 

0.01746 

0.03492 

0.05242 

0.06995 

0.08753 

0.1052 

0. 1229 

0.1407 

0.1586 

Geometric Error 
Sensitivity, km/km 

0 

0.01746 

0.03492 

0.05242 

0.06995 

0.08753 

0.1052 

0.1229 

0.1407 

0. 1586 
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Table 6 . Geometric Errors and Error Sensitivity 
Due to Roll Attitude Error 


• Nominal Earth Radius “ 6356.785 km (Polar) 
« Nominal Altitude ■* 400 km 


View Angle, deg. 

-9 

-8 

-7 

-6 

-5 

-4 

-3 

-2 

-1 

0 

Geometric Error Induced 
by +1° Roll Error, km 

7.153 

7.114 

7.081 

7.052 

7.029 

7.010 

6.996 

6.987 

6.982 

6.982 

Geometric Error 
Sensitivity, km/deg. 

7. 173 

7.133 

7.097 

7.066 

7.040 

7.019 

7.002 

6.991 

6.984 

6.981 


View Angle, deg. 


Geometric Error Induced 
by +1° Roll Error, km 

6.987 

6.996 

7.010 

7.029 

7.052 

7.081 

7.114 

7.153 

7.196 

Geometric Error 
Sensitivity, km/deg. 

6.984 

6.991 

■ 7.002 

ii 

7.019 

7.040 

7.066 

7.097 

7. 133 

7.173 


' J ' . ' V* 

.X* .•. ... 

'.'it; 
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Figure 32 . 


Geometric Error due to Roll Uncertainty (£) About 20° Roll Offset 
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Table 7. Geometric Errors and Erior Sensitivity 
Due to Pitch Attitude Error 


• Nominal Earth Radius ■ 6356.785 km (Polar) 

• Nominal Altitude « 400 Um 


View Angle, deg. 

-9 

-8 

-7 

-6 

-5 

-4 

-i 

-2 

-1 

0 

Geometric Error Induced 
by 1° Pitch Error, km 

6.988 

6.986 

6.985 

6.985 

6.984 

6.983 

6.983 

6.982 

6.982 

6.982 

Geometric Error 
Senate ivlty , km/deg. 

6.987 

6.986 

6.985 

6.984 

6.983 

6.982 

6.982 

6.982 

6.981 

6.981 


View Angle, deg. 

1 

2 

3 

4 

5 

6 

7 

8 

9 

Geometric Error Induced 
by 1° Pitch Error, km 

6.982 

6.982 

6.983 

6.983 

6.984 

6.985 

6.985 

6.986 

6.983 

Geometric Error 
Sensitivity, km/deg. 

6.981 

6.982 

6.982 

6.982 

6.983 

6.984 

6.985 

6.986 

6.987 
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ERROR 



-9° -8° -7° -6° -5° -4° -3° -2° -1° 0° 1° 2* 3° 4° 5° 6° 7° 8° 9° 

VIEW ANGLE 


Figure 33. Geometric Error due to Pitch Attitude Error ( 8 ) 
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Tabic 8. Geometric Errors m.l Error Sensitivity 
Due to Yaw Attii Error 


• Nominal Earth Radius “ 6356.785 km (Polai) 

• Nominal Altitude m 400 km 


View Angle, deg. 

0 

1 

2 

3 

4 

5 

6 

7 

8 

9 

Geometric Error Induced 
by 1 Yaw Error, km 

0 

0. 1219 

0.2438 

0.3659 

0.4883 

0.6109 

0.7340 

0.8576 

0.9818 

1.107 

Geometric Error 
Sensitivity, km/dcg. 

0 

0. 1219 

0.2438 

0.3659 

0.4883 

0.6109 

0.7340 

0.8576 

0.9818 

1.107 
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Table 9 . Geometric Errors Induced by 
Yaw, Pitch, and Roll Errors 


• Nominal Earth Radius = 6356.785 km (Polar) 

• Nominal Altitude <* 400 km 

• Yaw Error - +1° 

• Pitch Error - +1° 

• Roll Error - +1° 


View Angle, deg. 

-9 

-8 

-7 

-6 

-5 

-4 

-3 

-2 

-1 

0 

Error, km 

9.245 

9.296 

9.353 

9.414 

9.480 

9.551 

9.625 

9.704 

9.788 

9.875 

View Angle, deg. 

l 

2 

3 

4 

5 

6 

7 

8 

9 


Error, km 

9.966 

10.06 

10.16 

10.27 

10.37 

10.48 

10.60 

10.72 

10.85 



S2\ 




ERROR 



VIEW ANGIE 


Figure 37. Geometric Error due to Roll (#), Pitch ($), and Yaw ( ^) Errors 
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Figure 38. Geometric Error due to Pitch ( 0 ) and Yaw (40 Uncertainties About 
20° Roll Offset 
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figure 39. Geometric Error due to Roll (<f>) end Yaw (</0 Uncertainties About 
22.5° Pitch Offset 
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Figure 40. Geometric Error due to Roll (<£) and Yaw Uncertainties About 4 

45° Pitch Offset 1 
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Table 10. Ground Point Shift Induced by 
Shuttle Motion and Earth Rotation 


• Nominal Earth Radius - 6356.785 km (Polar) 

• Nominal Orbit Inclination ■ 85° 

• Nominal Shuttle Ground Speed ■ 7.202 km/sec 


Latitude, deg. 

0 

20 

40 

— 

60 

80 

Ground Point Shift Velocity 
due to Earth Rotation, km/sec 

0,4623 

0.4344 

0.3541 

0.2311 

0.08027 

Ground Point Shift Velocity 
Relative to Nadir Point, 
km/sec 

7.177 

7.175 

7.171 

7.166 

7.162 

Angle of Velocity Vector 
with Y-axis, deg. 

-86.32 

-86.54 

-87.19 

-88.18 

-89.44 

Shift in 0. 1 sec, km 

0.7177 

0.7175 

0.7171 

0.7166 

0.7162 
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H. 3UKMAHY cr KAJCS. FINDIKSS 

In this section geometric errors caused by the "direct error sources" 
are analyzed. Equations that map the direct errors onto the geometric errors are 
derived in Subsections B through G. Extensive illustrations of geometry are used 
to assist the derivations. Direct error effects Included here are those of pointing 
errors and angular rates, ephemeris prediction errors including altitude, intrack, 
and crosstrack uncertainties , earth rotation and curvature, and shuttle and instru- 
ment ml8aligment. For quick reference the key error mapping functions are listed 
in Appendix A. A program list is given in Appendix E. 

The numerical results for both the geometric errors and error sensi- 
tivities are tabulated in Tables 5 through 10 and also plotted in Figs. 30 through 
41. To illustrate the effects of view angles and linearity, an extended range of 
+ 9 ° is used. For the push broom Imaging spectrometer studied here, the field of 
view is limited by +.825° (see Fig. 4) . The following are the specific findings 
resulting from further analysis: 

1. The effects of earth curvature are very small for the applications 
here. For instance, the linear displacement of a 10 km arc is 9. 9999... km. Fig. 



'S' = 10 km -?*= 413.752 ... km 

a = 9.99999 ... km c = 413.440 ... km 

b = 9.99994 ... km a = 400 Ian 


I 

* 

* 


Figure 42. Effects of Earth Curvature 
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2. Altitude uncertainties cause only moderate geometric errors. Table 11 
shows that the worst la geometric error is 11.71 m and the corresponding rate is 

. 133 m/sec. The worst case associates with the STDN system and large unmodeled 
perturbation and 200 km orbit. The least la errors are .22 m and .0032 m/sec 
corresponding-to. the TDRSS system and small uimodeled perturbation and 300 km orbit. 
The effects of other navigation uncertainties such as dcwntrack and crosstrack 
position errors are not small, however. Table 12 shows the 3a uncertainties in 
feet. Downtrack and crosstrack position errors are mapped directly to the 
geometric errors. 

3. Geometric errors caused by altitude uncertainties are, within the range 
of interest here, proportional to the view angle and altitude error (see Table 5 
and Fig. 30) . 

4. The effects of roll and pitch attitude errors are relatively large 
compared with those of altitude errors and yaw attitude errors. Table 13 shows 
thac, for nominal flight conditions, the error sensitivity is approximately 
1.94 m/arc sec, or , v-7000 m/ degree. The yaw sensitivity is very small, from 0 
for 0° view angle to .0279 m/arc sec for maximum IS view angle ±.825°. 

The effects of attitude errors increase significantly for large attitude off- 
sets. For instance, for 20° roll offset (side looking), the roll sensitivity increases 
to 2.28 m/arc sec from 1.94 and yaw sensitivity increases to .75 m/arc sec from 
.0279. The pitch sensitivity, in this case, is nearly unaffected. By offsetting 
the pitch angle to 45° (forward looking), the roll sensitivity increases to 
2.83 m/arc sec, the pitch to 4.36 m/arc sec, and the yaw to 1.97 m/arc sec. 

Table 13 shows the error sensitivities and 1° geometric errors for many combina- 
tions of interest. 
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5. Without attitude offsets and within the IS field of view, the geometric 
errors are almost proportional to roll errors (Fig. 31) , pitch errors (Fig. 33), and 

1 

yaw errors (Fig. 36), and are independent of the view angle except for the yaw cases 
(for which the error is proportional to the view angle). For the cases of large roll 
or pitch angular offset, most of the properties change only slightly (see Figs. 

32, 34, 35, 38, 39, and 40). 

6. The performance of the Imaging Spectrometer is limited by the shuttle 
Inertial Measuring Unit accuracy, the shuttle reference misalignment, and the 
misalignment between the shuttle and the Imaging Spectrometer, unless means of 

- error reduction, such as using ground control points and a precision point mount 

between the shuttle and the IS instruments, are employed. 

j On top of these uncertainties, the shuttle deadband is another source of 

i 

error that can cause gross geometric errors. This problem, of course, can be 
resolved by using a precision point mount. 

I ° 

, Table 14 shows geometric errors caused by the IMU/shuttle misalignment, ' ' ’ 

IMU resolution, IMU noise, rate gyro drift, and combined shuttle/IS misalignment, 
for nominal attitude and attitude offsets. Single axis geometric errors corres- 
ponding to the combined shuttle/IS misalignment are 169 m, 198 m, and 379 m, 
respectively, for nominal, 20° roll offset, and 45° pitch offset, for instance. 

Other cases are shown in detail in Table 14 . 

7. Earth rotation causes images to shift toward the direction of rotation. 

The magnitude of this shift depends on the latitude of the object (see Fig. 41). 

For instance, at the equator, the object moves approximately 462 m in 1 second, 
while at 60° latitude, it moves only 231 m in one second. 
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Table 11. Geometric Errors Due to Expected 
On-orbit Navigation Uncertainty 


Navigation 


Minimum 

la Altitude 

Max. Geom. 

la Radial 

Max. Geom. 

Tracking 

Un mode led 

Perigee, 

Error, 

Error in 

Velocity, 

Error in 

System 

Perturbation 

km 

m 

FOV, m 

m/s 

FOV, m/s 

STDN 

Nominal 

200 

508 

7.32 

3.87 

.056 

TDRSS 

200 

366 

5.27 

2.25 

.032 

STDN 

• - 

200 

305 

mm 

2.50 

m 

Small 

300 

122 

— 

.59 

BS9 

TDRSS 

200 

183 

2.63 

1.12 

.016 



300 

81 

1.17 

.22 

.0032 

STDN 

Large 

200 

813 

11.71 

9.27 

.133 

TDRSS 

200 

610 ! 

8.78 

5.69 

.082 


STDN — Spaceflight Tracking and Data Network 
TDRSS — Tracking Data Relay Satellite System 
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Table 12. Expected On-Orbit Navigation Accuracies (Ref. 12) 


Navigation 

Unmodeled 

Minimum 

3cr Position, 

Feet 

3c Velocity, 

FI’S 

Tracking 

System 

Pertur- 

bation 

Perigee, 

n.ml. 

Radial 

Down 

Track 

Cross 

Track 

Radial 

Down 

Track 

Cross 

Track 

STDN 

Nominal 

‘ 105 

5,000 

34,500 

5,000 

38.4 

5.9 

9.8 

TDRSS 

105 

3,600 

18.500 

4,000 

22.1 

3.3 

4.0 

STDN 

•• 

105 

3,000 

22,000 

1,500 

24.6 

3.3 

2.0 


Small 

150 

1,200 

5,500 

1,000 

5.8 

1.3 

1.8 

TDRSS 

105 

1,800 

10,000 

2,000 

11.0 

1.8 

2.9 



150 

800 

2,000 

1,500 

2.2 

0.9 

2.1 

STDN 

Large 

105 

8,000 

80,000 

5,000 

91.2 

8.0 

8.0 

TDRSS 

105 

6,000 

50,000 

5,000 

56. C 

6.2 

7.8 


NOTE: 

The correlation between dovntrack position and radial velocity is -0.95. 
The correlation between radial position and downtrack velocity is -0.80. 


STDN - Spaceflight Tracking and Data Network 
TDBSS - Tracking Data Relay Satellite System 
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Table 13. Geometric Error Sensitivity 
Nominal Orbit: 400 km Circular 


Error 

Source 

Attitude 

Offset 

Error Sensitivity, m/arc sec 

Geometric Error 
for 1° Angular 
Error, m 

For View 
Angle “ 0° 

For View 
Angle ■ ±.825° 

Roll angle 

No 

offset 

1.939 

1.940 

6,987 

Pitch angle 

1.939 

1.939 

6,982 

Yaw angle 

0 

.0279 

101 

Roll angle 

20° 

roll 

offset 

2.25 

2.28 

2.22 

8,200 

Pitch angle 

1.94 

1.94 

7,000 

Yaw angle 

.72 

.75 

.70 

2,683 

Roll angle 

22.5° 

pitch 

offset 

2.11 

2.11 

7,600 

Pitch angle 

2.31 

2.31 

8,300 

Yaw angle 

.83 

.83 

3,000 

Roll angle 

45° 

2.82 

2.33 

10,200 

Pitch angle 

pitch 

offset 

4.33 

4.36 

15,700 

Yaw angle 

1.94 

1.97 

7,100 

Altitude uncertainty 

0 

.0144 m/m 

14.4 m per 1 km 
altitude change 
max in F0V 
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Table 14, Geometric Error Induced by Shuttle and 
Imaging Spectrometer Measurement Uncertainties 

Nominal Orbit: 400 km Circular 


Source - - - - 

Attitude 

Offset 

Per Axis* 
arc • sec 

Per Axis 
Geometric 
Error, 
m 

2 Axes.t 
Geometric 
Error, 
m 

IMU/Shuttle Misalignment 

No Offset 

82 

159 

225 

IMU Resolution* 

20 

39 

55 

IMT Noise 

20 

39 

55 

Rate Gyro 

60 arc sec/sec 

116 m/s 

165 m/s 

Combined Shuttle/IS 
Misalignment 

87 

169 

239 

IMU/Shuttle Misalignment 

20° 

Roll 

Offset 

82 

187 

245 

IMU Resolution 

20 

46 

60 

IMU Noise 

20 

46 

60 

Rate Gyro 

60 arc sec/sec 

137 m/s 

180 m/s 

Combined Shuttle/IS 
Misaxignment 

87 

198 

260 

IMU/Shuttle Misalignment 

45° 

Pitch 

Offset 

82 

358 

426 

IMU Resolution 

20 

87 

104 

IMU Noise 

20 

87 

104 

Rate Gyro 

60 arc sec/sec 

262 m/s 

312 m/’s 

Combined Shuttle/IS 
Misalignment 

87 

379 

452 

•im 

‘Combined roll and pitch axes. 


100 




V. GROUND PATTERNS AND IMAGE DISTORTIONS 


The extent that direct errors and earth geometric properties affect 
the performance of the Shuttle Imaging Spectrometer are analyzed and extensively 
Illustrated with charts and tables in section IV. The advantages of tables and 
charts are that they carry precise and specific Information which is invaluable 
for design and performance prediction. However, it is difficult to relate this 
information directly to^images of ground objects by human brains. To compensate 
for this and to make direct observation of how errors in the system affect ground 
Images, selected ground patterns as seen through this optical system are studied 
here. 

The patterns that are readily generated and, most Importantly, suitable 
for exhibiting imaging distortions are square grid patterns. 

To cover a large ground that is comparable to the field of view of 
the Imaging spectrometer, a 10 km field is selected and the field is evenly divided 
into 10 segments of 1 km each. Since the imaging spectrometer employs a push broom 
principle, each field is then 10 km wide (cross track) and 30 m "long" (along track). 
However, to see a large ground area, a 10 km x 10 km field is employed with the 
along track grids tagged with time. Considering a perfectly spherical earth, an 
imaginary 10 x 10 grid is painted on the ground. As the satellite flies through 
this region, in the direction from the bottom to the top of the paper, a push broom 
camera should see a pattern just like that of Fig. 43(a). This pattern is referred 
to as the nominal ground pattern. 

Fig. 43 shows how the ground pattern changes when simple attitude error 
occurs without considering the effect of Earth rotation. To show the effect of 
errors, the solid pattern (actual) is overlayed with the "+" pattern (the undis- 
torted nominal pattern). Fig. 43(b) shows that when a 0.1° roll attitude error 
occurs, the IS instrument points 0.1° left of the object field and the image on 
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che "film" has shifted Co Che right; therefore, the solid pattern has shifted to 
the right of the "+" pattern. The same explanation applies also to the 0.1° pitch 
error as shown in Fig. 43(d), except that, for pitch, the pattern image has 
shifted to the "-along track" direction. A 5° yaw error (clockwise rotation) 
will make the pattern appear skewed as shown in Fig. 43(c). 

Fig. 44 shews the effects of simple attitude errors with 20° roll 
offset. With the large roll angle, the nominal pattern Itself has suffered 
cross track distortions, l.e., the pattern image appears shrunk cross-trackwise 
as shown in Fig. 44(a). In Fig. 44(b), (c) , and (d), the "+" patterns represent 
the distorted nominal patterns due to roll offset alone and the solid patterns 
represent those due to attitude errors on top of roll offset. The distortions 
are similar to those found in Fig. 43 except that with the 20° roll offset the 
images appear narrower. 

The effects of simple attitude errors with 45° pitch offset are 
shown in Fig. 45. The effect of pitch offset is similar to roll offset, l.e., 
the nominal pattern Image has shrunk cross-trackwise. One might expect also 
along-track shrinkage. The reason that there is no along-track shrinkage 
is because of the push-broom effect. Shrinkage will occur for a frame camera, 
however. It is noted that the pattern images associated with the 45° pitch off- 
set are narrower than those of the 20° roll offset. This is due to the fact that 
20° is a lot less than 45°. 

The effects of constant attitude rate errors are shown in Fig. 46. 

Rate errors cause accumulation of attitude errors, i.e. , the attitude errors grow 
with time t. Take the 0.1 deg/sec rate error for the roll axis for instance, the 
roll error is 0.1 t; whereas, the 0.1° roll error (Fig. 43(b)) at t ■ 0 is 0; 
therefore, the "+" pattern image coincides with the solid one at a very ] 



103 


a 







105 



(c) WITH 5° YAW ERROR 



(b) WITH 0.1° ROLL ERROR 




Fig. 45. Ground Pattern Distortions With 45° Pitch Oftset and No Earth 
Rotation Effect 
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small t. Shown in Fig. 46(b) Is the solid pattern shift to the right at a constant 
slope with time, which is quite different from Fig. 43(b). Similar comment 
applies to Fig. 46(c) and (d). 

The effects of sinusoidal attitude rates on the pattern Images are 
shown in Fig. 47. For instance, the sinusoidal rate .causing a 0.2 sin (4.51378t) 
degree roll error will make the image appear as the shape "S" , as shown in 
Fig. 47(b). The same function applied to the pitch axis will cause along-track 
distortions with densely packed grid lines followed by loosely packed ones, etc., 
as Illustrated in Fig. 47(d). 

Altitude change causes images to vary in size. Again, it affects the 
cross-track much more than the along-track due to the push-broom, principle. 

As shown in Fig. 48, the solid Image has shrunk cross-trackwise but expanded 
along- trackwise. The former is due to the altitude increase of 40 km which has no 
amplification effect on push-broom camera; the latter is due to the fact that, at 
higher altitude, the orbital rate is decreased, i.e., it will take a longer period 
of time to fly through the same ground area. 

Finally, the effects of earth rotation are shown in Fig. 49. Earth 
rotation will cause the image to skew toward the direction of earth rotation. As dis- 
cussed in section IV, the effect of the image shift is most pronounced when the 
satellite flies through the equator, and reduces as the latitude increases. Earth 
rotation will not cause along-track distortions. 

Computer programs that are used to generate these patterns and distor- 
tions are listed in Appendix F. 
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(a) NOMINAL GROUND PATTERN 
WITH 20 ROLL OFFSET 


(b) NOMINAL GROUND PATTERN 
WITH NO OFFSETS 


(c) NOMINAL GROUND PATTERN 
WITH 45° PITCH OFFSET 
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Ftg. 48. Ground Pattern Distortions With Effect of Altitude Change (increase) of ' 
40 km and No Earth Rotation Effect 
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(a) NOMINAL GROUND PATTERN 
UITH 20 ROLL OFFSET 
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(b) NOMINAL GROUND PATTERN 
WITH NO OFFSETS 





(c) NOMINAL, GROUND PATTERN 
WITH 45 PITCH OFFSET 




Fig. 49. Ground Pattern Distortions gith Effect of Earth Rotation for 40 Latitude 
and Orbit Inclination of 85 










VI. CONCLUSIONS 


1. The effects of earth curvature are very small for the application here (see 
Fig. 42). 

2. Altitude uncertainties cause only moderate geometric errors. The worst lo 
geometric errors are 11.71 m in position and C.133 m/sec in rate with STDN 
and large unmodeled perturbations at 200 'on orbit. The performance improves 
with TDRSS. For the 300 km orbit and with small unmodeled' perturbations, 
the la geometric errors will reduce to 0.22 m and 0.0032 m/sec (see Table 
11 ). 

3. The effects of other navigation errors are significantly greater. The la 
dovntrack errors range from 203 m (300 km orbit) to 8128 m (200 km orbit) ; 
and those for the cross track are 152 m to 508 m (see Table 12). 

4. The effects of roll and pitch attitude errors are relatively large compared 
with, for instance, those caused by yaw errors and altitude uncertainties. 

The error sensitivity is 1.94 m/arc sec or approximately 7000 m/degree. The 
yaw sensitivity is 0 for a 0° view angle and 0.028 m/arc sec for the maximum 
view angle of + 0.825° (see Table 12). 

5. The error sensitivities of attitude errors increase significantly for large 
attitude offsets. For 20° side looking, the sensitivity Increases to 2.28 
m/arc sec for roll errors and to 0.75 m/arc sec for yaw errors, and the pitch 
error sensitivity is almost unaffected. For the 45° forward looking case, 

Che sensitivities for the pitch, roll, and yaw errors increase to 2.83, 4.36, 
and 1.97 m/ arc sec, respectively (see Table 13). 

6. The performance of the Shuttle Imaging Spectrometer is limited by the Shuttle 
IMU (Inertial Measuring Unit) accuracy, instrument misalignment. Shuttle RCS 
(Reaction Control Subsystem) deadband, etc. unless some means of error reduc- 
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cion are employed. Far instance, ground control points may be used to reduce 
navigation prediction errors; and precision point mounts, such as AGS (ASPS + 
Gimbal System) and IPS (Instrument Pointing System), may be used to reduce 
the attitude errors. The single axis geometric errors due -to -the combined 
Shuttle/IS misalignment, for Instance, for normal nadir pointing, 20° side 
looking, and 45° forward looking are 169 m, 198 m, and 379 m, respectively 

e - 

(refer to Table 14). 

7. Earth rotation causes shifts of Images toward the direction of rotation. 

The magnitude of these shifts depends on tha latitude of the object. For 
Instance, at the equator the object moves approximately 462 In one second 
(for 400 km orbit), and at 60° latitude It moves only 231 m In one second. 


8. Distorted Images for selected ground patterns provide revealing information 
of the effects of system errors to the Images of ground objects. 



Annula r Suspension Pointing System 
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APPENDIX A 

GEOMETRIC ERROR MAPPING FUNCTION TABLE 


Pure A: Given the attitude offsets and a view angle, find the coordinates 

of the corresponding ground points before and after the errors 
are introduced. 


Error 

Sources 


Geometric Error Mapping Function 


Altitude At nominal altitude I), a view angle X is given which sights a ground point 
^ ror ’ P(0, Y q ) . After elevating Ah, the same view angle will sight another ground 
point P'(0, TIu. following equations ure for obtaining Y q and Y^: 


Y = -R a in X 
o 


Y' = -R sin X 

o 


(, + CO, 1 - Jl - (t + £) sin 2 

(, ♦ !L±-S!l) cos X -Jl - (l *4^)' sio 2 X_ 


Comments 


P - Earth 
Radius 


Yaw GivenSa yaw offset i (» and a view angle X which sights a ground point P(X*, Y*) , 

Error, Tlie same view angle will sight another ground point P'(X*', Y*') after yaw 

error 1 1> is introduced. The following equations are for obtaining X*, Y*, X*', 
and Y*' : 


j^(l + r) cos ^ - yji ~ + r) 2 ^ R s 

- [(• + 1) “» » -J' - ("'if s1 " 2 R 


R sin X sin ip 


sin X R sin X cos ij> 


R » Earth 
Radius 
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Error 

Sources 


Roll 

Error, 


Pitch 

Error, 

0 


Geometric Error Mapping Function 


X*» = 


Y*‘ 


[(■ • s) 

-M 


cos X 




.2 

>-( 

h \ 2 

1 + sin X 

F 

/ h\ 2 2 

(1 + r) sln x 


R sin X sin (ij) + <(.•) 
o 


]■ 


Given! a roll offset and a view angle X which sights a ground point P( 0 , Y*) . 
The same view angle will sight another ground point P'( 0 , Y*') after roll error 
$ is introduced. The following equations are for obtaining Y* and Y*' : 


Y* - -R sin (X 


Y*’ - -R sin (X + <j> + 4.) 

o 


-Fhl? 


+ *0* [(' + «) “* (1 + »o> " •/' ■ (' + I) s1 ” 2 (> + *0’] 

(‘ ♦ 5 ) 


- j cos (X + 4 > o f <}>) 


sin (X + <J> o + 4 >) 


Given! a pitch offset and a view angle X which sights a ground point 
P(X*, Y*). The same view angle will sight another ground point P'(X*', Y*') 
after pitch error 6 is introduced. The following equations are for obtain- 
ing X*, Y*, X*', and Y*’: 


X* - R 


Y* 


•-( 

'F 

sin 2 tJ ■ 
> 


( 1+ s) 

2 2 1 
sin r 


sin 6 


[sin 2 6 + tan 2 X 

o 


] s in t 

m~ 


tan X 


J 


2 2 
ain 0 + tan X 

o 


Comments 


R - Earth 
Radius 

h ■ Nominal 
Altitude 


R - Earth 
Radius 

h ■ Nominal 
Altitude 

For t and 5, 
refer to 
Fig. 23 . 


n 

if. 

V# 


<r 


• / 
. J 


• | 


.‘■i 


f 


A 


\ \ 



Error 

Sources 


Geometric Error Mapping Function 


Comments 


**' - * [(> + 1) <»• « - J‘ ' (• + if 

" -« [(> + !) “* « -J' - (> + if »*» J t] 

-1 ( I 2 2 \ 

t - tan lsec 8 ^sin 0 + tan X/ 

? “ tan * j^sec (0 q + 6) Jsln^ (0 q + 9) + tan^ xj 


2 I sin £ sin (6 f 0) 

sin 2 ? 1 ° 


J sin 2 (6 + 0 ) + tan 2 X 

* O 


sin C tan X 


Jain 2 (0 + 0 ) + tan^ X 


Yaw 
Error 

*. 

Followed 
by pitch 
error 0, 
then 
roll 
error d> 


Given* a yaw offset ij 1 , followed by a pitch offset 0 , then a roll offset $ . 

o o o 

A view angle X is also given which sights a ground point P o (X o ,Y q ). After 
the errors ij/, 0, and <J> are introduced, the same view angle will sight 
another ground point P^fX^, Y^). The following equations are for obtaining 
X o . Y o , X 3% and Y 3 * 


X' 

o 


R sin 0 sin £ 
o o 


/, h\ „ I 

r~ t 

f h \ : 

> 

* . 2 

j 1 + r) cos 5 o ■ V 

1 - 1 

+ r) 

sin ? 

o_ 


J sin 2 0 + tan 2 (X + 4 ) 

’op o 


R tan (X + <J> ) sin ? 
o 




sin 2 ? 




i 


sin 2 0 + tan 2 (X + 4 ) 

o T o 


R “ Earth 
Radius 

h ■ Nominal 
Altitude| 

For ? 0 and 
C3, refer tol 
Fig. 23 . 
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I 

I 


I 












t'art U: Given the coordinates of a ground point, find the attitude 

offsets, the corresponding view angle and the coordinates 
of the point with the same view angle after the errors are 
introduced. 


Error 

Sources 


Geometric Error Mapping Function 


Altitude At nominal altitude h, a ground point P(X q , Y q ) corresponding to view angle 
j^ r ° r ’ X is given. After elevating Ah, the same view angle X will aim at another 

ground point P'(X' , Y'). The following equations are for obtaining X' and Y': 

„\ 


Comments 


R - Earth 
Radius 


RX sin X[ A h + Ah \ f /, . h + Ah\ 2 . 2 J 

K - 7=== l 1 + -r—J cos x “V 1 " l 1 + —^l sln X J 

JX + Y 
\ o o 

RY sin X Jr , + Ah \ f / . . h + Ah\ 2 , 2 . 

{ o r J ===[ l+ — — ) cos x V 1 “l 1 +_ R“) sln X J 

/X + Y L 







Error 

Sources 


Geometric Error Mapping Function 


Comments 


Roll 

Error, 

♦ 


A ground point P(0, Y*) corresponding to roll offset and view angle A is 
given. After roll error is introduced, the same view angle A will aim at 
another ground point P'(0, Y*'). Tire following equations are for obtaining 
Y*' : 


A + ij> -» -sgn(Y*) tan 
o 


-1 




(r + h -Jr 2 - Y* 2 


) 


Y*' - -R sin - (A + <t> o + <j>) 


[R) 

“7 1 " ( ! + r) 811,2 (X + *o + 


— ) cos (A + + ♦) 


R = Earth 
Radius 

h * Nominal j 
Altitude 






Error 

Sources 


Geometric Error Mapping Function 


Comments 


Pitch 

Error, 

0 


A ground point P(X*, Y*) corresponding to pitch offset 0 q and view anjle X is 
given. After pitch error 0 is introduced, the same view angle X will aim at 
another ground point P'(X*', Y*'). The following equations are for obtaining 
X*' and Y*' : 



R ■ Earth 
Radius 

h ■ Nominal 
Altitude 

For C, refer 
to Fig. 23. 









VJ' 


I , / 


t 



Error 

Sources 

Geometric Error Mapping Function 

Comments 

Yaw 

Error, 

A ground point P(X*,Y*) corresponding to yaw offset and a certain view 

angle is given. After yaw error tj> is introduced, the same view angle will 
aim at another ground point P'(X*',Y*'). The following equations are for 
obtaining X*' and Y * 1 : 

*o “ - tan_1 (-S) 

^* 2+V * 2 \ 

\R + h - </r 2 - x* 2 - Y * 2 / 

X*' ■ + -j| j cos A -^jl - ^1 + ain^ > R sin A sin (i/< o + ip) 

Y*' » - cos A - 1 - |l + "j sin 2 A jj R sin X cos (<l* o + </>) 

R « Earth 
Radius 

h - Nominal 
Altitude 

Earth 
Rotat ion 

Let T be the shuttle orbital inclination at the equator, and ? be the latitude 
of the nadir point. A given point P(X,Y), after time period At, will 
move to P'(X' ,Y’>. The following equations are for obtaining X' and Y' : 

V » uR cos t 
e 

“ “ 86,400 f adl ans/sec 
x ’ + 

V '- Y -/-(coac r ) V '=“ ' 

V e “ Linear 
Velocity of 
Ground Point 
at Latitude 
C 

u ■ Angular 
Velocity of 
Earth Rota- 
tion 

At is 

assumed 

small. 

i 


W"' 





■MW 


Error 

Sources 


Geometric Error Mapping Function 


Comments 


to 

u> 


|Yaw 

Error ip 
followed 
[by Pitch 
Erroi 0 
then 
Roll 
Error <(> 


A ground point P(0,Y*) corresponding to view angle X (with no attitude offset) 
is given. A yaw error ip is introduced first, then a pitch ei ror 0, finally a 
roll error <j>. The same view angle X will then aim at another ground point 
P* <X** ,Y**> . The following equations are for obtaining X*' and Y*' : 


X *» - sgn(Y*) tan 


X 


*( - 'T- ) 

\R + h - VR - Y* / 

i * £sec 9 J sin^ 9 + tan^ (X + <p)j 

[(l -4- |) cos C -J 1 - (l + if siTT 


5 “ tan 

R sin 9 sin C 


J sin^ 9 + tan^ (X + <t>) 


-R tan (X + 0) sin £ 


(l +|) cos C -Ji - (lt|) sin 2 J 


2 • 2 


fsln 0 + tan (X + <p) 


X*' *• X cos ip - Y sin ip 
Y*' “ Xsin ip + Y cos ip 


R «• Earth 
Radius 

h » Nominal 
Altitude] 

For £, refer 
to Fig. 23. 
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APPENDIX B 


^jsas^s&^t 


i 


NAPIER'S ROLES FOR RIGHT SPHERICAL TRIANGLES* 


XL* 


• A right spherical triangle has five variable parts. If these components and 

t 

their complements (complement of T = 90° - I*) are arranged in a circle, as illus- 
trated below 



The sine of any component equa' - the product of either: 

1. The tangents of the adjacent components, or 

2. The cosines of the opposite components 
For our case, from 2: 

sin (90° - 7 ) '■os C cos (90* - A) 

i.e. , 

cos r * >.03 ? sin A 


i *See Ref. 13 
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APPENDIX C 


COMPUTER PROGRAM FOR SHUTTLE IMAGING SPECTROMETER POINTING ERROR 
POWER SPECTRAL DENSITIES FOR CONFIGURATION A — PAYLOAD-BAY NADIR POINTING 


CUT WE0-10/'20>'82-10:23:Z8-<32, J 
12337»IS(U.SRMA^PR0GI321 

1 : PROGRAM IMAGING SPECTROMETER SHUTTLE RIGID MOUNT - R 

2:INITIRL 

3 : VARIABLE T-0.0 

4 : integer fnuh.npoint 

3: 

6 : 

7: comment initialize some parameters 

8: CONSTRNT NPOINT-641 

3: CONSTANT FNUM-39 

10: CALL FOPEN(FNUM) 

11: TFINAL-1 .0»NP0INT-1 .0 

12: CONSTANT FACTOR-O.O 

13: CONSTANT SCAFAC-1 .OE+5 

14: SCFAC2-SCAFAC»SCAF' I C 

15: 3CF2D8-10.0*ALOG10<SCFRCZ> 

16: 

17:C0MM£NT SET VARIOUS MATHMATICAL CONSTANTS 

10: PI-3.14159265 

19: RSRCC-PI'1180.0»3600.0> 

20: RSHZCC-1 ,0si2.0*Pl > 

21 : 

22: COMMENT DEFINE MOMENTS OF INERTIA 

23: CONSTANT IX-1.38E+6 

24: CONSTANT IY-1.00E+7 

23: CONSTANT IZ-1.03E*7 

26: IX2-IX0IX 

27: IY2-IY«IY 

28: IZ2-IZ*IZ 

29: IXYZ-IX-IY4IZ 

30: IXYZ2-IXYZ*IXYZ 

3i: 

32: COMMENT DEFINE WO 

33: CONSTANT WO-O. 001 1313 

34: uo2-wo«wo 

35: W04-W02«W02 

36: 

37: COMMENT DEFINE INITIAL ANGLES AND RATES 

38: CONSTANT THO-O.O 

39: CONSTANT PHI 0-0.0 

40: CONSTANT PSIO-O.O 

41 : CONSTANT thod-o.o 

42: CONSTANT PHI0D-0.0 

43: CONSTANT PSIOD-O.O 

44: THO2-TH0«THO 

45: TH002 - THCD»THOD 

46: PHI02-PHI0»PHI0 

47: PHI 002 -PHI CD*PHI OD 

48: PSI02-PSI0«PSIO 

49: PSI0D2 -PSI OD*PSI OD 

30: 

31 : COMMENT DEFINE O'S 
32: CONSTANT 0X-2.0I4E-6 

33" CONSTANT CY- 1.25 IE -3 

34: CONSTANT 0Z-5.438E-4 

33: 

36-.C0MMENT DEFINE RATE FILTER PARAMETERS 

37: CONSTRNT HD-0.2513 

38: CONSTANT ZETRD-0.8 

39: UD2-MD*W0 

60: U04-WD2»WDZ 

61 : ZD2-2ETAD»2ETRD 

62: 

63: COMMENT DEFINE VARIOUS CONSTANTS 
64: CONSTANT R-3.35E-'0 
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S5: CONSTANT S1G30-87.32 

EE: SIGBO-SIGBO»ASRCC 

E7: SIGB02-SIGBO«SIGBO 

69: CONSTANT TAU-O.OZ3 

S3: TRU2-TAU«TAU 

70: DELH-1.0 

75: 

72: COMMENT DEFINE U CONSTANTS AND INITIALIZE U 

73: CONSTANT WL0-S.2e318531E-7 

74: CONSTANT WFACT-1 .029200527 

75: H-NLO/WFACT 

7E: 

77:EN0 
7B:0YNAMIC 
79:CINTERVAL CI*1 .0 
80: 

81: COMMENT COMPUTE U 

fl9* U«U«UTO^T 

03: PROCEOURBL(DELM-M) 

84: DELW‘0.0 

85: IFCM.LE.O.) DELW-1.0 

85: END 

87: WHZ»M»RSMZCC 

88: FREO«AL0G1O(MHZ) 

89: HZ-N«W 

30: H4-H2»U2 

91 : 

92:C0MMEtT COMPUTE HW2 AND DW2 
93: HU2-t.0xCTRU2»W2+1.0> 

94: HDW2-WD4»W2/( (U2-WD2 )»*2+4 . 0»ZD?*W02 *W2 5 

95: DW2-IX»IZ*W4-W02«<IX«IZ*IY2+3.P IY»IZ-3.0*IZ1>*H2 

98 : COMMENT 0W2-DW244.0»W04»CIY-IX)»(IY-IZ> 

97: DWZ-DWZ/'SCAFAC 

98: DWZ-0U2»0W2 

99: 

100. COMMENT COMPUTE F'S 

101 : COMMENT FTHW2«HW2»SCFRC2X((1Y»W2+3.0»U02»(IZ-IXJ)«»2> 
102:C0hMENT FRU2->(IZ»U2-U02«<IY-IX)>«»2»HW2/DU2 
103: CO WENT F YW2 ■ C 1 X4W2+4 . 0»W02« CIZ-I Y) >*»2«HU2/D142 
104: FThW2«HW2*SCFAC2^t <; y»M2J»*2> 

105: FRW2* CIZ»W2 )*»2»HU2-'DU2 

108: FRYU2*UC2»IXYZ24W2*HW2xDW2 

107: F YW2 • t I X«W2 ) «*2 *HW2xDW2 

108: 

109: COMMENT COMPUTE PSD'S 
110 : 

1 11 :COMMENT PITCH PSD COMPUTATIONS 
112: PPU1-FTHW2»IY2*TH0D2 

113: PPW2-FTHW2*! Y2*TH02»W2 

114: PPW3-FTHH2»OY 

115: PP14-R+SI GB02«DELU+CPPU1 +PPW2+PPW3 1 /’SCFAC2 

11E: PPMD-PPW«HDW2 

117: 

1 18: COMMENT ROLL PSD COMPUTATIONS 
119: PRN1*FRw2»IX2»PhI0D2 

120: PRH7»FRWZ«IX2»PHI02«W2 

121* PRW3-FRU2»U02»IXYZ2»PSI02 

122: PRI44*FRU2»0X 

123; PRW5-FRYW2»IZ2»PSI0D2 

124: PRUG.FRYU2»IZ2»PSI02»U2 

125: PPM7-FRYW2»U02»IXYZ2»PHI02 

128: PRa8-FRYW2«02 

127: RS-H*SIGB02«DELU 

128: PRU-RS* (PRU11 ♦PRW24pfiN34PRW4*PRW54PRWE4PR147+PRW8 J/-SCFACZ 

123; PRWD-PRW*MDW2 

130: 
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131 : COMMENT YAW PSO COMPUTATIONS 
132: PYW1-FYW2*IZ2»PSI0D2 

133: PYW2-FYW2»IZ2»PSI02»K2 

134: PYW3 -FYW2».J02*I XYZ2»PHI 02 

133: PYW4-FYK2»0Z 

136. PYW5»FRYW2«IXZ*PHI0D2 

137: PTW6»FRYW2«IX2*P4I02«J2 

138: PYW7*FRYW2»WO2«IXYZ2»PSI02 

139: PYW8-FRYW2»OX 

140: PYW»RS+ 1 PYW1 +PYM2+PYH3+PYW4+PYW5+PYW6+PYW7+PYW8 ) XSCFAC2 

141: PYW0-PYW«HDW2 

142: 

143: 

144:C0MM£NT PREPARE VAPIA3LES FOR OUTPUT 

143: RSIG-P+SIG802»DELW 

146: RSIGDB-10.0»RLnG10(RSIG)+FBCTOR 

147: PPW30B* 1 O. 0« AL0G1 0( PPM3 ) -SCF2DB+FACT0R 

148: PPWDB-10.0»AL0G10(PPW;+FACT0R 

1 43 : PRW4DB* 1 0.0» AL0G1O(PRW4 > -SCF2DB+FACT0R 

1 SO: PRW8DB* 10. 0*RL0G1O(PRW8 > -SCF2DB+FPCT0R 

151 : PRWDB*10.0«ALOG10(PRW)+FACTOR 

152: PYW4DB-10.0*RLOG10<PYW4>-SCF2DB+FACTOR 

133: PYW8DB* 1 O. 0»AL0G1 0( PYW8 ) -SCF2DB+F ACTOR 

154: PYWDB-10.0»ALOG10(PYW)+FACTOR 

133: 

136: PPWDD8-10.C»ALOG10<PPWC»+FBCTOR 

137: PRWDD3- 1 0. 0* AL0G1 0( PPWD ) +F ACTOR 

138: PYWDDB • 1 0 . 0» ALO G 1 0 1 PYMC ) ♦FACTOR 

1 39: 

16O:C0HM£NT SRVE NUMBERS IN FILE 
• 61: CALL FSfiVE«FREO,PPWDB,PRWDB,PYWOB,FNUm 

162: 

163: 

164: TERMT(T.GE.TFINAL) 

165: DERIVATIVE 

166: ALGORI THM IALG*3 

167: OO-INTEG(I.O.O.O) 

168: END 
169:END 
170:TERMINAL 
171: 

172:C0MMENT CLOSE FILE 
173: CALL FCLOSE(FNUM) 

174: END 
173: END 


EOF: 175 
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0:>ELT HED-10'20^82-10:30:02-C3,) 

1 2337»ZS( 1 1 . SRnft'FSUBSO I 
IS #fOP,IS FF.F0 D EN 
2: SUBROUTINE FOPEN(N) 

3: REMIND N 

« S RETURN 

5: END 

6:#F0R,IS FF.FSRVE 

7; SUBROUTINE FSfivE <M,P1 ,P2,P3,N) 

8: WRITE(N,100)M,P1,P2,P3 

9: 100 FCRHflTMGH.8) 

10: RETURN 

11: END 

12:*F0R,IS FF.FC'.OSE 
13: SUBROUTINE FCLOSEINJ 

14S ENDFILE N 

IS: RETURN 

IS: END 

EOF: 16 
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APPENDIX D 


COMPUTER PROGRAM TOR SHUTTLE IMAGING SPECTROMETER POINTING ERROR 
POWER SPECTRAL DENSITIES FOR CONFIGURATION B — NOSE-DOWN NADIR POINTING 


ELT WED-1<W<V02-1O:36:«t-C1«, J 
12337*IS(1 J.ERMB-/PR0GC14) 

1 : PROGRAM IMAGING SPECTROMETER SHUTTLE RIGID MOUNT - B 

2:INITIRL 

3: VARIABLE T-0.0 

4: INTEGER FNUM,NPOINT 


COMMENT INITIALIZE SOrE PARAMETERS 
CONSTW'T NPOI NT-641 
CONSTANT FNUM-39 
lALL FOT-SNtFNUM) 

TFINRL-1 .0»NP0INT-1 .0 
CONSTANT FACTOR-O.O 
CONSTANT SCAFAC-1.0E+10 
SCFHC2-SCAFHC«SCAFfiC 
SCF20B-10.0»ALOG10<SCFACZ> 

COMMENT SET VARIOUS MATHMATICRL CONSTANTS 
PI -3. 14159265 
RSRCC-PI/l 180. 0*3600.0) 

RSHZCC • 1 . Cs C 2 . 0*PI ) 

COMMENT BEFIK MOMENTS OF INERTIA 
CONSTANT IX-1.00E-7 
CONSTANT IY-1.05E+7 
CONSTANT IZ-1.38E+6 
IX2-IX»IX 
IY2-IY*IY 
IZZ-IZ*IZ 
ir-Z-IX-IY+IZ 
IXYZ2-IXYZ«IXYZ 

COMMENT DEFINE NO 

CONSTANT N0-0.001 1315 
U02 -U0*U0 
U04-U02*U02 

COMMENT DEFINE. INITIAL ANGLES AND RATES 
CONSTANT THO-O.O 
CONSTANT PHIO-O.O 
CONSTANT PSI 0-0.0 
CONSTANT THOD-O.O 
CONSTANT PHIOD-O.O 
CONSTANT PSIOD-O.O 
TH02-THO*THO 
TH0D2-TH0D*TH0D 
PHI 02-PHI 0*PHI0 
PHI0D2-PHI0D*PHI0D 
PSI02-PSI0-PSI0 
PSI 0D2 -PSI ODftPSIOD 

COMMENT DEFINE 0 # S 

CONSTANT 0X-1.970E-4 
CONSTANT GY-1.2B6E-3 
CONSTANT OZ-6.3Z0E-4 

COMMENT DEFINE RATE FILTER PARAMETERS 
CONSTANT WD-0.2513 
CONSTANT ZETAO-O.B 
U02-UD*WD 
UD4-UD2*U02 
ZD2-ZETAD*ZETAD 

COMMENT DEFINE VARIOUS CONSTANTS 
CONSTANT R-9.55E-I0 


63: CONSTANT SIGBO'87.32 

68: 5IGB0-SIGB0»ASRCC 

6?: SI GB02 *SI CBO«SZ GEO 

66: CONSTANT TAU»0.023 

69: TAU2-TAU»TfiO 

70: DCLM-1.0 

71: 

72:C01*'ENT DEFINE u CONSTANTS AND INITIALIZE N 

73: CONSTANT NL0-6.2S318531E-7 

74: CONSTANT Wf ACT* 1 .029200327 

75: H-MLO'KFACT 

76: 

77. END 

7d:DYNRMIC 

79 : Cl NTE AVAL CI*1.0 

80 : 

81:C0r**NT COMPUTE W 

82: H-H»WFRCT 

83: PROCEDURALCDELW-M) 

84: OELW'O.O 

85: IFCM.LE.O.) DELM-1.0 

86: END 

97: WMZ-W»RSH2CC 

88: freo*alogio:whz) 

89: W2»M*M 

90: M4*U2»H2 

91: 

92: COMMENT COMPUTE HM2 AND 0H2 
93: HW2*1.0' , CTRU2»N2+1.0) 

94: HDU2*UC4*W2/( (M2-Sffl2)*»2+4.0»2D2*ND2»H2) 

95 : DM2 •! X«I Z»M4-W02»<IX»I Z*I Y2*3 . 0«I Y»I Z-3 . 0»I Z2 ) «W2 

36: Dw2*DW2-»4.0*N04««IY-IXJ4CIY-IZ> 

37: DU2-0M2/SCAFAC 

98: DM2-DW2»DM2 

99: 

100: COMMENT COMPUTE F'S 

101 : FTHH2*HM2«SCFAC2x(tIY«M2»3.0*N02«aZ-IX)>*«2} 

102: FRM2-UZ»M2-W02*UY-IX>>m2«HW2xDW2 

103: FYW2*(IX»H2»4.0«W02«(IZ-IY1)*»2»HN2/DW2 

104: FRyH2*WC2»IXYZ2»M2»MH2xDH2 

105: 

106 .'COMMENT COMPUTE PSD'S 
107: 

108: COMMENT PITCH PSD COMPUTATIONS 
109: PPW1-FTHW2«IY24ThOD2 

110: PPW2-FTHU2»IY2«TH02«W2 

in: ppw3*ftf442*oy 

112: PPWR*SIGB02»DELH4(PPW1 ♦PPW2+PPW3)/SCFAC2 

113: PPWD*PPU«HDM2 

114: 

115: COMMENT ROLL PSD COMPUTATIONS 
116- PRH1*FRM2«IX2»PHI0D2 

117: PRW2»FRW2»IX2»PHI02»H2 

118: PRM3-FRM2«HG2»IXYZ2«PSI02 

119: PRW4-FRM2»0X 

120*. PRW3-FRYW2»IZ2«PSI0D2 

121: PRH6-FRYH2«IZ2»PSI0Z»W2 

122: PRH7-FRYW2»MO2«IXYZ2»PHI0Z 

123: PRM8-FRYM2«GZ 

124 : RS*R*SICB02«DELM 

125: PRH*RS*CPRW1+PRU24PRM3*PRV»4+PRH3+PRW6*PRM7+PRM8)^SCFfiC2 

126: PRWD-PRM»HDW2 

127: 

128: COMMENT YAM PSD COMPUTATIONS 
129: PYM1 •FYH2«I22«PSI0D2 

130: PYM2-FYH2»IZ2«PSI02»H2 
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131S PYW3*FYW2»W02»IXY22»PHI02 

132: PYW4-FYW2*OZ 

133: PYWS*FRYW2«IX2*PMI0D2 

13C FYW8-FRYW2»IX2«PHI02«W2 

135: PYW7-FRYW2»W02«IXYZ2»PSI02 

133: PYWS-FRYW2*0X 

137: PYU-RS+1PYW1 +P V W2+PYW3+PYW4+PYW5+PYW6+PY147+PYW81/SCFAC2 

138: PYWD*PYW»HDW2 

13S: 

140: 

141:C0MMENT PREPARE VARIABLES FOR OUTPUT 

142: RSIG*R+SIGB02«DELU 

143: RSIGDB-10.0»PlOG101RSIS>+FACTOR 

144: PPW3DB- 1 0. OvALOGI 0(PPU3 ) -SCF20B+F ACTOR 

145: PPWDB*10.0*ALOG10(PPM)4FACTOR 

146: PRW4DB-10. 0*AL0G1CCPRM4 1-SCF2DB+F ACTOR 

147: PRU8DB- 1 0 . 0» AL0G1 0 ( PRW8 } - SCF2DB+F ACTOR 

148: PRWD8-l0.0*AL0G10tPRW>+FACT0R 

149: PYW4 DB-10.0»BLOG10( PYW4 > -SCF2DB+FACT0R 

150: PYU8DB- 1 0. 0«ALOG1 0CPYH8 ) -SCF2DB+F ACTOR 

151 : PYUDB* 1 0. 04RLOG1 0(PYW J+FHCTOR 

152: 

153: PPWDDB -10. 0»RLOG1 0 ( PPWD > +F ACTOR 

154: PRUDDB- 1 0 . 0* AL0G1 0 ( PRWD ) *F ACTOP 

155: P YWDOB - 1 0 . 0» ALO G 1 0 « P YWO > +F ACTOR 

156: 

157:C0MMENT SAVE NUMBERS IN FILE 

158: CALL FSAVE CFREO, PPWDB, PRWDB, PYUDB, FNUfl) 

159: 

1E0: 

181: TERMTtT.GE.TFINAL) 

162: DERIVATIVE 

163: algorithm ialg-3 

164: 0Q-INTEG<1. 0,0.0) 

165: END 
166:END 
1671TERMINAL 
168: 

169: COMMENT CLOSE FILE 
170: CALL FCLOSEIFNUM) 

171 :END 
172: END 
EOF: 172 





ELT WED-1CV20^82-10:S0M2-<C, ) 

12337*IS<1 >.SRHBxFSUBSCO> 

1: •FOR, IS FF.FOPCN 
2: SUBROUTINE FOPENtNl 

3: BEHIND N 

4: RETURN 

3: END 

e:«for,is ff.fsrve 

7; SUBROUTINE FSRVECk,P1,P2,P?,N) 

a: NRITEIN,100>H,P1,P2,P3 

9; 100 rORHRT(4G14.8> 

to: RETURN 

11: END 

12:«F0R,IS FF.FCLQSC 
13: SUBROUTINE F CLOSE (N) 

14: encfile n 

13: RETURN 

16: END 

EOF: 16 
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APPENDIX E 




GEOMETRIC ERROR ANALYSIS PROGRAM LISTINGS 


100 

D 

000552 

DOLtA 

0000 

R 

000263 

DPITCH 

OOOO 

R 

000156 

DPSA 

OOOO 

R 

000067 

C 

>00 

R 

000064 

DPSI2 

oooo 

R 

000244 

droll 

oooo 

D 

002262 

DSXY 

OOOO 

R 

000170 

l 

>00 

R 

000125 

DIET A i 

0000 

R 

000130 

DTETA2 

oooo 

R 

000213 

dtmeta 

OOOO 

D 

003214 

t 

>00 

D 

001240 

OV 

0000 

R 

000225 

OVAN 

oooo 

D 

02267a 

ERR 

oooo 

D 

023102 

1 

>00 

D 

013626 

FI 

oooo 

D 

000352 

Fll 

oooo 

D 

006360 

FI2 

oooo 

R 

000000 

f 

>00 

I 

030711 

J 

oooo 

I 

03q706 

K 

oooo 

D 

000302 

land* 

oooo 

0 

000626 

l 

>00 

D 

001S3o 

PDXV 

oooo 

D 

013640 

PDV 

oooo 

R 

030710 

PI 

oooo 

D 

017166 

I 

>00 

D 

004252 

PSI 

oooo 

D 

001776 

PSXY 

oooo 

D 

003264 

PXT 

oooo 

D 

0007*2 

( 

>00 

D 

004266 

RA 

oooo 

D 

006520 

Rho 

oooo 

D 

010254 

RH01 

0900 

0 

010436 * 

>00 

D 

026522 

R01 

oooo 

0 

002522 

8 

oooo 

R 

0000*7 

SDH 

oooo 

D 

013350 

« 

>00 

0 

013166 

8 MFT 

oooo 

D 

022532 

SLUN0A 

oooo 

D 

015576 

SOURT 

oooo 

D 

010620 


>00 

D 

027614 

SOUT 1 

oooo 

D 

002236 

SXY 

oooo 

D 

0032*0 

T 

0009 

D 

022502 


>00 

R 

000052 

T0H2 

OOOO 

0 

007314 

tetai 

oooo 

D 

007322 

TETA2 

oooo 

D 

014732 


>00 

D 

007526 

TII2 

OOOO 

D 

007330 

Tit 

oooo 

D 

007336 

712 

0900 

D 

00314* 


>00 

0 

015552 

T 1 

oooo 

D 

025372 

Til 

oooo 

D 

005130 

T2 

oooo 

D 

023206 


>00 

0 

022750 

VP 

oooo 

0 

022734 

VPX 

oooo 

D 

027734 

VS 

OOOO 

0 022556 


>00 

D 

014756 

XI 

oooo 

R 

000075 

X0 

oooo 

0 

000326 

XY 

OOOO 

R 

000133 


>00 

D 

01207* 

X2 

oooo 

0*0*1530 

*22 

oooo 

D 

017224 

X3 

oooo 

D 

000674 


>00 

0 

014136 

YI 

oooo 

D 

005062 

YO 

oooo 

D 

0231«0 

YO 

oooo 

D 

0113*6 


>00 

0 

011712 

Y 22 

oooo 

0 

020316 

Y3 










I* 

2* 

3* 

a* 

5* 

a* 

7* 

a* 

a* 

10 * 

ii* 

12* 

13* 

14* 

15* 

16* 

17* 

ia* 

JO* 

20 * 

21 * 

22* 

23* 

24* 

25* 

26* 

27* 

26* 

20* 

30* 

31* 

32* 

13* 


C THJ5 PROGRAM CO m PUTES The ERROR SENSITIVITIES AnD GEOMETRIC 
C ERRORS OF GROUND POINTS INDUCED BT ALTITUDE*R0LL»YAN AND 
C PITCH ERRORS AS NELL AS EARTH ROTATION AND SHUTTLE MOTION 


PARAMETER NLAN0A*10tNL0MDA*lO*NDMaS»NFIaSfNP81*3»NTETA*3 
S«NDTa3*NDETAaS«NOF!as«NDPSAalOtNDTHEat0tNCOMalS(NSOHS2 
REAL H#R tDH(NDH) «OLAmDA(nl*mOA) *DLDHOA(nlohDA) .|OH(nIoh) • 
STOMl,TOH2»DFIlCNFI5»DPI2(NPl)iOPSll(N*Sn»OPSI2(NPSl), 
SDPSI(NPSI) tOBETACNPJI) tXO( NLOHDA) tOFI(NCPl) t 
SOTETa t (NTETA)*0TET*2(NTETA) • XO(NLOmDA) •OPSA(NDPsA) 
S.DT<NOT).DTAU.OETA(NOETA).OaLFA<NC*SA) iOTHETA(NdTHE) 

S t DTA n(NCOH) f OROLL(NCOH) .OPITCH(NCO m) 

DOUBLE PRECISION L*MOA tNLAMOA) .XY(NLA mDA) ,FI t (NFl) «FI2(NFI) 
SiOXY (NDH.NLAMDA) , lOmd*CNLOMd*} .Y(NlOMba) .PY tNDFt. NLOHDA) 
S.DYlNOFI. NLOMOA) »PDXY(NDH,NLAMO*J.PSXr(NLAMOA«NoH). 

»8XV (NLAHDA) »08XY (NLAMOi*NOH) »8 (NL»hOA,NOH) »BY1 (nF I* NLOMOA) ( 
STXYKNLAHDA) .T*Y2(NU*MOA),0T*TtNLAHD*l •T(NLAMDA)tPXY(NDH«N|.AHD 
I A) •By2(NF1«nL0MDA) iOBYINF I t NLOMOA) . B( nF I.nLOmOA) .P tl (NPSI) » 
S8ETA(NP81 ) «RA (NLAMDA) (DELTA (NPSI iNLAMHA) «D(NPS1 .NLAMOA) • 
SPSA(NOPSA) « ALFA(NDPSA) •OOLT*(NOPSA,NLAMO*} (YOINtOMDA) « 
ST2fNOTME.NtOMDA),DEL(NDTHE»NLOMOA),RHO(NOTNEtNLOMOA) 
»«TET»j(NTETA),TETA2CNTETA),tI1(NTETA),TI2(NTETA), 
STIU(NTETA«NLOMDA)«TlI2(NTETAfNLOHOA)tDELl(NTETA>NLOHDA)t 
SDEL2CNTET A.NLOMDA) t RHO 1 ( NTET A • NLOHDA) « RMQ2(NTET A. NLQMDA) • 
SSOURTl (NTETAtNLOMDA) .S0URT2{NTETA,NL0-0»1 *X11 (NtETA,n L 0 MDA> , 
»Y1 ltNTfTA, NLOHDA). X22{NTETA,NU0HDAJ»Y22<NTETA,N L 0H0A) 

DOUBLE PRECISION X2(NC0M,NL0MDA), 
lSHFT(NTETA.NL0M0AJ.SEN(NTETA*NL0M0A}.YAM{NC0M)*R0LLtNCOM), 
SFItNDFI).POY(NDFI,NLOMOAJ .YKnOThE. NLOMOA). THETA InDTHE). 
»Xl(MOTME. NLOMOA) ,T1 (nOTHE) .SOURT{N oTHe,NLO“DA) . 

SSHFlNOTHE. NLOHDA). PITCH(NCOH). 
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ORIGINAL PAGE IS 
OF POOR QUALITY 


sa* 

H< 

J** 

ST* 

s«* 

jw 

oo* 

«l* 

•2* 

as* 

at* 

*5* 

80 * 

87 * 

as* 

as* 

so* 

SI* 

52* 

5S* 

Sa* 

55 * 

50* 

57* 

58* 

5** 

so* 

*»• 

*2* 

03 * 

•a* 

05* 

00 * 

or* 

os* 

00 * 

70 * 

71* 

72* 

7S* 

7a* 

75* 

70 * 

77* 

78* 

70* 

80* 

81* 

82* 

8S* 

sa* 

85* 

86 * 

87* 

88 * 

80* 

* 0 * 


IIIHCCM'NLOMOA) •Yl(KCOMtNlOMOA) 

tt DX88 (NCOflt NLONDA) tTAUtONCOAf ETAfNOf 78) • 

8V(NCETA)»tLUflDA(NDCTA)«Cl.UNOA(NOfTA)tVX(NOITA)*vr(NOfTA)i 

tDEX(NDCTA.NOT) tDCY(NOETAtMDT)tEM<*MTAtNDT)f V8.V8l(N0EVA)t 

»VRCNDETA).ANGLChDETA>tDANCl(NDCTA).DCl*tNDeTAtNoT) 

l*3ELY(WDET*.NDT)*E8*0(MDETa,t(DT>»T0(NL0x0A) 

S(T33(NC0««N|.0hd*) fDEl (NCOHtNLOHOA) «Tll (ICON) t 
•72 (NCOM« NLOHOAl • 801 (NCOMt NIONDA) ,|QUT| (NCO*f UOmOA) 

DATA COH(K) l K«l»NDN)/l.,.l. t 2«2.,|,,.J,,a.»»a./ 

DATA HtR«OTAU/aOO.t*S5*.785>85./ 

DATA (Dri(K),KBl«NDrn/18.tJ*.*20.«tl.f22./ 

DATA (DP8A(K) t Kot«»D»8A)/.5'l,«!.St2'*2.5t3.t3.5«*<t8.S*S./ 
DATA (DTHETA(K)tX B l t NDtHE)/20.3t21.5*22.5*23.S«:a.5 f as.<ea.» 
•aS,i«a.*a7./ 

data CDVA*(K)'K*:«NCO«)/O.»0*«0.tl.t>l.*0.f8.»0.t|..*|.«0" 
•O.tO.tl.t-1./ 

DATA (DRXTCrf(K)tKBlfNCOK)/0.tl. ••!•*•, t0.t22.St22.5tI8.1«>l.* 
122.5. aS.taSitAS.taS. .as./ 

DATA (OROU.(K)*K*ltNCON)/20,«20.t20.tt0.i20.tOtfl.t-l.tOtfO.» 
•O.t 1. t-l.f O.tO./ 

data (DL* N DA(X)«!*l«NlAMOA)/0.tl.t2 ( tS,ta, t S»fO ( i7. ( 8.««./ 
DATA (0|O'*IIA(!liIll«NlO*OA)MM*8 l l*7«i>liiilMca(«*)i(>] f i 
l>l»te.fl.«2.'3.>a.tS.tO. ( 7.«8.tO./ 

DATA f|DHCl)tX*l.N8DH)/,01.,02/ 

DATA T0KltTDH2/.01t-.01/ 

DATA (Dril(K)tKBltN8J)t<O8I2(K) l IUHtN8I)/.Olt.02t.OSt-.Oit 
l-.02t-.0S/ 

DATA (DPIXl(K)tKaltNMI).(DMX2(«)tKattNPSX)/.Oit.02t.05t 
l-.01ta.02i-.05/ 

OATA (DTETAl (K) tKalt NTETA) t (0TETA2(K)<K*t t4TETA)/.0] t .02* 
I.OSt-.Ol t-.02t-.05/ 

DATA (XO(K)tKaltNLOMDA)/l«*o./ 

DATA (XO(K) f K8ltNLOMDA)/t«*o./ 

DATA (DT(K)tMl»NDT)/|.t.lt.01/ 

data (OETA(K) t**l tODETA)/0»t20.*aD.t80».80»/ 

81*1.1815*265 


DO 100 KBltNDH 
00 90 XaltNLAMD* 

IAMDA(T)*DLAM0A(1) *81/180. 

*7(I)aR*C { I6H/RJ *DC08(LA*DA(X> )-0808T{l-(l*H/8) ttltH/RJ • 
<D8I9(LAK0A(Xn«D8X9(LA-PA(l))n*0SIN(LA-0*(Xn 
8XY(K«l)aR*(( It (H.OHfk; j/R)i0C08 (L**9Al I J )-DIORT t l-( l.lMtDH t * 
S>)/«)*{l*(H*0H(KJ)/«)*03X*lCLAH0ACn)*08TH(LA*DA(X))))»08Xl» 
SCtAHOA(XJ) 

DXV(N.Z)*8XY(KtX)-XV(X) 

X8(DXY(RtX).EQ.O.) SO TO 80 
RDxY(K.n*DAB8(DxY(K.X>)*l00./DABS(XY(X)) 

60 TO 90 
60 8DXT(KtX>*0. 

90 CONTINUE 
100 CONTINUE 


DO 102 JaitNtDH 
00 101 laltNLAMOA 
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91* 

92* 

93* 

999 

95* 

96* 

97* 

9«« 
99* 
00 * 
01 * 
02 * 
03* 
0«* 
05* 
06* 
07* 
00 * 
09* 
10* 
11* 
12* 
13 * 
U* 
15* 
10* 
17* 
10 * 
19 * 
20* 
21* 
22* 
23 * 
2«« 
25* 
26* 
27 * 
20* 
29 * 
30* 
31 * 
32 * 
33 * 
3 «* 
33 * 
3 »« 
, 37 * 
130 * 
1 39* 
1*0* 
IM* 
1*2* 
103* 
:*o* 
103* 
1*0* 
1 * 7 * 


l.OMOAa)aOl.A’‘DO<Il**!/10O. 

SXY(I>*R*{{l+M/R)»OCOS(l.AMDA(in-O0URT(1*(l*H/R)*(Wh/R)* 
1OS1N{l*mO*(I)1*D0IM(l*MO*Cj)))*O31H(L*RI>*(IJJ 
PSXY(XfJ)aR*( { 1* tH*0OH( J) l /R)*0C03(LA*3A(I) )«O0ORT( 1*( 1* J) j 
»/Ri*(i*t“*soH(jn/«)*ot:M(i.*9o»(:i)*O0tN(uHo*{in))*osiMc 
tlA*OA(X)) 

DS5f(IfJ)*PS*TtI.J)-3*VCll 

3(l»JJ«O0Xr(I.J)/SOHCJ) 

101 CONTINUE 
1q2 CONTINUE 


00 103 I*ltN|.A"Oi 
UA«0A(1)*0LAHDA(I)»PI/160. 

TX71 (!)*■•( (l*(H»TOH|)/q)*DCOI(L*HOA(tn-O0ORT(t>(U(H*?DMl )/«)• 
S(t«(M*TOHl)/R)*O0IN(uAMOA(Z))*O0ZN(l.AMD*fI)}))*D0XN(LAMDA(J)) 
TXV2(I)»B*(tl*(H*TOH2)/«l*OCO0(tAPOA(:n«OIQR7<l»(U(H*TOH2)/RJ* 
}(l*:H«TOH2)/R)*OSZR(L* M OA(Z))*O|lN(UAHO*(Xm)*O0ZH(LAHOA(Z)) 
0TxYtX)*TXYl(X)-TX»2(IJ 
T(Z)ooTXY<X)/(7DH1»TOM25 
103 CONTINUE 


00 in KaltNDPI 
ri(K)a0FI(K)«PX/t09. 

00 110 J*1 *N',ONOA 
*.OHOA(jj«OLONOA(J)*Ri/iao. 

Y ( J)a»R*C (1 *m/R)«0C03(L0n 0A(J) 1*0 SORT ( I* t 1*H/R)a ( 1 *h/R)* 
J03ZN(t0H0A(J))*0IZN(L0H0A{J)))l*0SlN(|.0N0A(J)) 
PY(K«J)»-R*C(l*H/R)*OCOS(LOHOA(J)*rXt*n«0$aRT(l»(UH/R)*(UM/»)* 

SOSZN(LONOA(J)*Fl(K))*DSXN(UOMOAU}*FX(K)n)*OSlN(l.OMOA(J)*PZ(K)) 

OY(XfJJ«PY(X*JI-Y(JJ 
XP(V(J).£Q.O.) 60 TO 109 
POYCKt J)*OA0S(OV (Kt J) )*100«/0A9S(Y{Jl 1 
60 TO 110 

109 POV(K«J)*10.0100 

110 CONTINUE 

111 CONTINUE 


00 130 k*1«n*I 
PZ1{*)*0PX1C'}*PI/1S0. 

FX2<*>*0FX2(K>»PX/1®0. 

00 120 JatiNLOPOA 
LONOA(J)aOLONDA(J)*Pl/l0O. 

8Yl(K,Jl*-R«(Cl*H/R)*OCOJ(lO*OA(J}*Fl!(KJ)-O0aRTU«<UM/R)* 

S<l«H/R)*OSlN(LONOA(J)*PZt(K))*OtZN(LONOA(J)*PXt(a)))>* 

S0SXN(L0H0A(J)«PX1(X)) 

BY2(A«J)a.R*((l*N/R)«DC0|(L0P0A{J)«Fl2(K))-0SQRT(l*(t*H/R)« 

S(1«H/r}*OSZN(LOHOAUi«PZ2(X))*O0ZN(LOnOA(J)*PX2(K))))* 

10SINtL0H0A(J)*ri2{K)) 

S0Y(K,J)*0Y1(K.J)-BY2(X.J) 

B(KtJ)aSBY(X«J)/(0PZl(Kl«0FZ2(xn 
120 CONTINUE 
130 CONTINUE 
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- ** 


148* 

ta«* 

no* 

HI* 

152* 

193* 

13«* 

t5S* 

198* 

197* 

138* 

139* 

1 * 0 * 

1 * 1 * 

1*2* 

1*3* 

1*4* 

1*3* 

1**« 

1*7* 

1 * 8 * 

1*4* 

170* 

171* 

172* 

173* 

174* 

173* 

176* 

177* 

178* 

177* 

180* 

181* 

182* 

183* 

188 * 

189* 

186 * 

187* 

t88* 

187* 

1*0* 

1 * 1 * 

192* 

193* 

19** 

193* 

19** 

197* 

148* 

1*4* 

200 * 

201* 

202* 

203* 

20** 


00 132 KaitNOPSA 

DAi.rA(K)«U80.>0*tA(K))/2. 

P8ACKJa0P8A<«)*PX/180» 

ALFA(K)«0ALPA(K)*PI/l80, 

DO 131 Jal«ML*»08 
LA**0*( J)»DLA M OAf J)**I/180. 

RACJlNa*({l*M/R)*0CQS(CA*0A{J))«»C2QRTf|»Cl*H/R)all*M/R) 

t*oaxN (lahoa <2) ) *osiNa ahoa ( jn ) i *D8IM(land* tj) ) 

OOlT*(K.J)a(*A<J}«0|!N(P|A(Kn>/OIXN(lirA(K)] 

13t CONTINUE 
132 CONTINUE 


00 13a Kal'NPfX 
0PSX(K)a0P8Xl^K)-09Sl2(<) 

D8ETA(H)a(lB0,a0PS!(Kn/2. 

PSI(K)aOPSI(K)*PX/IS«« 

BETA(«J«O8CTA<KJ*PX/l80, 

00 139 JaliNLAMOA 
LAHDA(J)aOLA"OA{J)«Pt/ie8. 

BA( J)»R*( ( l*N/B)*OC08(t*POA(JJ J»08Q8T(1»C t*N/8)a(l*M/9)a 
I081 N((.AMOA(J);*OSXN(LAnOA(J])))*OSIN(LAMOA( 2)) 
0euTA(K,J}a(RA(J)*0SXN(P8X(K)n/08XN(B8TA(xn 
OCit J)»DtLTA(«* JJ/OPBl (*1 
139 CONTINUE 
13* CONTINUE 


00 29* KaltNOTMC 
THETA(K)aOTHETA(K)4PI/tBO. 

90 203 jattNLOMOA 
UCHOA(J)aOLC M DA( J) *PI/180> 

r0(J)a-R*ai*H/R)*0C0t(l0M0A(an*08M7(>*(l*N/«)*(UH/R)*0fIN 

s(L0N0A(jn*0iiN(L040A(j)m*08XN(ipH0*(jn 

T1 (KJ«N*OTAn(TNETA(KJ J • 

T2(K * JJb»OA 0B (M/0C0J(THETA(K) 1 )*OTANitONOA( JJ J 
DEL '<• J)a08QRT(Tl (K)«Tl (K]*72(KtJ)*T2(Kt J)) 
RHO(*t«J)aOATAN(0EL(K«J)/N) 

8euRr(AtJ)*R*((l«H/R)*0C0t(RH0(K,2))«0SaRT(l S (UH/R)*(l9N/p)*0BlN 
3CRH0(Kf jn* 0 *IM(RHO{K»J)l) J*08INCRH0(K. J) ) 
XI(K ( J} a (S0URT(K>J)*Tl(K))/0CL(K«2) 

Yl(K« jji(80uRT(R» J)*T2<*. J) j/OEKXf J) 

SHP (X • J) *OSQRT (XI (KiJ)*XI(K«J)*(YI(K?J)«VO(J))*(YI(X>J)aYO(J])) 
203 CONTINUE 
20* CONTINUE 


00 20B X* 1 iNTCT A 
TETA1(K)*0Y£TA1 <*!•*!/ ISO. 
T£TA2(n)«OTETA2fA)*PI/lBO. 

00 207 jaitNLOMOA 
LONOA(J)*OLOHOACJ1*PI/1BO. 

TIl(X)«H*OTAN<T£TAl<*» 

TI2(X)*H*DTAN(TCTA2{x) ) 

TIIl(K«J)a-0AB8(M/DC08(TCTAl(R)))«0TAN(L0N0A(J)) 

TII2(K«Jla-OAB8(H/DCOS(TKTA2(Km*OTAN(LONOA(j>) 

DELnK«J)a080RT(TIl(R)*TIt{R)*TIIl(KtJ)«TIll(K>jn 
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.ORiGWAL PAGE IS 
QE POOR QUALITY 


2*9 

20 * 

207 

208 
20 « 
210 
211 
21 2 
211 
21 « 
215 
21 * 

217 

218 
210 
220 
221 
222 
221 
22 * 
225 
220 

227 

228 
220 
210 
211 
212 
211 
21 * 
235 
210 

217 

218 
210 
2*0 
2*1 
2*2 
2*1 
2 «« 
?*5 
2*0 
2*7 
2*8 
2*0 

250 

251 

252 
251 
25* 
255 
250 

257 

258 
250 
200 
201 


DE12(k»J)bD3QRT(7X2(K)*7I2(K)oT 1X2(K« J)*TXI2(K»,|)) 
RH01(KtJ]a0AT*N(DELl{K«J)/H) 

RH02(«.J)a0A7AN(0CL2(K.J]/H) 

S0URTl(K« J)mR*({l*M/R)»0COB(SH0l(K.Jl)»n$ORTCl»(l*M/R)* 
»U*H/R)*D3IN(RHai<*»J))»08MfRH01CKfJ)n)*08IM(RHOt(**J)) 
SQURT2(K.J)afl*tCl«H/R}*DC08(RH02(K ( J)l«O10RT(l*(l*H/*)* 

*( l*H/R)«D3lN(RM02(Rt J))*DSlN(RH02(Kt J1)))*D8IN(RHQ2(*«<I)) 
Xll(K,J)R(80URTl(K»J)*TIl(R))/0ELl(KtJ) 
Vll(K,J)a(SOURTl(RtJ)*Tm(i(tjn/DeLUK*J) 
x22(R.J)B(SauRT2fK»Jl«7l?(R))/0Et2(K*Jl 
722C*tJ)o(80URT2(K»j)*7li2tiC*4J)/0EL2(R»JJ 
3 MOTU.J)* 030 RTC{Xllf*.J}-X 22 (K,J»*(xIIU,J)-X 22 (*,J)>* 
Mril(KtJ)-v22(R*J))«(rtt(K*4).r22(Kijni 
SEN(K,J)R8HP7(K«J)/(oTETM(K)-0TET*2(xn 

207 CONTINUE 

208 CONTINUE 


00 2tt Kal'NCON 
V *R (KlaOV«N(K) * 01 / 180 . 

R0LL(K)a0R0LL(K)*PI/l80. 

0ITCN(K}*0OITCH(X)*RI/180. 

00 210 JaltNLOHOA 
U0NDA(J]a0L0 H 0A(J)*OI/180. 

YOCJ)e-R*(( t*M/RJ*C,C05CLOHOA(J))«OSORTtl»{t*M/R)*(l*N/RJ* 

!O3INCLQHOA(jn*08IN(LOH0A(J))))*D3tN<LO«DA(Jl) 

Til (K )aH*OT AM (PITCH (K ) ) 

T33(K,J)a-OAB3{H/DC03(OITCH(Xin*OTANrLONOA(J)**OUL(R)J 
0E1 (Kt JJaOSQRT CTll(K)*Tlt (K}*T31(K.J)*T11(K t 4)1 
R01(K«J)a0ATAN(0EltX«J) /H) 

SOUTt (X« J)aR*((l*M/R;*0C03(R0l(K.J) J»N30RT(I»(l*H/Rj*tl*H/R) 

1*08 IN(R01 (K.4))*OSlN(ROt(K«4))))*OIIN(R01(Kf4)) 

X2(K«J)a(30UTl(K.J)*Tli(K})/0El(KfJ) 

y2(K.j)a(80U7i (KtJ)*T3liiC(j))/0Ct (Ri45 

X3(K.J)aX2(K«J)*DC03(VAH*K})-V2(K(3)*0SIN(yAR(K)) 

Yl(K.J)ax2(K.J)*03lN(VA.(K))*T2(X,j)«aC08(YAa(X): 

OlOFlKt J)a03QRT((X3(*. J)aX0(J))*(X3(XiJ1«X0>JlI:( v l(Kf J)*Y0(J))O 
S(YKKtJ)-VOCJ))) 

210 CONTINUE 

211 CONTINUE 


TAuaOTAU*Pl/180. 

0HECAS2. *PI/8o*00. 

DO 220 K~i'NOETA 

ET A (K) ■ JETA C*5 *01/180. 

V (K} *C m ECA*R*OCOI(ETa(*) ) 
SUINDA(K}aOC03(TAU)/OC03(CTA(Kn 
CUJNOA(X)a030RT(l-«UUNOA(Kl*8LUHOA(*n 
YX{K)aY(K ) *3LU M 0A(K) 

VY(K1«V(I<)*CLUH0ACKJ 

00 215 Jat'NDT 
0CX(K.J)avx(K)*0T(J) 

OEy (K l 4)aVY(X)*0T( J) 

ERR(RfJiaOSQRT(DCX(KtJ)*DCX(XtJ)«OEY(KtJ}*OEY(K,J)) 
215 CONTINUE 
220 CONTINUE 
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262 * 

263* 

26** 

265* 

266* 

267* 

2*a* 

26 «« 

270* 

271* 

272* 

273* 

276* 

273* 

276* 

277* 

279* 

274* 

280* 

281* 

292* 

293* 

286* 

289* 

286* 

287* 

288* 

284* 

240* 

241* 

242* 

243* 

24«* 

245* 

246* 

247* 

249* 

244* 

300* 

301* 

302* 

303* 

308* 

305* 

306* 

307* 

308* 

304* 

310* 

311* 

312* 

313* 

31** 

315* 

316* 

317* 

318* 


V9*2.*Pl*R/59*9.*8« 

T AU«DT AU*P 1/188* 

0»CS*a2.*PI/86*00. 

00 222 KaltNDETA 
ETa(K)>DETA(K)*P1/180. 

VC*) bO m CCA*R* 0C09(ET* cm ) 
9LU N O*CK)aOCO8(T*U)/DC0|(CT*(K)) 

CUUH 0 *OO* 0 SaRTCt- 8 L 0 M 0 *CK)* 8 LUM 0 *CKn 

vx(*)aVCm*SLUNDACN) 

vv(K}aWCK}*CLU"0*C*) 

VPX(i)aVX(K)*V8 

VP(msOSQRT(VPX(K)*VPX(K)*WCK)*VV(K)) 

ANCL(*)sOATAN(VPX(K)/VV(K)) 

OANCL(K)«*NCL(X)*180./PZ 
CO 221 Jali NOT 
oeuxca.j)avP*:m*oTcjj 
OELVd, Jlavv(K)*OTC J) 

ERR0(i«J)*D3QRT(0ELX(K«J)*0ELX(k, J)*OEL v (K,J)*OELV(k< J) ) 

221 CONTINUE 

222 CONTINUE 


DO 237 «ll»NOM 

«RITEC6»236) R»N,DH(k) , (OLANOAC 1) • jal , nla« 0 A) t CxT ( I ) , I* 1 «nla h O 
**) • (PxT(K»X) flat tNtAaOA) ( (OxT(K«l) ( lai «NtAaOA) •(POxv(KiI)«Xalf 
SNLANOA) 

236 FORHATCtHt* '•»•••••••*•••••••••• GEOMETRIC ERROR XNOUCEO BY'* 

S' ALTXTUOE ERROR a*******************!,//////////, 

Jl EARTH RA0XU8 » a» ,F9 .«,ikm' *///«! NOMINAL ALTITUDE 1 , 

9' HB»,F*,0» * KN »,///, • ALTXTUOE ERROR OH *I,FI.«,'KMI» 
*/////////»» VlENi ,1x, ' ANCLE *«/• • OLAROA (DECREES) '» 
S12X,F2.0,9(9X,F2.0),///,I NOMINAL LOCATION (KM) t f 100> l.«* 
>///•• SHIFTED LOCATION (IN) • , 1001 t.*,///»7x» 'ERROR* ,6*. 

S '(KH), ,10011. «»///,**» 'ABS(£RROR) I, 3X, ’ IS) i, 2X, rS.3, 
S8x«F*.3i8(6XiF9>3) ) 

237 CONTINUE 


00 234 Jat.NSOM 

*RITE(6,23S) RtH,SDH(«*) t (OLANOA (11,1*1 *nlamOA).( 8XV(1) f la 
S1 ,NLAhOA),(P 8XY(I, J},Isl,NLANDA) t (DSXV(If J),:*I,NLA mDa),(S(I,J) 
tfIal,NLAMOA) 

238 FORHATCIHI, '•*••»*•*••»«*•••*•*» GEOMETRIC ERROR SENSITIVITY* 

S' *ITH RESPECT TO AltITUOE ERROR M******************, ,//////// 
S//,' EAR’h RADIUS » a*,F4.A,'IMi,///,i nominal ALTITUOE H 
SM.OtUM' ,///, I ALTITUOE ERROR SOm a* ,F«,2«'*M' ,/////////, • VIE*' 
S*5X. I ANGLE' ,/* ' OLAmdA(OEGREES) ' , t2X,'2.0«9C9X«F2.0) .///, 

S' NOMINAL LOCATION t |H) * , 10 q 1 1 •*,/// • • SHIFTEO LOCATION (IN)', 
S100 11. *•///, 7X< 'ERROR ' ,*X,'(*M)i,ioOtl.A,//.6X, i ERROR ' ,/• 
si SENSITIVITY (IN/*N) I, 1001 1 •*) 

23* CONTINUE 


«RITE(6,240) R,H,TOHi,TOM2 t( OLAHOA(I),Ial,NLAMOA)*(TXVt(I),ia 
1 1 , NLAmOA) «(TXY 2(I), Xat fNLAMQ A) , (DTxYCI) « Tal, NLA mQA] « (TCI), la 
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ORIGINAL rAG'i .3 
OE FCOa QUALITY 



Sl.NLAMOA) 

240 FORMAT! IM l , i ••**•••• •****••••••• GEOMETRIC ERROR SENSITIVITY* 

S' mITH RESPECT TO ALTlTUOE ERROR 

$//////////» • EARTH RAOIUS R R'tFO, 4, "»"»///,• NOMINAL' 

S' 4LTITU0E M «',F4.o,'KM>,///,l ALTITUDE ERROR 1 TOhI 
SF4.2 ( ikm',/// ( < ALTITUDE ERROR 2 T0H2 «F5.2, 'km, , 

!////////» 1 VIE"' .SX. « Angle * »/t ' DLAmDa(DEGREES) • > 12x< 
$F2.0,9(9X,F2.o),///,' SHIFTED LOCATION' ./, I OuE TO ERROR' 

S' 1 (F**) • ■ 1301 • .4,///, I SHIFTED LOCATION, »/» » OUE TO' 

if ERROR 2 (hN)«. 10011,4. ////.' DIFFERENCE (*"),, 

JIOHII. •«///. «*.'ER»00'|/.' SENSITIVITY (KH/nHj |,| 9 

SOI ! .4) 


DO 24b KBltNOFI 

■ RITE(6.2«S) R.H.OFI(X), (OLOMDA( J) , J»l , 10) , (Y( J) , J»1 . 10) . 
S(PY(K. J).jBi,tO)«(0Y(K.J),Jal.i4),(PDY{*«J)tJal,10),(DL0H0A(J). 
Sj41l.NL0MDA),(Y(J),J.Sl,NL0'«0A),tPY(K,J),J«ll,N L CMDA),{0T(K,J), 
S Joi 1 «NL0"0*) »(P0 Y(K»J)»Jb 11 f NLOHQA) 

243 F0RHATUHIt'»**»44*»*»»*4*»««A** GEOMETRIC ERROr INOUCED By*. 

S' ROLL ERROR •*••*•»•***•*»»•»•** i ,///////, 

S' EARTH RAOIUS R ■'»F9,4,ikm»»///»I NOMINAL *LTlTUOE h a*. 
SF4,0»'KM',///,' ROLL ERROR OFI « • ,F4.n, ' OESREE ',////// , 

S' VlEM'»5Xt 'ANGLE'*/.' OLOM0A(DEGREESt'.12X,F3.0,9(ftX,F3,0). 
»///*• nCInal LOCATION (KMjifiODtt .a*///* » ShjFtEO * f 
S' LOCATION («MJ lODt t. a*///. 7X. 'ERROR'. bx. ' (KN) ' .toOU.Ai 
S///.4X. 'A0S( ERROR) ' »5X» ' (X) ' .FtO.a.SFl 1 .4, JX .F5.2, 

S//////.' VIE*' »3X» ' ANGLE ' » 0L0M0AC0EGREES)'.jJX.F2.0. 
S8(9X,f2.0) • ///. ' nominal LOCATION (KM) ' ,901 1 ,4.///, « SHIFTED' 

S' LOCATION (KM)',90ll.4,///,7X,'ERR0R'.»X.'(KM)I.R0ll.4, 
S///.4X. • A8S (ERROR) ',3*,'(X)'.FtQ,4«8F1|.4) 

24b CONTINUE 


00 300 KMl.NFI 

MRITEfb.250) R.H,0Fll(K),0FI2(K),'(0L0»0i I) » Jbl . 10) , (SY1 (K 
S» J).J«t.lO).OY2(K.J).J«l, 10).COBY(N«J).. If 10) t (B(K.J)tj4 
SI . to) f (OLOMOA(J) t JR 1 1 1 NLQMDA) , (BY 1 (K, J) . J»1 1 1 NLOMOA) « (BY2 
S(K.J)fjRllfNLONDA)t(OBY(KtJ).JRlltNLOMOA)f (B(Kfj)tJ B ll.NLOM 
IDA) 

230 FORMAT ( 1HJ « • ••*4444»***»*»*A»b»4 GEOMETRIC ERROR SEnSITIVI' 
S'TY -ITH RESPECT TO ROLL ERROR • »- *Mb44*«a**a*«a»b « ,///// 
St* EARTH RAOIUS R »'.F9,«,i*Mt,///,l WOMINAL ALTITUDE M 4', 
SFa.O.'KM, ,///,■ ROLL ERROR l DFIl ■' tF3.2f 'DEGREE '»///» 

S' ROLL ERROR 2 0FI2 ■' .F5. 2. 'OEGREE './////, ' VlE«'.3X t 
S' ANGLE'./ 1* OLOHOA< DEGREES) 12X.F3.0.9(BX.F3.0) •//, 
s* Shifted location'. /.* cue to error i (kmji.iooh.r 
$//.' shifted location',/, i due to errmr 2 (kh)I, 

SIODtl.4 ,///,' DIFFERENCE (KM)', 10011. 4.//, bX. 

s 'E rror*,/, 1 sensiti yity (km/deo ',iodii. a,/////. 

S' YiEt". 3X. 'ANGLE', /,' DLOMOA(OEGREES) ',12 XiF3.0,B(SX, 
SP3.0) ,//, ' SHIFTEO LOCATION', /•' DUE TO ERROR 1 (kM)', 

S9D| 1 . d, //« 1 SHIFTEO LOCATION',/,' CUE TO ERROR 2 
SRDll.a,///, I DIFFERENCE (KM) , • «n 1 1 .a,//,»x, 

S 'ERROR',/, • SENSITIVI TV (KH/OEG) ',9011.4) 

300 CONTINUE 
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WA l.c«a.3.\r* 


378* 


1T«* 

380* 

m« 

382* 

383* 

»M 

383* 

186* 

387* 

388* 

38*9 

3*0* 

3*1* 

3*2* 

3*3* 

3*«* 

3*3* 
3*6* 
3*7* 
3*8* 
3**« 
*00* 
*01* 
*03* 
80S* 
It)* 
*09* 
•08* 
•07* 
408* 
•98* 
• 10* 
«tl* 
412* 
*13* 
at«* 
4t3* 
Mil 
417* 
418* 
41*4 
*20* 
*21* 
422* 
423* 
424* 
429* 
428* 
427* 
428* 
•2** 
•30* 
431* 
432* 


00 32t K«1*N0PSA 

nRITE (8*320) R.M,DP8A(K),(DL»HDA(»,J*1.NLANDA), 
S(D0LTA(K»3)»3*1 *NLAmdA) 

320 FORHatUNJ, '••••••••••«••••••••• GEOMETRIC ERROR INDUCED' 

S' BY yll E»R0R ••••«••••••••••*•••••«////////(• EARTH • 

S ■ RADIUS fl »i,F9.4*'KH',///*i NOMINAL altitude h • ••M.O* 
S'KN',///,' TAR ERROR 0R8A •» *73.1. 'DECREE* */.'//////, 

»• VlEN'.SX* • ANCLE ••/«• OLAMOA(OCCRZES)'*I2X*F2.9,*(«X,P2.0)* 
S.////.7X. 'ERROR '.OX* i «H) • *10011.4) 

321 CONTINUE 


DO 400 xeliNPSI 

«RITC(**330) RiHfO*SIl(K) »CR8I2 (K) • (DL* m OA(1)«ImI«NL*MQA) • 

S (DELTA (K tj) «J* I INLAND A) • (0(K*3l I J*l *Nt A»IJA J 
390 FORMATtl"!, ••*•••*••••**»*•*•••• 6E0RETRIC ERROR SEnSITIVITV' 

*• »ITh RESPECT TO YAn ERROR •••••••••••••••••••••*////////« 

S' EARTH RAOluB R «'tM.4.l«|l,///,l NOMINAL ALTlRuOE H 
SF4.0«1KMt,///,i TAN ERROR > DPSlt *P9.2» 'OECREt'*///* 

S' YAm ERROR 3 0P0I2 •' «P9*2»' DECREE* •////////* • YIEm'«9I* 

S' ANCLE' t/t' DLAMOA(OCCREES) ' , 12X«F2.8*8(*X*F2«0) »////, 

S' TOTAL SH1PT (*H) ' «10ol 1 «4*///*8X* 'ERROR* ,/, 

S' SCN3IT!V!TY(KM/OCC)MOOtl.«) 

400 CONTINUE 


DO 43 t «*t»N0TME 

■RITE(8»*30t R*M*OTMCTA(K)*(OLOMOA(3)*3*t«tO)*(xO(3)'3«l*l8)# 
S(VO(3),2Mt>lO).(Xt(K.J),3*l*I0)*(VI(K.31'3*lilR)«(SHF(K«3)*3*l'10) 
8. (DLOMOA (3) • I : (NLOMOA) . (XOC J) «3Rt t iNLONOA) t (YO(3) * Jal I • 

SNLONOA) • (XI (Xf 3) f 3«lt • NLOMOA) » (TI(K*3l «3R1 1 'NLOHDA) » (8MB (K«3J *3B 

SlItNLOMOA) 

•30 FORMAT ( IHI » ' a******************* GEOMETRIC ERROR INOUCED BY' 

S' PITCH ERROR •••••••••••*••••••••'*/////, • EARTH RaOIUB R •' 

S*P*.4.'KN',///,I NOMINAL ALTITUOE rf a'*P4,0«'KMl*///,l FITCH' 

S' ERROR OTHETA a '*F9.1' 'DECREE ••/////»• VU*'t9x« • ANCLE •*/» 

S' OLOiOA (DECREES) * t 1 2X.F3.0**(BX*F3,0) *///• ' NOMINAL X COORO' 

S'. (KM) *«7S ,F2.0.*(*X,F2.0) «///* ' NOMINAL Y COORO. (KM)'* 

SlOOlt. •*///*' SHIFTED X COORO. (KM) •* t OOtt.4*///* ' 8«IFTE0' 

S' Y COORO. (KM) ' * tOOl t .4 1 ///*?X* 'ERROR' »4X« ' (KM) '.1 OOtt.4* 

S/////*i VICMi*9X«'AN6LC'«/t' OLOmda(CICREES) '*12X.F3.0» 

S8(8X.F3.0) «///•' NOMINAL X COORO. (KMj ' *7x.F2.0.8(8y. F2.0) * 

S///*l NOMINAL Y COORD. (KM) l*4Dlt,4(///*' SHIFTeO X COORO, • 

S' (KN)',*Olt. •*///♦' SHIFTED Y COORD. (KM) • ,*DI 1 .4. ///. 71* 

S 'ERROR »*8X* '(KM) *,*011.4) 

491 CONTINUE 


90 900 Kmi.ntET* 

mR:TE( 8*480) R*H*OTETAl(K}*oTCTA2(R)*(DL0MDA(3)*3Rl.tO)« 
S(Xll(K«3)f3Rl*10)i(Y)l(Ki3)*3Rl*t0)*(x22(Kt3)*3Bltl0). 
S(Y22(k.3)*3r1«19)*(ShFT(X*3] »3*1.10)».SEn(k,3).J*1.|0j. 
S(0l0M0A(3)*3Rtl'NL0M0A)*(Xll(«*3)i3«11'NL0M0A)* 
S(Yll(K*3)f3Mll«NLOMDA)*(X22(K*3)*3*U*MLOMOA), 

S(T22(K. 3)*3*1 1*ML0 M DA) • (|MFT(K*3)*3*1 1 * n LOMDA) * 
S(SCN(K«3) ,3*1 1 •NLQMOA) 
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433 

434 
alS 

436 

437 

438 
434 

440 

441 

442 
403 

444 

9 45 

446 

447 

448 

449 

450 

451 

452 
45J 

454 

455 

456 

457 

458 

450 

460 

461 

462 

463 

464 

465 

466 

467 

468 
466 

470 

471 

472 

473 

474 

475 

476 

477 

478 
474 

480 

481 

482 

483 

484 
465 

486 

487 
468 
464 


460 706*47 { INI* »*••*•*»»«»**»»****»» GEOMETRIC ERROR SENSITIVE 
S'TY WITH RESPECT TO PITCH ERROR **M44M«4«M60m*»,//// # 
S' E4RTH RADIUS R • • *64.4* 1 4N • 1 //* • NOMINAL ALTITUDE H a'* 
SFa.o* »kmi ,//, 1 PITCH ERROR 1 DTETAi *75.2* 'DEGREE'*//* 

S' PITCH ERROR 2 DTETA2 •• *75.2* 'DEGREE' *////* • VlER'tSXt 
S' ANGLE' •/• ' DLOHOA(DC6 REES)'i12X*P3.0*4(8x*P3.0)i//« 

S' SHI7TE0 X COORD.'./*' DUE TO ERROR 1 (KM) • MODI). 4*//* 

S' SHIFTED V COORD. './*• OUE TO ERROR 1 (KM) 1*10011. 4.//. 

S' SHIFTED X COORD.'*/*' DUE TO ERROR 2 (KN) • *10011.4*//* 

*' SHIFTED T COORD.'*/.' .OUE TO ERROR 2 (KM) « * t 001 1 .4.///. 

»• DIFFERENCE C*H)i,iooil.4*//«6*« 'ERROR' •/» • SENSI' 

S'TIVIty(km/uEG) '. 10011. 4.////*' VIEM'.5K* 'ANGLE'*/* 

S' 0L0H0A(DEGREE3)'*12X*F3.0*6(8X.F3.0)*//* 

S' SHIFTED X COORD. ••/«' OUE TO ERROR 1 (KM) » *4011.4,//* 

S' SHIFTED v COORD.'*/*' OUE TO ERROR 1 (KH) ' *9011.4,//, 

S' SHIFTED x COORD.',/,' OUE TO ERROR 2 (KM) I ,9011 .4,//* 

S' SHIFTED 7 COORO.'*/,' OUE TO E»RO« 2 (KM) I *9011 .4*///* 

S' DIFFERENCE (*M) 1 *901 1 .4,// * 6X * 'ERROR »*/, ' SENSI' 

S'TIVITT(km/0EG) ',4011.4) 

500 CONTINUE 


00 601 K6l*NC0M 

mHITE( 6*600) R, H, DTA*. (*), DRITCH(X), DROLL (K)*(DL0HDA(J),J41, 10), 
S(X0(J)*jal«lO)*(V0(J)*Jat •10)*(X3(K*J)*jMl*10)*(V3(K*J)*j4t,10}* 
S(OIFF(K,J),j a i,10)«(OLOHDA(J)*Jatl,NLOMOA),(xO(J)*Jall*NL9MCA) 
S.(VO(J),J»ll.NLOHDA),(X3(K.J),J«ll,NLOMOA),(»3(K*J)*j4ll*NLOMDA) 
S*(OIFF(K*J)*Jal 1 *NLOmDA) 

600 FORMAT! 1 * 1 , 'aaMaaaaaaaaaaaaaaaa geOmfTRIc ERROR INDUCED SV 
S' yam, FITCH AND ROLL ERRORS •6«a*#6»A44*a*a4*6»a' *////, 

S' EARTH RADIUS R a ' *F4.4 * 'km'*///,' NOMINAL ALTITUOC H a', 

SF4.0* 'KM' •///« ' YAM ERROR OYAM ai *F4* t •• DEGREE ' •///* 

S' FITCH ERROR OPITCH at * F4. 1 * ' OEGREE '•///* ' ROLL ERROR' 

S' OROLL a* «F4.1* 'DEGREE' *////*' VIE"' ,5X, ' ANGLE'*/*' DLOMDA' 

S' (OEGREES) '*12X*F3.0*4(8x*F3*0) *///•' NOMINAL X COORO, (KM) I 
S*7X*F2.0*6(9X*F2.0) «///*• NOMINAL V COORO. (KM) • , |0D1 1 .4* 

S///*' SHIFTED X COORD. (KM) • * 10D 1 1.4*///, l SHIFTED Y COORO.' 

S' (KM) * *1001 1.4*///*7X* 'ERROR' *6X* • (|») ' * 1001 1.4*////* 

S' VIEW *5X* • ANGLE' */• ' OLOMOA(DEGREES)',12X*F3,o*8(6X,F3*0)* 
S///.I NOMINAL X COORO. (KM) l*7X*F2.0*A(4X*F2.0)*///* 

S' NOMINAL Y COORO. («M) '*4011.4*///*' SHIFTED X COORO. (KM) I 
S.4011,4*///*' SHIFTED Y COORD. (KM ) 1 *401 1 ,4*///*7X* 'ERROR' • 

S6X* ' (KM) • *4D1 1 .4) 

601 CONTINUE 


00 7C5 Kal.NOETA 

"RITE (6, TOO) R,0TAU*0ETA(K) ,W(H) . (oT(J).Jal • NOT) * (DEX(K* J) * 
SJal «NDT) * (OEY(K«J) *Jal*NOT)*(ERR(K*J)*Jal*NDT) 

700 FORMAT! 1*1* <a*a««a»»***» •***•*»« GEOMETRIC ERROR INDUCEO BY' 
S' EARTH ROTATION aaaa6aa«««a«««***6«*' */////////* ' EARTH* 

S' RADIUS B'*F4,4* 'KM 1 *////,' ORBIT INCLINATION a'.Fa.O* 
S'OEGREEl «////«• LATITUDE a' *F4.0* 1 DEGREE • *////////* 

S' GROUND POINT SHIFT VELOCITY OUE TO EARTH ROTATION a'* 

SO 1 0.4* • KM/SEC ' •////» 

S' TIME INTERVAL (SEC) I»3F15.2,///*' SHIFT IN',/, 

S' X DIRECTION (kM)',3X,3015.4*///*' SHIFT IN*,/* 
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490* 

491* 

492 * 

493* 

4949 

4999 

496* 

4979 

496* 

4999 

500* 
5019 
502* 
50J* 
504* 
505* 
504* 
507* 
SOM 
5099 
510* 
511* 
512* 
51 J* 
514* 
515* 
514* 


»' v DIRECTION (*M) », 3X.3015. *.////.' TOTAL IHIPT' 
*' (KH) 1 *31*3015, *3 
709 CONTINUE 


00 707 KaltNOCTA 

-RITE (4*7061 R*DTAU»DETA(K),V(K).VS*VK«),DANGL(K)»(0T{J)* 
t J*l (NOT) • (OELX(K* JJ * JBt .NOT) • (OELY(KtJ) • J*1 .NOT) * (ERROCK* J) « 
»J*1»N0T) 

706 pornatcihi* ■ ••*•••*«••**•*•«•••• ground point imipt induced' 

!• M OHUTTLE MOTION AND EARTH ROTATION •*••••••••*•••••••*•• 

If//////////*' EARTH RAOIUS a'*F9. 4 *'KMi. ////,' ORBIT 1 
I' INCLINATION 9'*F*.0* 'DEGREE' *////»• LATITUDE s'*F*.Q* 
••DECREE' .////*• GROUND POINT IHIPT VELOCITY OUE To EARTH* 

• » ROTATION a' *0I0«4* IKM/iEC •////»• SHUTTLE GROUNO |PEEO *•* 
SOlO.Af'KM/tEC'f////*' ground POINT IHIPT VELOCITY RELATIVE' 
!• TO NAOIR POINT a< *010.4* 

| ' RH/IEC • *////*' ANGLE OP VELOCITY VECTOR nITH Y.Avli ,i, 
1010.4* 'DEGREE '*////* ' T>t INTERVAL (IEC) ' »3FlS.2*///« 

*' IHIPT IN'*/*' X DIRECTION (RN) • *3x* 3015.4*///* 

I* IHIPT IN'*/*' Y DIRECTION (KM) I *3X *5015.4*////* 

I* TOTAL IHIPT (kH) • *3X«3D15.4) 

707 CONTINUE 


•TOP 

CO 


•1 CTPl.134 IUP|t 10*719 

«A ihutexp.image/abi 

2IJ-JB 04/27/13 10131121 

2 SUPb. 734 CPU*. 090 109,302 CC»CRa.*33 
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cr;ov,l Fs/iZ ’ o 

OF POOR QUAUTT 


6IC"tT*IC t**b* 1'IOUUB ■» illlTittl t*P« •••••••••••••«••••♦• 


’MTN "tOIUt O •IHt.UHl 

KJHINH H.TITUO* 
llTJTUOt 1**00 ON ■ J.«K 


»«■ wcu 

nwimmui 


t. 

t. 

2. 

*• 

*0«JNAi. LOCATION 

(IN) 

• AMI 

.«t*2**«1 

*110T*««2 


IhXOTCO IK*TI0« 

(*«) 

•MM 



•iriMii 

CMON 

(KM) 

•lit* 




Alt (CN*ON) 

(« 

•III 

•Ml 

.20* 

•w 


•*•«*•***•■■•••••*** OlOOfTOlC lOOO* OlX$ITl*f TO »JtH ItlHC* <0 


■• 5. 


7* 

•• 





•S»25«M2 

,*)*0**02 


■ •IIMH! 
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•siinttc 

••HMU 

.mj'lll 

• t*0**M* 

•t22«**00 

•imNim 

.10*0**** 
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.20* 

.2*2 

•291 

.111 


•LttTuDt t»*OI 


(MIH RtOIU* * nHt.HIMK 
xo*iN«t n'ltuoc n ••••.<■ 
■ITITUOC to*0* I TONI • .•ton 
■ltxtuoc uooo i too* » •••!«" 


VIC* ANtU 

OtAMOA(KSNCCt) 


o. 

1. 

2. 

9« 

«• 

9* 

A* 

7. 

«• 


Kimi LOCATION 

oue to t*oo* i 

(OKI 

• MAI 


• iumn 


.2T««*«*2 

,»<|.H2 

.«2*****2 

•■010*002 

.00*0*0*2 

•*i«i**> 

out TO 10*00 2 

(00) 

• •Ml 

.***2***1 

.IITMII 

.2AAMAA2 

•27«T*«o2 

• S)M««IS 



•9*2*«4«a 


01M(K«CI 

(AN) 

• MM 

,nn»us 


• I0aa«M2 

• 1 Joo**o 2 

.ITOt-ol* 

• 2IM«M2 

iZlSIMM 

•2tt«*««2 

■ttii*ib 

ill AON 

HNltTlVlfY (AN/ftM) 

• •Ml 

# |7AA«MI 

• )«*2.MI 



• •79)«MI 

.1IU.IM 

.1220**00 

.1** i«*** 

• I9IAMI 
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Wi’OrtAL PAG! *c 
OF POOR quaut} 


•••••••••••••••«•«•• t(ot(«ic t»«o« t«o<x(D ** mu. i""’* •••«••*••••••••••••* 


(ONTO •ioim * mm. »•**«■• 

000100), OtTITuOt « 0777.00 
«IH.W 7*707 on it.ttt'U 


Mi* •*«.« 

DVOMPA(PCft>(II) 


•«a 

••• 

•7. 

••• 

*9. 


•5a 

•2 a 

•l. a. 

N0*lN*L LOCATION 

(HU) 

•CSMTlM 



.19*****2 

.277***72 

•2*7***«2 *1177***2 

• ittMV 

loimo L0C4TI0O (OO) 

•MltCll 


•0277.772 .19*0. **2 

•2777*7*2 


.|)*7**o 2 .*7*2***l 

.*•*• ..**7*,7*| 

Utoa 

Coo) 

..nikiii 




..7717*7*1 




Ct) 

II.1MT 

12.0*71 

■•••17) 

. I*.T*7* 

2*. *7*2 

!9.tM« 

SM7M 

9*. 71*7 

••••• 

Vila *N&V« 

Ot0«U*(O«6UCt) 


w 

2a 

J. 

•• 

9. 

•• 

Ta 

• a 

*. 

oo"I«*l lOCOTlOo (0*1 

•••HllNI 

..!)* 7. ,02 

.•llTttlM ..277 0 *o«2 

..197**7*2 

•••2*7*772 



»»imo loC»tio* fooi 

.#11*7. *10 

..iitmim 

..27***772 ..S99*«**2 




••7***«0*2 

(MO* 

U"> 


a»t 



..7771**01 

••7ll****l ..7193*001 


•tSUMOft) 

Ci) 


90.0719 

)l.o)Tl 

29. <2*1 

7».«oT9 

l*.«l*9 


12.7*91 

11.100* 


IIMUHMHMmin CCOMCTOIC l»0> 7folITi»|T* >I(H «im(* tu .XL t"*'** »«•••••••••••••••«»• 

'Ml* iidiui « nlta.ltiio 

iooiooi oirtruot " mii.u 
IOLL i«*o* i one • .IINMU 
'Oil t»«0« 2 or n a ••UOfSaCI 


film 4NGU 

)t0»0A(uC6»UI) 
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••* 

*9. 

-A. 
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•2. 

•la 

0. 

INIITTCO LOCATION 
>U( TO (M0« 1 

COO) 

.7111*772 

.9717*772 


# «14'aO02 

• JMMK 

.27*7*072 

720IMMI 

iUHHOI 


• aMHlam 
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COO) 

•7)07*7)2 

,9712*772 
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,o21)»7*2 

,)9g7»((2 

,27*9*772 

.21)1*7)2 


.7)12*071 
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•toll*)** 

••1027*777 
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.*|0|)*07* 
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••1****7«* 
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1. 

2. 
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5. 

•. 

7. 

•. 
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out TO (UM I UK) ..7*92*771 .,t<M<Ni ..llihlll ..•Tlltlll ..MIKni 

I.ITTIO LOCtUOo 

Out TO 1*000 i C**l '.MIMII ..IkllMI ..lUTtlU ..1HMM1 •,)*7l**»2 ..•iM'tlC ..«»07 Mil 

oifnn»cc c«»i '.utiMM ••too**)** ..■•••*••« ..to****)* ..i*ti**«* •.iot««*** ••lom.o. 

»CO*lTX*IT»C*o/0t*l •■•HmIII ..IKCMI ..TIImHI ..T7o**77l «W>m<NI -.7W7»»o| ..7l))’0*t -.TIT1..Q1 
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•MMflftMlHMMM i|0*tTIK IMgfc IMPS tt$ fT PITCH 


ORJGff^AL FA*£ iZ 

OF POOR QUALITY 
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GROUND PATTERN AND IMAGING 
DISTORTIONS GENERATION PROGRAM LISTINGS 
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1 


00105 

ie« 


S <61EML1NES*NLINES)«YG1(NL1NES«NLJNES).XP6CNLINES) » 

00105 

19. 


STP61NLlNCS).X6BfMLlNES.XLINCS».Y0B(XLlNES.NLXNES)« 

001C5 

20* 


SAFLAG ( NL 1NES. ALI NISI. XP7INLINES) «YP7t KLINES) 

00106 

21* 



00106 

22* 


BATA RtPI.F.0.TSTAR«0T.VS,H<00.0S/63S6.7C5.3. 141592654. 

00106 

23* 


S100000B. »400»»-5»5..0C4*7.l3S9,440»»«0039».9/ 

00121 

24* 


BATA EPS1.BET.YSTAP.YSTAP1.PB. PS. FLAG/.02.1*. 

00121 

25* 


1-5., -4 34. 7.. 0055. .65. 2/ 

00131 

26* 


BATA YSTA Y. YAWP. TAWS/ -5. • .0037. .925/ 

00135 

27* 


BATA BROLL.VSTAO.RB.RS/0..-S.O.. 0042.1./ 

00142 

26* 


BATA BTAU.PETA.X STA. Y STA. YS4/8S..40.. 0 .. -5..7.202/ 

ceiso 

29* 



C01S0 

30* 

C 

INITIALIZE THE PLOTTER 

00150 

31* 


CALL P6NPLT 

00151 

32* 



00151 

33« 

C 

CENTER PLOT 

00151 

34* 


XCENTSC.O 

00152 

35« 


YCENTsl.O 

00153 

36* 


call ORI CINE YCENT.YCENT1 

00154 

37. 



CO 154 

36* 

c 

BEfINE THE TORN OF THE PLOT 

00154 

39* 


PTYPEs*L INLIN* 

00 155 

40. 


(LEAsp.O 

C0156 

41* 


YLENsp.O 

00157 

42* 


CALL PLFO CRIP TYPE. RLE N.YLENI 

00160 

43* 



00160 

44* 

c 

SCALE THE PLOT 

00160 

45. 


«Xll>:-a 

C0161 

46* 


>Il2i:l 

50162 

47. 


YY«l>r-fl 

00163 

4f • 


YYi2»*e 

C0164 

49* 


NXX=2 

00165 

50* 


NYY*2 

00166 

51* 


NGES-l 

C0167 

52* 


N6YS- 1 

C0170 

53* 


CALL PLSCALEXY.NXX.NGX.YV.NYY.NGV) 

00171 

54. 



0C171 

55* 

c 

TO START PLOTTING 

00171 

56* 


CALL PLGRAF 

'7172 

57. 



00172 

5a. 

c 

BEFINE GRIB LIMITS 

C0172 

5°* 


XLEFTs-5 

00173 

60* 


XRIGHTsS 

00174 

61* 


XOELs f XRICHT-XLEFTl/ENLINES-11 

0C175 

62* 


YBOT* -5 

OC 176 

63* 


YT0P=5 

00177 

64* 


YOEL* • f T OP-Y 80T1 /ENLINES-ll 

00200 

65. 



00200 

66* 

c 

CONSTRUCT BATA 

00200 

67* 


YCsYTOP 

00201 

63* 


BO 15B INORsl.NL IRES 

00204 

69* 


(CsXLEFT 

30205 

70* 


BO 100 IPTsl.NLINES 

S0210 

71* 


XNONsR.SIMXC/RI 

00210 

72* 

c 

XBstXC-XNCRSI.F 

00210 

73* 

c 

ENONsxC-xB 

00211 

74. 


Y NO Ns VC 



- < 4 


'I 
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0212 

75« 

X 1 1 HO R« IPTIsXNON 

0213 

76* 

YtlHOt.lPTIsVMM 

02X4 

77. 

<C3«C.v»CL 

02i: 

70. 

100 CONTI HUE 

0217 

70. 

YCsYC-YOCL 

0220 

80* 

150 CONTINUE 

0222 

61* 

IMFLXG.EO.il GO TO 200 

0224 

e2« 

1MFLXG.F6.2I 60 TO 208 

0226 

83* 

III FL A6.EU.3I 60 TO 231 

0230 

84. ' 

1FIFLX6.FB.4) 60 TO 251 

0232 

85* 

bO TO 261 

0233 

86* 


0233 

87. 

2g0 CONTINUE 

0234 

ae« 

9C 292 Isl.NLlNES 

0237 

87* 

00 201 Jsl.NLINES 

0242 

90* 

NFLAG ( 1. J 1=0 

0243 

91* 

201 CONTINUE 

0245 

92* 

202 CONTINUE 

0247 

93* 

YUsVSTAR 

0250 

94. 

00 206 Lsl.NOY 

0253 

95* 

OLANOAsDS 

0254 

96* 

00 205 KsltNLAPOA 

0257 

97. 

LAHORsOL AHOA.PI/180. 

0260 

98* 

XH3*R.|(1.h/R>.COSILAN»AI*$ORYI1*I1«H/RI 

0260 

99. 

tSIXIL AND A J.SINILA80AI II.SINILANBAl 

0261 

108* 

00 SO Isl.NLXNES 

0264 

101* 

DO 70 0=1*NL1NES 

0267 

102* 

<Z=XH-K< l.JI 

0270 

103* 

YZsVW-YT I.JI 

0271 

104. 

IR=1 

0272 

105* 

JR=J 

0273 

106* 

IM A8 SfXZI.LE.EPSI.AND. ABSIYZl.LE.EPS 11 

0275 

107. 

70 CONTINUE 

0277 

109. 

80 CONTINUE 

0301 

107* 

GO TO 99 

0302 

110* 

204 NFLAudR.JRlsl 

0303 

111. 

XFI1R. JNIsOX-TANILAHOAII 

0304 

112* 

YFIIR.JS ISYST AR*(L*1 ) - .0263 

3305 

113. 

WRITE I 6.9S1 tb.TW.lR, JR.DLAPDA.XFI1R.JRI 

C316 

114. 

98 FORXA I ( / « 2F1 0.5. 2 14.3 FI 0.51 

0217 

112* 

49 BLAND AsD LAHOA*OD 

0323 

116 * 

205 CONTINUE 

0322 

117. 

yw=tstar.l.ot*vs 

0323 

no* 

20E CONTINUE 

0325 

11°. 

60 TO 299 

0326 

12C* 


3326 

121* 

208 CONTINUE 

0327 

122* 

00 2 Isl.NLXNES 

C332 

123* 

DO 1 Jsl.NLINES 

3325 

124. 

NFLAG 11. J IsO 

0326 

125. 

1 CONTINUE 

0*46 

126* 

2 CONTINUE 

0342 

127. 

OTHETASC. 

33*3 

12*. 

VA3YSYAP 

3344 

12*. 

DO 217 Lsl.NDT 

03*7 

130* 

THETAsPYHCTA.Pl/lFO. 

0350 

131* 

OLPHOAsPS 


% - - Z* 

* 4 
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03S1 

132* 

DO 216 K=1«MLPH0A 

10354 

133* 

LPH0Ax0LPP0A*Pl/l*Q. 

>0355 

134* 

TlsO* TAN (THETA) 

0356 

135* 

T2=-AbS( t /COS (THETA) ) *T4M (LPNDA) 

0337 

136* ’ 

0EL=SQ8T(T1*T1*T2*T2> 

*0360 

137* 

RHOxAT AN (DEL/D) 

’0361 

13** 

S0U8T =t* ( (1*0/8) • COS (RPO ) -S0RT(1>(1*0/*) • (1*0/8 >*S1N 

0361 

137* 

S(tHO) *SlN(RHO)))*SIN(RH3) 

0362 

140* 

Tls(SOURT*Tl)/DEL 

0^63 

141* 

Xls(SQURT*T2)/0EL 

<0364 

•• 142* 

YBsYI *ta 

0365 

143* 

X6=(I 

0366 

144* 

00 22 Isl.NLlNES 

10371 

145* 

60 11 J:1«NL1NES 

'0374 

146* 

<C=KS -X( I • J ) 

0373 

147* 

TCxTa-T(I.J) 

•0376 

148* 

1631 

’0377 

149* 

J6=J 

'0400 

150* 

IFtABS(XC>.LE.EPSI.ANO.ABS(YC>.LE.EPSl> 60 TO 210 

0402 

131* 

11 CONTINUE 

0404 

132* 

22 CONTINUE 

=0406 

153* 

60 TO 215 

10407 

15** 

210 NFLA6 (I6« J6>3) 

’0410 

155* 

X6(IS , J6 > =0* f-TAN(LP“DA > > 

•0411 

156* 

T6( IS, J6»stB0T*(L-1)* .0280 

<0412 

157* 

•/RITE ( 6* 212) XB,YB,16,J6,DLPHDA.XG(I6«16).Y6(16.J6> 

0 423 

138* 

212 FORHAT(/*2F10.3*214«3F10«S) 

0424 

15«* 

215 DLPHDA=DLPHDA-PD 

■0425 

160* 

216 CONTINUE 

<0427 

161* 

Y AsYST AP»L«» 0288 

0430 

162* 

DTHETAsL/2500* 

0431 

163* 

217 CONTINUE 

•0433 

164* 

60 TO 299 

0434 

165* 


0434 

166* 

00 220 1 = 1. NUNES 

•0437 

167* 

DO 214 Jsl.NLlNES 

<0442 

143* 

NFLAG1(I.J)=0 

0443 

167* 

217 CONTINUE 

0445 

170* 

220 CONTINUE 

• C447 

171* 

TEHTAxbTHETA.pl/lfO. 

’0450 

172* 

ETsDET*Pl/l!0* 

C451 

173* 

YAlsYSTAPl 

0452 

174* 

DO 227 L=1.M0T 

3455 

175* 

DLPNOA=PS 

•0456 

176* 

00 226 K=1«NLPH0 A 

C461 

177* 

LPNOAsOLPNOA*PI/160* 

C462 

17** 

T11=0*TAN( THETA* ET) 

0463 

179* 

T21=- ABSC O/COSIT HETA* ET 1 )*TAN(LPH0 A) 

0464 

180* 

D EL 1= SORT (T11*T11*T21*T21> 

0465 

131* 

RHO 1= ATAN (DELI/D > 

0466 

182* 

SBJRT l:R* ((1*0/S )*C0S(RH01)~SQRT(1~(1«0/R) *(1*D/R )• 

C466 

US- 

TS IN (R h01»*S1H(RH01)))*S IN (R HOI) 

0467 

164* 

111=($CURT1*T11)/K£U 

0470 

145* 

<I1=(S6URT1*T21) /6EL1 

C471 

186* 

YB13YI1*YA1 

£47? 

187* 

XB1SX11 



l 




XLJ 


0473 

193* 


00 222 1*1*NLINES 

0476 

189* 


00 221 J*1*NL1NES 

0301 

190* 


(CKtBWdiJS 

0302 

191* 


YClsTBl-tll.J* 

0303 

132* . 


Ifil*l 

0504 

193* 


J£1=J 

0303 

194* 


IFCABSIYC13.LE.EPSI.ANO.ABStTC13.LE.EPSl) 60 70 223 

0307 

193* 

221 

CONTINUE 

osu 

196* 

222 

CONTINUE 

0313 

197.' 


60 TO 223 

0314 

198* 

223 

NFLA61 (1 61* J 61) * 1 

0313 

— 199* 


«61 <1C1. J61)=0*C-TAN(IPN0A)> 

0316 

200* 


Y61 (I 61* J61)*TSTAPl*t L”1 )».0288*Yll 

0317 

201* 


WRITE 16*2241 *81 • YB1 • 161 • J61* OLPHOA* X61I 161* J61) • 

0317 

202* 

SY61(I61*J61) 

3330 

203* 

224 

FORMA T|/* 2 11 0*3*214*3 FI 0*5) 

3331 

204* 

225 

OLPMO A*OtPMOA-PO 

0332 

203* 

226 

CONTINUE 

9334 

206* 


YA1*YSTAP1*L*.02?3 

3333 

207* 

227 

CONTINUE 

3337 

209* 


60 TO SOS 

3340 

209* 



9S*G 

210* 

231 

CONTINUE 

1341 

211* 


00 4 1*1 *NLINES 

3344 

212* 


00 3 JsltNLINES 

3347 

213* 


KFLA6II* JtsO 

3330 

214* 

3 

CONTI NUE 

3332 

213* 

4 

CONTINUE 

J534 

216* 


OYAW* 0* 

1555 

217* 


YDsVSTAY 

3356 

218* 


00 24 0 L*1*N0T 

J561 

219* 


VAW*OYAb*Pl/180. 

1362 

220* 


DLVHO AsVANS 

1363 

221* 


00 238 Ksl.NLYNDA 

1366 

222* 


LYMDA*DLYPDA*P1/1E0* 

367 

223* 


RAs-R • (11*0 /R)*C OS ILY ”0 A )-SQATf 1-1 1*0/8)* 11* D/R>* SIN! LYNDA) 

367 

224* 

**SINt LYNDA) ))*SIN(LYHOA> 

370 

225* 


(JsRA*CCS(YAU> 

371 

226* 


YJs.RA*SlMVAbl 

572 

227* 


<E*XJ 

573 

228* 


YEsYJ ♦ YO 

574 

229* 


00 44 1*1 «NL INES 

577 

230* 


00 33 J*1«NL1NES 

602 

231* 


*T=YE-»1 1*J1 

603 

232* 


YT**E-YI1*J) 

6C4 

233* 


IH*t 

603 

234* 


JN*J 

606 

235* 


Xfl A8$<YT>.LE.EPSI.AND.ABS(YT).LE.EPS1) 60 TO 233 

610 

236* 

33 

CONTINUE 

612 

237* 

44 

CONTINUE 

614 

233* 


60 T3 233 

613 

23«* 

233 

XFLA6 IIN*JH)*1 

616 

243* 


<MIIM.JH)*D*f-TANILYNDA)) 

617 

241* 


YMIlM,J»O*YB0T*IL-l)*.02?B 

620 

242* 


WRITE 16*2341 YE.VE*1N*JN.DLYMDA*RH(IN.WMI*YMI1H*JH) 

631 

243* 

234 

FORMA I |/.2F10*3*2I4*3 F10 *51 

632 

244* 

235 

OLYMO AsOLYMDA-YAnO 


l 
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0633 

16S* 

238 

CONTINUE 

0633 

266* 


Y03V$TAV«l«.02E9 

0636 

267* 


0YAN=3.*S1M IL*PJ>'176.> 

0637 

263* 

260 

CONTINUE 

0661 

" 111* 


60 TO 299 

0662 

230* 



0662 

231* 

261 

CONTINUE 

0663 

252* . 


00 6 Ul. MINES 

0666 

233* 


00 5 Jsl. KLINES 

0CS1 

236* 


LfL 66 1 I* J )s0 

0652 

233* 

3 

CONTINUE 

0636 

236* 

6 

CONTINUE 

0656 

237* 


ROLLsDROLL*Pl/16C. 

0637 

238* 


ASTRO sR* I (1*0 /R) • COS ( ROLL )~ SORT (1~(1*D/R ) • (1*D/R) *S IN (ROLL ) 

0637 

237* 

t.SINf ROLL) >)*SIN (ROLL) 

9660 

260* 


VKsVSTAO 

0661 

261* 


00 250 LsltNOT 

0666 

262* 


OLRNDAxRS 

3665 

263* 


00 265 K3l«NLR90A 

3670 

266* 


LRHOR=OLRPOA*P1/1<0« 

3671 

263* 


0IS*“R*H1*0/R1* COSELRHD A*60LL)*SQRT ( 1~(1*0/R)«(1«D/R 1*S1N 

3671 

266* 

SILRND A*ROLL)*SIN (LRHO A.ROLL) ) )*S1N(LRHDA*R0LL) 

3672 

267* 


AKsOI S*Y$TAO 

3673 

266* 


DO 66 lsl,NLlNES 

3676 

260* 


00 S3 Jxl, MIKES 

3701 

270* 


XLxXK-X( 1*J) 

3702 

271« 


YL=YiC-Y(I.J5 

3703 

272* 


INsX 

3706 

273* 


JH=J 

1703 

276* 


IMA8S((L).LE.EPSI.AN0.6BS(YL).LE.EPS1> 60 TO 262 

3707 

275* 

53 

CONTINUE 

1711 

276* 

66 

CONTINUE 

3713 

277* 


60 TO 266 

1716 

27S* 

262 

LELAKIN, JH):l 

:715 

279* 


j(N( IN, JN| sO*(-TAN(LRHOA) ) 

1716 

280* 


YNI IN «jn > SYS T AO* (L~1)*«02E8 

1717 

261* 


UR1TE (6*263) (K*YK*lP*jn*OLRPDA«XIUln«JN)*Vn(IH*JII) 

730 

292* 

263 

70R6AT (/ *2F 1 0*5*21 6 *3(10 .3) 

.731 

2A3* 

266 

OLRNOAsOLRNOA-RD 

•732 

206* 

263 

CONTINUE 

736 

233* 


YKsYST A0*L*»0289 

733 

296* 

250 

CONTINUE 

737 

2tT 



737 

283* 

251 

continue 

760 

299* 


00 9 in, NUNES 

763 

290* 


00 7 J = l. MINES 

766 

291* 


LELA61C1* J)sO 

767 

292* 

7 

CONTINUE 

751 

293* 

8 

CONTINUE 

733 

296* 


ROLL* 0R0LL*P1/180. 

736 

293* 


OERsO. 

753 

296* 


XSTAO = R*((l*(>/R>.COSCROLL)-SfcRT(l-(l«D/R><(l*0/R>*SIN(ROLL> 

733 

297* 

S.SINf ROLL >1 1 *SIN (kOLL 7 

756 

29** 


YKlsYSTAO 

737 

299* 


DO 260 Lsl*N07 

762 

300* 


ER30ER.PI/190* 


I 


) 



0763 

301* 


DLRHOAsRS 

0764 

302* 


00 235 Ksl.NLRNOA 

0767 

303* 


LANDA=olrpda.pi/1(o. 

0773 

304* 


0 IS Is >8* 1 1 1«0/R 1 *C3S I LR4D A*ROLL*C Rl-SORT (1*1 1*0/81* (1*0/81 

0770 

333* 

S*S1NI LRPDA.R0LL*ERI*S1N( LR8BA.|0LL*E8 1 11 « $ 1NILRN0 A* 80LL*ER 1 

0771 

306* 


<81*0 1S1* XST 40 

0772 

307* 


00 98 1*1.NL1NES 

0773 

303* 


00 77 Jsl.NLlNES 

1000 

309* 


XL1=XK1-VII*J1 

1001 

319* 


YllsYKl-VlI.Jl 

1002 

311* 


181*1 

1003 

312*- 


J*1=J 

1004 

313* 


IFtABStXLll.LE.EPSI.AND.ABS'YLll.LE.EPSll GO TO 232 

1006 

314* 

77 

CONTINUt 

1010 

313* 

83 

CONTINUE 

101*- 

316* 


60 TO 254 

1013 

317* 

232 

LILA61 (1 PI* J811s 1 

1014 

318* 


<811181, JP11*D*«-TANILR8DA11 

1013 

319* 


Y8HI81,J81lsYB3T*tL-ll*.026» 

1016 

320* 


481 TE <6*2331 <K1 *TK1 1 181 • J81 1 0L8H04* X81( 181* J81) * 

1016 

321* 


ST81UP1* JP11 

1027 

322* 

233 

rOR84T</«2riO*5«?14*3riO.Sl 

1030 

323* 

234 

0LR89 AsbLRPOA-RO 

1031 

324* 

233 

CONTINUE 

1033 

323* 


Y81*YSTAC*L*.C238 

1034 

326* 


D£R*0. 2*S INI IL*PI 1/174.1 

1033 

327* 

260 

CONTI NUE 

1037 

322* 


GO TO 299 

1040 

229* 



1040 

330* 

261 

CONTINUE 

1041 

331* 


DO 263 I sl.NLIKE S 

1044 

332* 


DO 262 Jsl.NLlNES 

1047 

333* 


AFLAC I 1* J 1*0 

1030 

334* 

262 

CONTINUE 

1032 

333* 

263 

CONTINUE 

1034 

336* 


T4U*0T4U*ri/lE0. 

toss 

337* 


08EG4 =2. *T 1/(6400. 

1036 

33S* 


ET6sDtT6.Pl/lSO. 

1037 

339* 


V*ON£ G4<R<C0S<ET 4) 

1060 

340* 


SLUNDAsC^SITAUl/CGSIETAl 

1061 

341* 


CLUNO 4SSCRT 1 1-SLUPDA*SLU80A1 

1062 

342* 


VX*V*CLUPDA 

1063 

343* 


VT*V. SLUhOA 

106* 

344* 


VTXs-VX 

1063 

345* 


VTY*VS4-VY 

1066 

346* 


X64SXST4 

1067 

347* 


YAAsTSTA 

1070 

346* 


00 270 Lsl.NOl 

1073 

349* 


0LT8D4*Y4hS 

1074 

330* 


00 269 P = 1«NLYN0A 

1077 

351* 


LVND6sOLtPOA.Pl/ 140. 

1100 

332* 


*AAs. R. 1 ( 1. D/Rl • COSI LY804). SORT 11-11*0/8 }• 11* D/RI*S IN 

1100 

353* 


4ILY80 4 1*S INI LYNDA } 1 1 . SIN ILY80A1 

1101 

334* 


<Sds< 44*4 44 

1102 

353* 


VBdxYAA 

1103 

336* 


00 26S 1*1. MINES 

1106 

337* 


OO 264 JsI.NLINES 


J 

/ 


kil 

L* 


♦ ; 

£ 
! ' 
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« * 


I 


1111 

35s* 



XCCs*BB» »( I* J) 

1112 

359* 



YCC=YBB-Y«I, J> 

:1113 

360* 



100 = 1 

1114 

361* 



JOOsJ 

1115 

362* 



1F<ABS(*CCI.IE.EPSI.AND.ABS(VCC).LC.EPSI> 60 TO 266 

HIT 

363* 


26* 

CONTINUE 

•1121 

364* _ 


265 

CONTINUE 

1123 

365* 



CO 73 268 

112* 

366* 


266 

ATLAS ( 1, J > = 1 

1120 

367* 



*00(100, J 001 s0 *( -tan (LYNO All 

•1126 

366* 



VDO(IOO«JCO)sYSTA*<L>l)*.026E 

1127 

369* 



WRITE (6,2671 »BB,YUB, 103 • JOD.ULYYOA, *00(100. JOOI, 

1127 

373* 


1700(1 DO* J D0J 

11*0 

3T> 


267 

FOR HAT (7, 2110, 5, 214, 3 FI 0,5) 

11*1 

372* 


268 

DLYHD As0LYH0A»D0 

11*2 

373* 


269 

CONTI ‘.UE 

11** 

374* 



AAAsX STA*l*DT* VT * 

11*3 

375* 



TAAS7 STA*L«OT*VTT 

1146 

376* 


270 

CONTINUE 

1150 

377* 




115C 

37e* 

C 


PLOT HORIZONTAL LINES 

1150 

379* 




1150 

230* 


299 

CONTINUE 

1151 

3*1* 



T 1 Cs* ♦ • 

1132 

362* 



NTICs-1 

1153 

383- 



NLINsNLlNES 

113* 

394* 



00 400 IHOR=l,NLINES 

1137 

345* 



DO 390 IPTS1. CLINTS 

116? 

386* 



<P( IPTIsX ( IHGR, IPTI 

1163 

297* 



TP(IPT|S7(IH0P,IPT) 

116* 

3S3* 


300 

CONTINUE 

1166 

359* 



CALL rLCURV( XP, YP.NL1 N,NTIC,T1C> 

1167 

390* 


409 

CONTINUE 

1171 

391* 




1171 

392* 

c 


PLOT VERTICAL LINES 

1171 

293* 



DO 500 IVCkT=l.NLINES 

117* 

3*4* 



CALL TLCUKV 1*11,1 VERT >,VI1«1VERT>,NL!N,NTIC,T1C> 

1173 

3»5* 


SCO 

CONTINUE 

1177 

396s 



IFlTLAC.EC.l) CO TO 655 

1201 

397* 



1F(FL*G.E0,2) GO TO 80S 

1203 

*9»« 



I F t FL AG. CC.3F 60 TO 901 

1200 

JQ«# 



IFKL AG.EC.4I GO TO *09 

1207 

*00* 



uU T3 921 

1210 

401* 




121b 

402* 

c 


PLOT If OR 17CNTAL LINES 

1213 

403* 




1210 

404* 


655 

CONTINUE 

1211 

*05* 



TIC1*0 

1212 

406* 



NTIC1=0 

1213 

407* 



NLIN1 =NL IKES 

121* 

405* 



DO 730 1 HORsl.NL INES 

1217 

40®* 



DO 700 I®Tsl,NLINES 

1222 

410* 



<HlCPT> = *r(IHOR.IPT> 

122* 

411* 



YP1(IPT»=YF(IH0A.1PT> 

122* 

*12* 


700 

continue 

122* 

*13* 



CALL PLCUNV(*P1,TP1,NL1N1.NT1C1,TIC1) 

1227 

41** 


750 

CONTINUE 



L 



i 

t 

< 

i 

i 

i 
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3 . 



1231 

*15* 




1231 

*16* 

C 


PLOT VERTICAL LIMES 

1231 

*17. 



00 30? 1 VERT- l.MLINES 

1234 

*1 6« 



CALL P«.CURV(Xf(l,lVEKT>.vr<l,lVtRT>,NLINl.MTlCl*TICll 

1239 

*19* 


600 

CONTINUE 

1237 

*20* 



00 TO 9999 

12*0 

*21. 




12*0 

*22* 

• c 


PLOT HORIZONTAL LINES 

12*0 

*23* 


605 

CONTINUE 

12*1 

*2*» 



TIC2S0 

12*2 

*25* 



NT1C2*0 

12*3 

*26* 



NLIN2:NLlhES 

12** 

*27. 



OO 550 IMORsl.NLINES 

12*7 

.»28» 



DO 320 IPTsl. NLINES 

1292 

*2*. 



XP2(IPT>sV6(IH0R.IPT> 

1253 

*3C» 



TP2lirT>sTC(INOR.JPTl 

125* 

*31* 


620 

CONTINUE 

1256 

*32* 



CALL PlCURVI VP2* VP2.NL1N2.NT1C2.T1C2I 

1257 

*33* 


650 

CONTINUE 

1261 

*3*. 




1261 

*35* 

c 


PLOT VERTICAL LINES 

1261 

*36* 



00 *51 1 VERT-1. NL1NES 

126* 

*37. 



CALL PLCURVI vGI1«IVERT).YGI1«1VEIT>»NL1N2.NT1C2.T1C2» 

1265 

♦36* 


651 

CONTINUE 

1267 

*33. 



GO TO 9999 

1270 

**0* 




1270 

**!• 

c 


PLOT HORIZONTAL LINES 

1270 

♦*;» 



TIC6=0 

1271 

**3» 



NTIC6=0 

1272 

***. 



NLINGsNLINES 

1273 

♦*;. 



00 *5* IHCRsl.NLINES 

1276 

*46« 



00 653 1PT=1.NLINES 

1301 

4*7. 



<P6 ( X PT) S7G1 IIH0R.1PT) 

1302 

446* 



YP6<IPI>sV61IIH0R.IPT> 

1303 

4*9* 


652 

CONTINUE 

1305 

*50. 



CALL PLCURV CXP6. YP6.NLIN6.NT1C6. TIC6) 

13C6 

*51. 


554 

CONTINUE 

1313 

*52* 




1310 

*53* 

c 


PLOT VERTICAL LINES 

1310 

*5*. 



00 555 IVERT-l.NLINES 

1213 

*55. 



CALL PLCURV|VG1I1.IVLRT>*YG1(1.IVERT).NL1N6.NT1C6*TIC6> 

■m* 

*56* 


655 

CONTINUE 

1316 

*57. 



GO TO «*«' 

1317 

*5‘. 




1317 

*5 a « 

c 


PLOT HORIZONTAL LINES 

1317 

*60. 


901 

CONTINUE 

<132C 

461* 



T1C3=3 

1321 

462* 



NTIC3*0 

1322 

*63* 



NL1M3SNL16ES 

1323 

46*. 



DO 905 1 HCRsl.NL INES 

*1326 

465* 



00 993 1 F T = 1 .NLINES 

13J1 

*66* 



<P3f 1 PTls.H4lHON.lPT7 

1232 

♦ 67* 



YP3HPTISYHT 1H0N.IPT1 

•1333 

*6*. 


903 

CONTINUE 

:i335 

46*. 



CALL FLCUNVI *P3. VP3.NL1M3.MT1C3.T1C3) 

.1326 

♦ 70* 


905 

CONTINUE 


158 


3L 





3WOMA*; FASita" 
° F POOR QUALrr^ 


:is*o 

471* 




31340 

472* 

C 


PLOT VERTICAL LINES 

31340 

473* 



00 907 IVERTsltNLTNES 

11343 

474* 



CALL PLCURVf XH(ltlVERT)tVH(ltIVERT)tNLI93tNTlC3tTIC3> 

->134. 

473* 


907 

CONTINUE 

11346 

476* 



SO TO 94 »9 

31347 

4 77* 




11347 

473* 

c 


PLOT HORIZONTAL LINES 

H347 

479* 

• • 

911 

CONTI N.'E 

11350 

440* 



T1C4S 

•1331 

461* 



NT1C4S-1 

•1352 

492* 



NL1N4 sALINES 

•1353 

483* 



00 915 IHORsltNLINES 

1356 

484* 



00 913 IPTsltNLINLS 

1361 

4«3* 



<P4( 1 PTI S7 NT iHORt 1PTI 

<1362 

486* 



YPNCIPTIsYM IHOR 1 1PT I 

1363 

487* 


913 

CONTINUE 

1365 

498* 



CALL PLCURV(XP4tVP4t*L114t*TlC4tTlC4) 

1366 

469* 


915 

CONTINUE 

1370 

490* 




1370 

491* 

c 


PLOT VERTICAL LINES 

137C 

492« 



DO 91 T IVERTsltNLlNES 

1373 

493* 



CALL PLC URVTXNEl t IVEMTI tVNXlt IVERT ItNLIV 4,NTIC4«T1C4> 

1374 

494* 


917 

CONTINUE 

1376 

495* 




1376 

*96* 

c 


PLOT HORIZONTAL LINES 

1376 

497* 


■09 

CONTINUE 

1377 

499* 



TIC5sfl 

1400 

499* 



NTIC5S0 

142} 

500* 



NLINbsNLlNES 

1402 

SCI* 



00 s l a IHORsltNLINES 

1405 

502* 



00 918 IPTsltNLINLS 

1410 

503* 



<FS(IPT)siHl<lHORtlPT) 

1411 

304* 



YPSlIPTIsTHHIHDRtlPTI 

1412 

505* 


■19 

CONTI HUE 

1414 

306* 



CALL PLCUNV < *PSt YPStNLlNStNTlC5tTlC3) 

1413 

507* 


919 

CONTINUE 

1417 

509* 




1417 

509* 

c 


PLOT VERTICAL LINES 

1417 

510* 



SO 92 C IVERTsltNLlNES 

1422 

511* 



CALL FLCUAVf <H1< 1 1 1VERT1 1 TNII It 1VE ATI t NLlNSt NTIC5tTIC5» 

.423 

512* 


920 

CONTINUE 

1 425 

513* 



SO TO 9993 

1426 

514* 




1426 

515* 

c 


PLOT HORIZONTAL LINES 

.426 

516* 


921 

CONTINUE 

.427 

517* 



TICTsfl 

.430 

518* 



HIIC7SQ 

:431 

319* 



NLIN7SNL1NES 

4 32 

520* 



00 923 IHORsl.NLlNES 

.435 

521* 



DO *2 2 1PT=1. KLINES 

.440 

322* 



<P7I 1 PTI s*0D(lH0Nt IPT 1 

441 

323* 



VP7CI PTIstDDIlHOR.lPT ) 

442 

324* 


922 

CONTINUE 

444 

523* 



CALL PLCURVf TP7t 7P7t NLlN7tNTlC7tTlC7l 

443 

526* 


923 

CONTINUE 
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M 


■i \ 


1447 

527 * 


1447 

52 *» 

C PLOT VERTICAL LIKES 

1447 

S 24 « 

DO 924 IVERTsl.NLINES 

1452 

530 * 

~ CALL PLCUAVf "DO t 1 .IVERTI » VOOf 1 * 1 VERTI«NLIN 7 »NTIC 7 »T IC 7 J 

1453 

531 * 

924 CONTINUE 

1455 

532 * 


1455 

533 * 

C TO FINISH PLOTTING 

1455 

534 * 

"999 CALL LNOPIT 

1456 

535 * 

# 

1456 

536 * 

•• STOP 

1457 

537 * 

ENO 
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